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coupled  transmission  lines  a development  is  provided  to  explicitly  relate  the  four-port,  terminal, 
electrical  performance  of  directional  couplers  to  the  modal  loss  coefficients. 

Losses  are  evaluated  for  examples  of  four  isolated  transmission  lines  and  one  coupled  trans- 
mission line.  For  microstrip  and  coplanar  waveguide  the  computed  loss  coefficients  are  in  reason- 
able agreement  with  experimental  data.  For  inverted  microstrip  and  trapped  inverted  microstrip, 
evaluations  presented  in  both  tables  and  graphs  provide  useful  design  information  for  circuit  ap- 
plications. A comparison  is  made  of  the  total  loss  characteristics  of  microstrip,  coplanar  wave- 
guide, inverted  microstrip,  and  trapped  inverted  microstrip.  The  utility  of  the  analysis  for 
coupled-transmission-line  losses  is  illustrated  for  the  example  of  edge-coupled  microstrip  with  a 
dielectric  overlay  by  comparing  computed  loss  characteristics  with  measured  values.c.The  accuracy 
of  the  loss  evaluations  is  quantitatively  assessed,  and  suggestions  are  made  for  additional  refine- 
ments in  the  solutions.^  — - ; 

The  five  computer  programs  employed  to  evaluate  dissipation  losses  for  microstrip,  coplanar 
waveguide,  inverted  microstrip,  trapped  inverted  microstrip,  and  edge-coupled  microstrip  witn  a 
dielectric  overlay  are  listed.  For  each  computer  program  information  is  provided  as  to  storage  re- 
quirement, execution  time,  compatibility  with  various  commercial  computer  systems,  and  input/ 
output  data  descriptions. 
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COMPUTER-AIDED  ANALYSIS  OF  DISSIPATION  LOSSES  IN  ISOLATED  AND 
COUPLED  TRANSMISSION  LINES  FOR  MICROWAVE  AND  MILLIMETER- 
WAVE  INTEGRATED-CIRCUIT  APPLICATIONS 


INTRODUCTION 

There  is  considerable  interest  in  investigating  and  exploiting  new  transmission  lines 
for  use  in  integrated  circuits  operating  at.  higher  microwave  and  millimeter-wave  frequen- 
cies. This  interest  has  been  spurred  by  the  success  in  reducing  circuit  cost,  size,  and  weight 
through  the  application  of  microstrip  at  lower  to  intermediate  microwave  frequencies.  Un- 
fortunately microstrip  is  discouragingly  lossy  and  more  difficult  to  fabricate  at  higher 
microwave  and  millimeter-wave  frequencies.  These  considerations  have  prompted  the  search 
for  transmission  lines  that  are  amenable  to  integrated-circuit  fabrication  methods  (thin  film 
and  photolithographic  technology)  and  that  have  better  loss  characteristics  than  microstrip. 

To  facilitate  the  investigation  of  transmission  lines  which  offer  potential  for  improve- 
ments over  microstrip  at  the  frequencies  of  interest,  a flexible  computer-aided  analysis  of 
transmission-line  losses  has  been  implemented.  This  analysis  is  suitable  for  application  to 
a wide  variety  of  transmission  lines.  This  report  is  a description  of  the  implementation  of 
that  analysis  as  it  applies  to  both  isolated  and  coupled  transmission  lines,  where  losses  due 
to  both  conductor  and  dielectric  dissipation  are  taken  into  account.  Various  examples  of 
loss  evaluations  using  the  analysis  are  presented  for  both  isolated  and  coupled  transmission 
lines  of  interest.  For  tire  examples  of  isolated  microstrip  and  coplanar  waveguide  [1]  cal- 
culated loss  characteristics  are  compared  to  experimentally  determined  loss  parameters, 
thereby  enabling  an  assessment  of  the  accuracy  of  the  analysis.  These  examples  are  fol- 
lowed by  a presentation  of  computed  loss  characteristics  for  inverted  microstrip  [2]  and 
trapped  inverted  microstrip.  Also  presented  are  evaluations  of  loss  parameters  for  edge- 
coupled  microstrip  with  a dielectric  overlay.  Calculated  and  experimentally  determined 
loss  characteristics  are  compared.  Explanations  are  given  for  the  use  of  computer  programs 
listed  in  Appendixes  A through  E which  are  useful  for  computing  loss  coefficients  due  to 
conductor  and  dielectric  dissipation  for  microstrip,  coplanar  waveguide,  inverted  microstrip, 
trapped  inverted  microstrip,  and  edge-coupled  microstrip  with  a dielectric  overlay. 


FORMULATION  OF  ANALYSIS  FOR 
EVALUATION  OF  LOSSES 

The  following  analysis  of  conductor-loss  and  dielectric-loss  characteristics  for  isolated 
and  coupled  uniform  transmission  lines  is  consistent  with  the  quasi-TEM  models  described 
in  Refs.  3 and  4.  For  isolated  transmission  lines  the  direction  of  propagation  is  taken  to 
be  along  the  z direction.  Consistent  with  the  quasi-TEM  model  and  the  transmission-line 
wave-approach  described  in  Ref.  5,  the  z dependence  for  voltage  and  current  along  the 
transmission  line  is  accounted  for  by  a factor  e~iz,  where 


Manuscript  submitted  March  22,  1976. 


1 


B.E.  SPIELMAN 


7 " « + J'l 3,  (1) 

in  which  a is  the  attenuation  constant  due  to  conductor  and  dielectric  losses  and  0 is  the 
phase  constant.  Following  the  development  set  forth  in  Ref.  5,  the  attenuation  constant 
a is  given  by 


where  Pd  is  the  time-averaged  power  dissipated  per  unit  length  and  Pf  is  the  time-averaged 
power  flow  along  the  line. 

In  the  following  subsection  explicit  expressions  are  developed  for  the  quantities  Pd 
and  Pf.  It  is  these  expressions  which  were  used  to  provide  the  results  presented  in  the 
third  section  after  appropriate  incorporation  within  the  computer  programs,  described  in 
the  fifth  section  and  listed  in  Appendixes  A through  E. 


Isolated  Transmission  Lines 

To  obtain  a useful  expression  for  Pf  in  Eq.  (2),  the  following  treatment  is  employed. 
By  virtue  of  Eq.  (1)  a +2-traveling  wave  on  the  transmission  line  is  of  the  form 

V - V0e-(<*+jP)z,  !=■¥-.  (3) 


The  complex  power  flow  is  given  by 

IVo  I2 

Pf  = VI*  = e~2az  , (< 

Z0 

where  J*  and  Z*  are  the  conjugates  of  / and  the  characteristic  impedance  respectively. 
Then  Pf  is  given  by 


Pf  = (Re  Z0) 


e-2az 


As  is  shown  in  Appendix  F,  for  the  type  of  transmission  lines  considered  here,  Re  ZQ 
and  \Z0 1 are  nearly  equal  to _(Zq)ll,  the  characteristic  impedance  of  the  lossless  line.  By 
virtue  of  this  consideration  Pf  can  be  expressed  as 

Pf  = \V0\2  vCe~2az  , (6) 

\ here  v is  the  phase  velocity  and  C is  the  electrostatic  capacitance  per  unit  length. 
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In  the  next  two  subsections  of  this  report  explicit  expressions  are  developed  for  eval- 
uating isolated-transmission-line  conductor-loss  and  dielectric-loss  coefficients  ac  and  o :d, 
respectively.  The  total-loss  coefficient  a is  obtained  from  these  by  summing  ac  and  ctd. 


Conductor  Losses 


To  obtain  an  expression  which  is  useful  for  evaluating  Pd  in  Eq.  (2)  for_losses  due  to 
imperfect  conductors,  the  approximation  described  in  Ref.  6 is  employed.  Pd  can  be  ex- 
pressed approximately  by 


conductor 

surfaces 


(7) 


where  the  additional  subscript  c on  Pd  c denotes  power  losses  due  to  imperfect  conductors, 
ltf0l2  is  the  amplitude  squared  of  the’magnetic  field  at  conducting  surfaces  for  the  loss- 
less case,  and  R is  the  surface  resistance  of  the  metals  in  the  system.  For  good  conductors 
R can  be  written  as 


R = 


(8) 


where  f is  the  frequency  of  the  propagating  wave,  p0  is  the  free-space  permeability,  and 
o is  the  conductor’s  dc  conductivity. 

Consistent  with  the  quasi-TEM-propagation  assumption  for  transmission  media  in- 
homogeneously  filled  with  dielectric, 


H0 


X E0  , 


(9) 


where  rjeff  is  the  effective  intrinsic  impedance,  uz  is  the  unit  vector  in  the  z direction,  and 
Eq  is  the  electric  field  for  the  lossless  case  at  points  just  outside  the  conductor  surfaces. 
By  virtue  of  Eq.  (9) 


IEqI2 

^eff 


(10) 


where  coff  is  the  effective  permittivity,  determined  as  described  in  Ref.  3.  Using  the  results 
of  Eqs.  (10)  and  (8)  in  Eq.  (7),  Pd  c can  be  written  ae 


— — eelT  -l 


J I E0?  dZ- 

conductor 

surfaces 


(11) 
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In  the  lossless  TEM  solutions  described  in  Refs.  3 and  4,  the  electric  field  at  a point 
along  the  surface  of  a perfect  conductor  is  given  by 


E0  = 2? :q,  (12) 

where  q is  the  equivalent  charge  density  residing  at  the  conductor  surface  (sum  of  free 
and  polarization  charge  densities).  In  Ref.  3 the  free-charge-density  distribution  along 
the  surface  of  the  Afyth  conductor,  in  a system  having  Nc  conductor  surfaces,  is  approxi- 
mated by  a pulse  expansion  as 


<NS>y 


N, 
Q 1 


= V NJ  N.=  1 
Zj  9,  p.or  ' ’ - * 

i»l 


IV, 


(13) 


This  representation  arises  by  subdividing  the  contour  defined  by  tire  surface  of  the  NjU\ 
conductor  into  (Ns  )j  segments.  Along  the  ith  segment  the  free-charge-density  distribution 
is  taken  to  be  Constant  at  the  value  qfi . P(i)  is  a pulse  function  defined  by 


PNJ(0 


|=  1 on  the  ith  section  of  Nj 
I = 0 on  all  other  sections  of  N; 


(14a) 

(14b) 


Then  Pd  c in  Eq.  (11)  can  be  rewritten,  using  Eqs.  (12)  and  (13), 


ld,c 


N (N  ) 

^ ^ 1 £ (•>”')'  »■' 


(15) 


7=1  <=1 


no- 


where 1 is  the  length  of  the  ith  segment  on  the  IVyth  conductor. 

Finally  the  loss  coefficient  due  to  conductor  losses  ac  can  be  written  using  Eqs.  (15), 
(6),  and  (2)  as 


47T2e, 


a. 


J*fPo 

F a N c Wi  / N \2  N 

2 V [oj 1 ) A 1 (neper/unit  length) 

2 IV0 1 v C / j / j 


Mo 


(16) 


It  is  the  expression  given  in  Eq.  (16)  which  is  embodied  in  the  computer  programs  de- 
scribed later  in  this  report.  This  expression  has  been  used  to  provide  the  design  informa- 
tion for  conductor  losses  presented  later  in  the  report.  To  obtain  the  coefficients  ac  (in 
dB  per  unit  length)  the  following  relationship  is  provided  for  completeness 
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uc  (dB / unit  length)  = — — ae  (neper/unit  length). 


Dielectric  Losses 


To  obtain  an  expression  which  is  useful  for  evaluating  P^  in  Eq.  (2)  for  losses  due  to 
imperfect  dielectrics,  this  quantity  is  initially  written  as 


D r'r 

^d,d=  ^ ff  OH?  IeP  dS, 


1=1  Al 


where  the  second  subscript  d on  P^d  denotes  that  the  dissipation  losses  are  due  to  im- 
perfect dielectrics,  ND  is  the  number  of  imperfect  dielectric  regions  in  which  the  ith  re- 
gion has  a complex  permittivity  given  by 

£ = */  -je",  (19 


Aj  represents  the  transmission-line  cross-sectional  area  spanned  by  the  ith  simply-connected 
homogeneous  lossy  dielectric  region.  Equation  (18)  can  be  rewritten  as 


•'O 

Pdtd  ~ ^ ' Jj  2o)(tan  8;)  WCI-, 

1 1 A, 

where 


tan  5;  = —~ 


and  the  time-averaged  energy  stored  in  the  electric  field  in  the  ith  dielectric  region  is 
given  by 


%i  = d/2) 


\E  1 dS. 


In  Appendix  G details  are  presented  of  a development  which  shows  that  t'ne_time  average  ofthe 
total  energy  stored  in  the  electric  field  per  unit  length  of  transmission  line  Wei  is  related  to  Wei  by 


w . = e.  — — 

ei  e‘  bei  ■ 


TO  «■ 
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The  development  in  Appendix  G provides  a general  result  for  transmission  structures  with 
many  conductors  and  dielectrics  over  the  cross  section,  but  for  isolated  lines  and  coupled 
lines  treated  by  an  even-  and  odd-mode  two-port  interpretation  1VC  can  be  expressed  as 

i?e  = (l/2)ClVl2,  (24) 

where  C is  the  electrostatic  capacitance  of  tire  transmission  line  in  tire  two-port  configura- 
tion (isolated  and  even-  or  odd-mode  line)  and  V is  tire  voltage  associated  with  a +z- 
traveling  wave  on  the  line. 

To  facilitate  the  evaluation  of  a dielectric-loss  coefficient  in  terms  of  readily  comput- 
able parameters,  use  is  made  of  tire  definition  of  effective  relative  permittivity;  hence 

ecrf  - Sr  , (25) 

c0 

where  C0  is  the  electrostatic  capacitance  of  tire  transmission  line  under  consideration  but 
with  all  dielectric  materials  fictitiously  removed.  Consequently  Pd  d in  Eq.  (20)  can  be 
written  as 


, i2  ”2<**  decff 

Pd,d  = <oC0  I Vo  r<  2^  e,.(tan  Si)  . (26) 

»'« 1 

By  use  of  Eqs.  (26),  (2),  and  (6)  tire  loss  coefficient  due  to  imperfect  dielectrics  can  be 
written  as 


j 


nf 

<Xd  = 


3e0ff 

e.-(tan  Sf)  — — (neper/unit  length) 

OS; 


(27) 


or  in  units  of  dB  per  unit  length, 


20irf 

ud  = 

c 10 


3eefr 

ei(tan  si) 


(dB/unit  length). 


(28) 


In  Eqs.  (27)  and  (28)  c is  the  speed  of  light  in  free  space. 


To  evaluate  the  partial  derivative  in  Eqs.  (27)  and  (28),  a “forward”  difference 
quotient  is  employed;  thus 


3ecff  _ coff  ~ ceff 
3 Gj  Gj  — 


(29) 


i 

t 

i 
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where  is  the  value  of  eoff  for  the  structure  under  consideration  when  the  value  of  the 
relative  permittivity  for  the  ith  homogeneous  dielectric  region  is  perturbed  to  a slightly 
different  (higher)  value  e/.  Incorporating  Eq.  (29)  into  Eq.  (28)  provides  the  computa- 
tionally useful  result 


ad  « 


20  nf 


c Sn  10n/^ 


C;(tan  5f) 


/ ^eff  ~ ecff\ 

U j 


(dB/unit  length). 


(30) 


It  is  this  expression  which  has  been  incorporated  into  the  computer  programs  which  are 
documented  in  Appendixes  A through  E.  For  these  programs  the  index  i takes  the  value 
1,  since  there  is  only  one  lossy  homogeneous  dielectric  region  in  the  overall  inhomogeneous 
structures  treated  by  these  computer  programs. 

In  the  next  subsection  a description  is  given  for  the  determination  of  conductor-  and 
dielectric-loss  effects  in  coupled  line  structures,  with  emphasis  on  directional  couplers  em- 
ploying such  configurations. 


Coupled -Transmission-Line  Structures 

For  the  purposes  of  this  subsection  the  coupled  transmission  lines  will  be  treated  by 
an  even-  and  odd-mode  interpretation  [7] . Although  two  transmission  lines  coupled  over 
a given  length  truly  represent  a four-port  structure,  the  even  and  odd  modes  associated 
with  this  structure  are  each  two-port  transmission  lines  which  can  be  treated  separately. 

To  properly  assess  the  effects  of  losses  in  such  structures,  the  problem  is  twofold.  One 
problem  is  to  determine  tire  loss  coefficients  otc  and  ad  for  each  of  the  even  and  odd 
modes.  The  second  problem  is  to  determine  the  effects  of  these  coefficients  on  the  four- 
port  terminal  electrical  performance  of  the  entire  structure.  In  the  material  to  follow  the 
problem  of  determining  the  loss  coefficients  for  each  mode  will  be  discussed  first  and  will 
be  followed  by  a treatment  of  the  problem  of  determining  their  effects  explicitly  for  four- 
port  directional  coupler  losses. 


Even-  and  Odd-Mode  Loss  Coefficients 

Since  the  even  and  odd  modes  for  a coupled  line  structure  can  each  be  depicted  as 
two-port  transmission  lines  whose  length  is  that  of  the  coupled-line  region,  tire  loss  co- 
efficients for  conductor  and  dielectric  losses  can  be  determined  using  Eqs.  (16)  (or  (17)) 
and  (30)  respectively.  The  even-  and  odd-mode  loss  coefficients  for  losses  due  to  im- 
perfect conductors  can  be  written  as 


a. 


ck 


4tt2(ee  ff)/n  A/fro 

Mp  r ° 

2 !V0 12  vkCk 


Ailfy  (neper/unit  length), 


(31) 
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N 

where  k = e for  the  odd  mode  and  k = 0 for  the  odd  mode.  The  quantities  Vf,,  C/n  qlkJ , 
(^sff)h  A2$  are  evaluated  as  described  in  Ref.  3.  The  expression  for  even-  and  odd- 
mode loss  coefficients  for  losses  due  to  imperfect  dielectrics  can  be  expressed  as 


adh 


ft# 


20nf 


Nr 


c fin  l(H/(ecff),. 


2 

i»i 


e,(tan  8,-) 


(eeff)fc  - (ecff)fe 


<7  - 


(dB/unit  length),  (32) 


where  k = e for  the  even  mode  and  h - 0 for  the  odd  mode.  In  Eq.  (32)  (ecff)/.  represents 
the  effective  relative  permittivity  for  the  even  or  odd  mode  in  the  coupled-line  configura- 
tion to  be  analyzed.  This  value  is  determined  using  the  method  described  explicitly  in 
Ref.  3.  The  value  (e^f);.  is  determined  the  same  way  for  the  coupled-line  structure  in 
which  the  relative  permittivity  of  the  ith  dielectric  region  e,-  is  perturbed  to  the  slightly 
higher  value  e/.  Equations  (31)  and  (32)  have  been  used  to  evaluate  the  even-  and  odd- 
mode conductor-  and  dielectric-loss  coefficients  for  the  coupled-line  structures  described 
later  in  tins  report.  The  total-loss  coefficients  ae  and  a0,  for  the  even  and  odd  modes, 
are  obtained  by  adding  the  respective  values  of  occJ.  and  otdlr 


Effects  of  Losses  on  Four-Port  Terminal  Performance 

The  effects  of  losses  on  four-port  terminal  performance  of  directional  couplers  can 
be  understood  by  considering  Fig.  1.  This  figure  portrays  the  problem  as  being  that  of 
determining  the  loss  effects  on  measurable  terminal  voltages  bt  (i  = 1,  2,  3,  4)  once  the 
loss  coefficients  (orce,  ade)  and  (aco,  (xdo)  have  been  determined  for  the  modal  transmission 
lines  in  Figs,  lb  and  lc  respectively.  These  loss  coefficients  are  assumed  to  be  known  for 
this  development,  having  been  computed  by  the  method  described  in  the  previous  section. 

The  approach  used  here  starts  with  a procedure  similar  to  that  described  in  Ref.  7 . 
The  voltages  (i  - 1.  2,  3,  4)  are  written  as 


°l  = 
b2  = 
b3  = 

fe4  = 


r + r 

1 oe  1 oo 


r - r 

1 oe  1 oo 


Toe  ~ T00 


Too  + T00 


— icflecled  signal  ~ br, 

(33) 

= coupled  signal  = bc, 

(34) 

= isolated  signal  = ty, 

(35) 

= transmitted  signal  = bt, 

(36) 
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where  (I^,  Toe)  and  (roo,  T00)  are  pairs  of  reflection  and  transmission  coefficients  for  the 
even-  and  odd-mode  two-port  transmission  lines  respectively.  The  reflection  coefficients  for 
the  even  or  odd  mode  can  be  expressed  in  terms  of  two-port-network  parameters  A,  B,  C, 
and  D as 


Bj 

A-,  + „ - Cj  Z0  - Dj 

r of  - f . (37) 

A;  + ■=*■  + C;Z0  + Di 
Aq 


where  i = e for  the  even  mode  and  i = 0 for  the  odd  mode.  Similarly,  the  transmission 
coefficients  can  be  expressed  as 


= % — . (38> 

Ai  + ~Z~  + CiZ  0 + Di 

where  again  i = e for  the  even  mode  and  i = 0 for  the  odd  mode.  In  Eqs.  (37)  and  (38) 
ZQ  represents  the  characteristic  impedance  of  the  input  and  output  ports  shown  in  Fig.  1 
( Z0  « V Z0CZ0O).  The  parameters  A,  B,  C,  and  D for  the  lossy  even-mode  two-port  are 
given  by 


Ad  = De  = cosh  yeB,  (39) 

Be  - Zoe  sinh  ye2,  (40) 

Ce  =^~-jsinh  7e2.  (41) 

Similarly,  for  the  odd  mode 

A0=  D0  = cosh  702,  (42) 

B0  = Z00  sinh  70e>  (43) 

C0  sinh  702.  (44) 

In  Eqs.  (39)  through  (44)  2 represents  die  physical  length  of  the  coupled-line  region,  Zoe 
and  Z00  are  the  characteristic  impedances  of  the  even-  and  odd-mode  transmission  lines 
respectively,  and  ye  and  y0  are  the  even-  and  odd-mode  propagation  constants,  taken  for 
this  development  to  be  given  by 


*> 

CO. 

+ 

e 

H 

(45) 

To  “ + $0  > 

(46) 
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where  cce  and  of0  are  the  total  (known)  even-  and  odd-mode  loss  coefficients  and  0C  and 
|S0  are  the  even-  and  odd-mode  phase  constants  (known)  respectively,  which  are  determined 
by  the  lossless  analysis  described  in  detail  in  Ref.  3.  The  generality  of  this  development 
permits  Pe  and  (30  to  be  determined  with  allowance  for  even  and  odd  modes  having  different 
phase  velocities. 

By  properly  combining  Eqs.  (33)  through  (44),  the  measurable  terminal  voltages  of 
the  four-port  coupler  under  analyses  can  be  expressed  as 


bt  = 


bi  ~ 


f tanh  7c2 

tanh  702 

l 2 + Fx  tanh  7C2 

2 + F-y  tanh  y0Z 

sech  7e2 

sech  702 

2 + F-y  tanh  ye% 

2 + Fj  tanh  7D2  ’ 

sech  7c2 

sech  702 

2 + Fl  tanh  7e2 

2 + Fj  tanh  7G2  ’ 

/ tanh  7e£ 

tanh  70C 

2 H-  Fx  tanh  je9.  2 + F1  tanh  7G2 
In  Eqs.  (47)  through  (50)  Fl  and  F2  are  quantities  given  by 


(47) 


(48) 


(49) 


(50) 


a/2 


, 1/2 


Fx  ~\  z, 


oe 

t 

'00  t 


JQO 


Joe 


(51) 


,1/2 


1/2] 


F2  = 1/2 


Joe 


voo  / 


uoo 


-'oe  t 


(52) 
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bi  = ~j 


2otc%  + F-y  2ot0Z  + F1 


(55) 


6r  = Fz  ( 2ac2  + F1 


2a02+F1 


(56) 


In  Eqs.  (53)  through  (56)  the  even-  and  odd-mode  loss  coefficients  appear  explicitly. 
Also,  these  relationships  reveal  that,  even  for  perfectly  matched  couplers,  for  otc  different 
from  a0  the  isolated  and  reflected  signals  are  not  zero.  For  couplers  with  other  sources 
of  impedance  match  and  isolation  degradation  (such  as  different  even-  and  odd-mode  phase 
velocities)  the  additional  degradation  of  isolation  and  reflected  signal  due  to  different  modal 
loss  characteristics  will  be  superimposed  on  the  other  effects.  The  relationships  developed 
in  this  section  have  been  tested  on  experimental  coupler  models,  and  the  results  are  de- 
scribed in  the  next  section. 


i 


EXAMPLES  OF  LOSS  EVALUATIONS  AND 
DESIGN  INFORMATION 

In  this  section  results  are  presented  for  specific  structures  of  isolated  and  coupled 
transmission  lines.  The  examples  of  microstrip  and  coplanar  waveguide  are  presented  to 
demonstrate  the  accuracy  of  the  computer-aided  analysis  from  the  preceding  section  and 
to  provide  useful  design  information  for  engineering  applications.  The  results  for  inverted 
microstrip  and  trapped  inverted  microstrip  serve  to  provide  useful,  currently  unavailable 
design  information  for  more  complicated  structures  which  offer  potential  for  circuit  ap- 
plications. Following  these  results  is  a comparison  of  the  loss  characteristics  for  the  four 
preceding  transmission  lines.  The  results  for  the  edge-coupled  microstrip  structure  with  a 
dielectric  overlay  confirm  the  utility  of  the  coupled-line  loss  analysis  described  in  the  pre- 
ceding section  and  provide  currently  unavailable  design  information  for  this  structure. 

This  structure  is  presently  the  most  viable  approach  for  providing  broadband  couplers,  fil- 
ters, and  Schiffman-ph^e-shift  sections  in  a microstrip-compatible  format. 


Isolated  Transmission  Lines 
Microstrip 

To  illustrate  the  approach  that  resulted  in  Eq.  (16),  the  conductor-loss  coefficient 
was  computed  for  a microstrip  line  with  an  impedance  of  50  ohms  (nominally).  A gener- 
ic cross  section  of  such  a line  is  portrayed  in  Fig.  2.  The  dc  conductivity,  for  Eq.  (16), 
was  selected  to  be  4.10  X 107  mho/m  based  on  data  provided  in  Ref.  8,  corresponding  to 
the  value  designated  for  gold.  The  values  of  ac  comp'  J for  this  configuration  are  repre- 
sented by  the  solid  curve  in  Fig.  3.  The  reference  values  shown,  represented  by  a broken 
curve,  are  due  to  a computer-program  version  of  Schneider’s  results  [2] . Schneider’s  re- 
sults were  selected  as  a reference  because  they  have  compared  favorably  with  the  experi- 
mental data  provided  in  Refs.  9 and  10.  The  agreement  here  is  reasonable,  with  the  values 
computed  by  the  method  proposed  in  this  report  being  only  a few  hundredths  dB  per 
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[wavelength  higher  than  the  reference.  This  discrepancy  can  be  attributed  to  the  discrete 

nature  of  the  equivalent  charge  densities  used  in  the  formulation.  This  effect  is  somewhat 
t equivalent  to  a surface-roughness  effect. 


j 

i 


METALLIZATION  THICKNESS 
Tm  = 6.35/i.m 


LINE  WIDTH 
J*W=0.6lmm— | 


SUBSTRATE  WITH  DIELECTRIC  CONSTANT  K=«,  =10.0 


T 


H=0.635mm 


i 


Fig.  2 — Generic  cross  section  of  nominal-50-ohm  microstrip  line 


Fig.  3 — Conductor  attenuation  cotfiicient  for  microstrip  as  configured  in  Fig.  2 
(a  configuration  with  Zq  = 50  ohms). 


! 

I 

t 


The  method  for  evaluating  dielectric-loss  effects  was  tested  for  the  same  microstrip 
configuration.  Equation  (30)  was  employed  with  ND  equal  to  1 and  the  loss  tangent  value, 
tan  Sj,  selected  to  be  6 X 10~4  (alumina).  The  results  computed  using  Eq.  (30)  are  repre- 
sented in  Fig.  4 by  the  solid  curve.  Shown  for  comparison  are  reference  values  due  to 
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Pucel  et  al.  [9] , which  are  in  excellent  agreement.  The  method  described  for  computing 
dielectric  losses  that  resulted  in  Eq.  (30)  is  similar  to  that  of  Schneider  [11] . In  Ref.  12 
Simpson  and  Tseng  show  additional  results  which  agree  within  1%  (typically)  compared  to 
results  computed  by  this  method. 

The  total-loss  coefficient  (sum  of  ac  and  ad)  for  the  microstrip  line  under  consider- 
ation will  later  be  compared  (in  Fig.  15)  with  similar  characteristics  for  nominally-50-ohm 
configurations  of  coplanar  waveguide,  inverted  microstrip,  and  trapped  inverted  microstrip. 


REFERENCE v ^ 


8 10  12 
FREQUENCY  (GHz) 


16  18 


Fig.  4 — 


Dielectric-loss  coefficient  for  microstrip  as  configured  in  Fig.  2 
(a  configuration  with  ZQ  = 50  ohms). 


Coplanar  Waveguide 

The  second  illustrative  example  chosen  to  demonstrate  the  effectiveness  of  the  meth- 
ods that  were  described  in  the  analysis  section  is  coplanar  waveguide  [1]  with  the  cross 
section  depicted  in  Fig.  5.  The  computations  here  were  for  an  approximately-5C-ohm  con- 
figuration. The  metallization  thickness  was  again  selected  to  be  6.35  pm  with  a dc  con- 
ductivity of  4.10  X 107  mho/m.  Figure  6 shows  the  computed  cc  for  this  case  plotted 
versus  frequency  as  a solid  curve.  The  computed  ctd  for  this  configuration  is  also  shown 
in  Fig.  6,  plotted  as  a broken  curve.  The  total-loss  coefficient,  obtained  by  summing  o:c 
and  ad,  is  shown  plotted  in  Fig.  7 (and  Fig.  15).  Also  shown  in  Fig.  7 are  measured 
values  [13]  for  comparison.  These  values  were  obtained  from  end-coupled  resonant-section 
transmission  measurements.  The  agreement  in  this  frequency  range  was  reasonable,  being 
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within  0.1  dB  per  wavelength  over  most  of  the  range  from  4 to  11  GHz.  Equally  good 
agreement  (about  0.1  dB/wavelength  or  better)  was  found  between  the  computed  values 
of  total  loss  and  values  obtained  by  a straight-through  transmission-line  measurement  [13] . 
This  correlation  in  values  was  obtained  from  4 GHz  through  about  10  GHz. 
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Fig-  6 — Conductor-  and  dielectric-loss  coefficients  for  coplanar  waveguide 
configured  (as  in  Fig.  5)  with  ZQ  = 50  ohms 
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Fig.  7 — Total  loss  coofficient  versus  frequency 
for  coplanar  waveguide 

Additional  design  information  developed  for  coplanar  waveguide,  using  the  methods 
described  in  the  analysis  section,  are  shown  in  Table  1.  The  information  in  this  table  re- 
lating to  configuration  geometry  and  material  characteristics  is  consistent  with  the  nota- 
tion defined  by  Fig.  5.  The  ratios  of  a1{b1  shown  were  selected  to  span  a range  of  values 
that  should  impact  on  a broad  range  of  applications  for  this  transmission  line. 


Inverted  Microstrip 

Inverted  microstrip  is  the  third  transmission  line  which  has  been  analyzed  using  tire 
method  described  in  the  analysis  section.  The  generic  cross  section  of  inverted  miciostrip 
is  depicted  in  Fig.  8.  Table  2 shows  the  analysis  results  for  aspect  ratios  (W/H)  spanning 
a range  of  values  from  1 to  8.  This  range  was  chosen  to  provide  design  information  for 
a broad  range  of  circuit  applications.  The  substrate  dielectric  constant  chosen  for  the  set 
of  data  corresponds  to  that  of  fused  silica,  which  is  deemed  to  be  a suitable  material  for 
this  transmission  line  for  applications  at  higher  microwave  and  millimeter-wave  frequencies. 
The  conductor-  and  dielectric-loss  coefficients  evaluated  for  this  sequence  of  transmission- 
line configurations  are  plotted  as  a function  of  aspect  ratio  W/H  in  Fig.  9 for  design  usage. 
The  characteristic  impedances  and  phase  velocities  for  configurations  spanning  the  range 
of  aspect  ratios  W/H  from  1 to  8 are  plotted  in  Fig.  10,  readily  usable  for  circuit  design 
applications. 

To  demonstrate  the  correlation  between  computed  and  measured  characteristics,  the 
computed  phase  velocity  is  compared  to  experimental  data  for  inverted  microstrip.  Figure 
11  shows  computed  values  of  Veeff  versus  aspect  ratio  W/H  as  a solid  curve.  Measured 
points  shown  in  this  figure  were  obtained  via  time-domain  reflectometer  measurements  [2] 
The  errors  between  measured  and  computed  values  are  within  3%  over  this  range  of  aspect 
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ratios.  (The  total-dissipation-loss  coefficient  for  a nominally-50-ohm  inverted-microstrip 
configuration  is  plotted  in  Fig.  15  for  comparison  with  the  other  transmission  lines.) 


Table  1 — Design  Information  for  Coplanar  Waveguide  as  Shown  in  Fig.  5 But  for  a 
Range  of  W and  S Values  Rather  Than  the  Nominal-50-Ohm  Values  (and  With  H 
= 0.635  mm,  = 10.0,  e'  = 10.1,  and  tan  Sj  = 0.0006) 


w 

(mm) 

S 

(mm) 

Zo 

(ohms)* 

V 

(108  m/s)* 

«C/VT 

(10-4  dB/mVHz  ) 

<*d/f 

(10-11  dB/m  Hz) 

0.118 

0.529 

102 

1.36 

1.32 

10.0 

0.353 

0.412 

72.2 

1.37 

1.22 

11.82 

0.588 

0.294 

57.5 

1.38 

1.51 

9.41 

0.823 

0.176 

46.2 

1.38 

2.29 

7.56 

Tabulated  electrical  characteristic  determined  by  the  method  described  in  Ref.  3.  This  characteristic  is 
furnished  in  the  output  of  computer  program  furnished  in  Appendix  B. 


Tm  — 6.35 


z — GROUND  PLANE 

Fig.  8 — Generic  cross  section  of  nominal-50-ohm  inverted  microstrip 
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1 Table  2 — Design  Information  for  Inverted  Microstrip  as  Shown  in  Fig.  8 But  for  a 

j Range  of  W/H  Values  (with  H,  T,  Tm,  and  ej  as  Shown  and  = 3.88,  tan  51  = 

0.0002,  and  a Metal  dc  Conductivity  of  4.10  X 107  mho/m.) 


W/H 

w 

(mm) 

*0 

(ohms)* 

V 

(108  m/s)* 

a c/vT 

(10-5  dB/m-s/Hz  ) 

<*dlf 

(10-12  dB/m  Hz) 

1 

0.381 

101.8 

2.41 

2.74 

8.40 

2 

0.762 

75.1 

2.52 

2.37 

6.03 

3 

1.143 

60.8 

2.60 

2.26 

5.02 

4 

1.524 

51.4 

2.66 

2.19 

4.24 

5 

1.905 

44.7 

2.70 

2.16 

3.69 

n 

2.286 

39.6 

2.73 

2.14 

3.27 

H 

2.667 

35.6 

2.76 

2.13 

2.94 

H 

3.048 

32.4 

2.78 

2.11 

2.67 

♦Tabulated  electrical  characteristic  determined  by  the  method  described  in  Ref.  3.  This  characteristic  is 
furnished  in  the  output  of  the  computer  program  furnished  in  Appendix  C. 


12  3 4 5 6 7 8 

W/H 

Fig.  9 — Loss  constants  vs  aspect  ratio  for  trapped 
inverted  microstrip  (Fig.  8) 
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Trapped  Inverte  d Microstrip 

This  section  provides  design  information  for  trapped  inverted  microstrip,  characterized 
by  tiie  generic  cross-section  configuration  shown  in  Fig.  12.  Table  3 provides  details  of 
configuration  geometry  and  material  parameters  for  the  cross  sections  analyzed  in  this  work. 
So  that  this  information  can  be  used  for  design  purposes,  the  conductor-  and  dielectric-loss 
coefficients  are  plotted  in  Fig.  13  for  aspect  ratios  W/H  ranging  from  0.25  to  6.0.  The 
somewhat  irregular  appearance  of  these  characteristics  can  be  attributed  to  the  effects  of 
channel-sidewall  interaction  with  the  conducting  strip,  which  is  not  explicitly  accounted 
for  by  the  “aspect  ratio”  definition  as  W/H.  A similar  situation  is  apparent  in  Fig.  14, 
where  the  characteristic  impedance  Z0  and  phase  velocity  t>_  are  shown  for  the  same  range 
of  W/H  values. 

The  total  dissipation-loss  coefficient  for  a nominally-50-ohm  configuration  of  trapped 
inverted  microstrip  is  plotted  in  Fig.  15  for  comparison  with  those  of  microstrip,  coplanar 
waveguide,  and  inverted  microstrip. 


Comparison  of  Loss  Characteristics 

This  section  presents  u comparison  of  the  loss  characteristics  for  specific  nominally- 
50-ohm  configurations  of  microstrip,  coplanar  waveguide,  inverted  microstrip,  and  trapped 
inverted  microstrip  transmission  lines.  The  total-loss  coefficients  (sum  of  the  conductor- 
and  dielectric-loss  coefficients)  for  these  lines  are  shown  in  Fig.  15  for  the  frequency  range 
from  0 to  60  GHz.  Effects  and  consideration  of  higher  order  moding  have  been  neglected, 
as  well  as  effects  due  to  radiation  losses.  Details  defining  the  configurations  of  microstrip 
and  coplanar  waveguide  analyzed  to  provide  the  corresponding  curves  in  this  figure  are 
shown  in  Figs.  2 and  5 respectively,  with  the  dc  conductivity  o selected  to  be  4.10  X 107 
mho/m  for  both  configurations.  The  configuration  of  the  inverted  microstrip  line  analyzed 
is  defined  by  Fig.  8,  with  W ~ 1,676  mm,  = 3.88,  a - 4.10  X 107  mho/m,  and  tan  5^ 

= 2 X 10-4 . The  configuration  for  the  trapped  inverted  microstrip  line  analyzed  is  specified 
by  the  following  values  in  conjunction  with  the  nomenclature  shown  in  Fig.  12:  W = 1.52 
mm,  S = 1.27  mm,  H - 0.508  mm,  Tm  = 6.35  pm,  ex  = 10.0,  ei  = 10.1,  o = 4.10  X 107  mho/m, 
and  tan  = 6 X 10”4. 
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Fig.  12  — Generic  cross  section  for  nominal-50-ohm  trapped  inverted  microstrip 


20 


NRL  REPORT  8009 


Table  3 — Design  Information  for  Trapped  Inverted  Microstrip  as  Shown  in  Fig.  12 
But  for  a Range  of  IF///  Values  (With  H,  T,  Tm , and  ej  as  Shown  and  With 
= 10.1,  tan  5j^  4 0.0006,  and  a Metal  dc  Conductivity  of  4.10  X 107  mho/m) 


4 2.03 


Zo 

(ohms)* 

(10"^  m/s)* 

«c/VT 

(10-5  dB/mv/Hz  ) 

(10“n  dB/m  Hz) 

122 

1.77 

3.61 

5.81 

102 

1.83 

2.72 

5.32 

80.7 

1.92 

2.25 

4.76 

61.6 

1.97 

2.22 

4.44 

49.0 

2.12 

1.97 

3.85 

41.6 

2.18 

2.00 

4.08 

35.0 

2.18 

2.19 

4.02 

29.5 

2.55 

4.35 

in  Appendix  D. 


Fig.  13  — Loss  constants  vs  aspect  ratio  for 
trapped  inverted  microstrip  (Fig.  12) 
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Fig.  14  — Characteristic  impedance  <£q  and  phase  ve- 
locity vp  vs  W/H  for  trapped  inverted  microstrip 
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- Comparison  of  total  dissipative  losses  for 
four  50-ohm  transmission  lines 


The  microstrip  configuration  chosen  here  represents  a standard  configuration  of  this 
transmission  line  as  it  has  been  used  in  many  applications  at  frequencies  up  through  about 
12  GHz.  This  line  is  representative  of  the  transmission  line  fabricated  using  an  alumina 
substrate  with  a predominantly  gold  metallization  system.  The  same  is  true  for  the  selected 
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coplanar  waveguide  configuration.  The  configuration  for  inverted  microstrip  was  chosen 
to  represent  a version  fabricated  on  a fused-silica  substrate,  again  with  gold  being  the  pre- 
dominant carrier  of  RF  current  in  the  metallization  system.  The  configuration  for  trapped 
inverted  microstrip  represents  a model  fabricated  on  a standard  alumina  substrate,  metal- 
lized with  a predominantly  gold  system. 


Coupled  Transmission  Lines 

Edge-Coupled  Microstrip  with  a Dielectric  Overlay 

This  section  serves  a twofold  purpose.  The  first  purpose  is  to  provide  an  illustrative 
example  demonstrating  the  applicability  of  the  coupled-line  loss  analysis  developed  in  the 
analysis  sections.  The  second  purpose  is  to  provide  a quantitative  assessment  of  the  losses 
encountered  in  quarter-wavelength-long  (at  midband)  sections  of  edge-coupled  microstrip 
line  with  a dielectric  overlay  for  various  coupling  levels.  This  coupled-line  structure  is  cur- 
rently the  most  viable  approach  for  providing  high-performance  broadband  couplers,  fil- 
ters, and  Schiffman-phase-shift  sections  in  a rnicrostrip-compatible  fbrmat  [14] . 

The  generic  cross  section  of  the  edge-coupled  microstrip  with  dielectric  overlay  is 
portrayed  in  Fig.  16.  Four  configurations  differing  in  geometric  parameters  were  analyzed 
(Table  4).  Table  5 lists  the  computed  electrical  parameters  for  these  sections. 

Sections  1 and  2 in  these  tables  were  fabricated  into  a two-section,  asymmetric  cou- 
pler similar  to  that  shown  in  Fig.  17.  The  circuit  pattern  was  etched  on  a portion  of  alu- 
mina substrate  with  dimensions  23  mm  by  10  mm  by  0.635  mm.  The  metallization  was 
chromium-gold  with  a thickness  of  6.35  pm.  The  overlays  were  made  of  alumina  pieces 
on  each  section  and  were  attached  to  the  substrate  using  Stycast  Hi  K epoxy  (K  = 10). 

The  two-section  coupler,  over  the  frequency  band  from  2 to  8.5  GHz,  has  a nominal  cou- 
pling value  of  6.7  dB.  Using  Eqs.  (53)  through  (56)  and  the  computed  values  of  loss  co- 
efficient shown  in  Table  5 the  calculated  dissipation  loss  for  the  two  coupled-line  sections 
was  found  to  be  0.2  dB,  obtained  for  the  coupler  midband  frequency  of  5.4  GHz.  When 
measured  losses  due  to  connectors  and  theoretical  lead-in-line  losses  are  added,  the  total 
dissipation  loss  determined  is  0.42  dB.  This  agrees  well  with  the  total  measured  coupler 
dissipation  of  0.4  dB. 

Sections  3 and  4 in  Tables  4 and  o were  combined  in  the  fabrication  of  a two-section 
asymmetric  coupler  similar  to  the  coupler  shown  in  Fig.  17.  The  experimental  model  had 
a nominal  coupling  characteristic  of  20  dB  over  the  frequency  band  from  2 to  8.5  GHz. 
Using  the  evaluation  scheme  described  for  the  coupler  in  Fig.  17,  the  dissipation  loss  for 
the  coupled  line  sections  was  determined  to  be  0.08  dB.  After  adding  contributions  for 
connector  and  lead-in  line  lengths,  the  total  dissipation  loss  was  computed  to  be  0.38  dB. 
This  agrees  well  with  the  measured  midband  (5.4  GHz)  diss>  >ation  loss  of  0.4  dB. 

Another  portrayal  of  the  loss  characteristics  computed  for  the  four  coupled-line  sec- 
tions of  Tables  4 and  5 is  shown  in  Figs.  18  and  19.  Figure  18  shows  plots  of  otco  and 
ado  versus  frequency  for  the  coupling  sections  used  in  the  6.7-dB  coupler.  Also  shown 
are  the  ratios  <xcol<xce  and  otdo/ade.  A large  disparity  exists  between  the  even-  and  odd-mode 
conductor  losses  for  the  tight  section.  Figure  19  shows  similar  data  for  the  coupling  sec- 
tions used  to  make  up  the  20-dB  coupler. 
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Fig.  16  — Generic  cross  section  of  nominal-50-ohm  edge-coupled  microstrip 
with  a dielectric  overlay 


Table  4 — Parameters  Which  Along  With  H>  L> 
Tm>  and  ej  as  Shown  in  Fig.  16  and  With  = 
10.1,  tan  Sjl  = 0.0006,  and  a = 4.10  X 107 
mho/m  Define  Four  Configurations  of  Edge- 
Coupled  Microstrip  With  Dielectric  Overlay 


Section 

W 

(mm) 

S 

(mm) 

1 

0.223 

0.0305 

2 

0.483 

0.541 

3 

0.483 

0.927 

4 

0.467 

2.31 

Table  5 — Computed  Electrical  Parameters  for  the  Coupled-Line 
Sections  Described  by  Table  4 


Section 

Lossless 

Miliband 

Coupling 

(<1B) 

Zo 

(ohms) 

^ oc 

(ohms) 

z 

*'oo 

(ohms) 

<VvT 

(10-5  OB/myTT) 

“co/vT" 

(10->  dB/m\/Hz  ) 

«dc'( 

U0->°  (lB/m  Hz) 

(1<H°  dB/m  Hz) 

1 

3.3 

50.0 

115 

21.7 

5.60 

76.9 

1.20 

i 59 

2 

12.8 

50.7 

64.1 

40.2 

6.69 

10.3 

1.36 

1 38 

3 

16.9 

50.7 

58.6 

43.9 

7.59 

9.31 

1.60 

1.38 

4 

26.8 

49.4 

51.7 

47.2 

9.77 

10.2 

1.90 

1.13 
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(b)  Loose  section  (section  4) 


Fig.  19  — Loss  coefficients  for  the  sections  of  the  20-dB 
two-section  coupler  that  combines  sections  3 and  4 of 
Tables  4 and  5 
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DISCUSSION 

This  report  details  and  illustrates  analyses,  amenable  to  computer  programming,  which 
can  evaluate  dissipation  losses  in  isolated  or  coupled  microwave  and  millimeter-wave  trans- 
mission lines.  These  lines  can  be  composed  of  both  imperfect  conductors  and  piecewise 
homogeneous  dielectric  materials.  The  analyses  described  here  have  been  implemented  in 
the  form  of  digital  computer  programs  for  treating  microstrip,  coplanar  waveguide,  inverted 
microstrip,  trapped  inverted  microstrip,  and  edge-coupled  microstrip  with  a dielectric  overlay, 

The  programs  for  microstrip  and  coplanar  waveguide  were  used  to  analyze  various  con- 
figurations for  these  transmission  lines.  The  computed  values  of  conductor  loss  for  micro- 
strip on  alumina  were  compared  with  well-accepted  values  due  to  Schneider  [2]  and  were 
found  to  agree  within  better  than  1.6  dB/m  for  frequencies  up  through  18  GHz.  The  com- 
puted values  for  dielectric  losses  in  microstrip  agree  with  Simpson  and  Tseng  [12]  and 
Schneider  [11]  and  agreed  to  within  3%  or  better  at  frequencies  up  to  18  GHz.  The  sum 
of  conductor  and  dielectric  losses  computed  for  coplanar  waveguide  agreed  with  the  ex- 
perimental values  due  to  McDadeand  Stockman  [13]  to  within  0.1  dB/wavelength.  The 
experimental  determinations  were  made  by  two  independent  techniques. 

Conductor-  and  dielectric-loss  coefficients  were  also  evaluated  for  inverted  microstrip 
and  trapped  inverted  microstrip.  This  information  was  presented  in  the  form  of  useful 
design  curves,  providing  loss  values  as  functions  of  the  aspect  ratios  for  these  lines.  There 
is  encouraging  agreement  (better  than  2%)  between  computed  phase  velocities  (v^ff ) for 
inverted  microstrip  and  measured  data  [2] . These  calculations  are  made  using  the  same 
equivalent  charge  densities  that  are  used  in  the  loss  calculations. 

The  computed  total  loss  coefficients  for  the  four  types  of  lines  are  compared  in  Fig. 

15.  The  losses  incurred  in  inverted  and  trapped  inverted  microstrip  at  frequencies  as  high 
as  60  GHz  are  comparable  to  those  incurred  in  microstrip  at  frequencies  in  the  range  5 
to  10  GHz.  Also,  coplanar  waveguide  appears  to  be  considerably  more  lossy  than  micro- 
strip. The  structures  compared  in  this  figure  were  selected  to  have  characteristic  impedance 
levels  of  nominally  50  ohms.  Consistent  with  this  characteristic  the  conducting  strip  widths 
required  for  microstrip  and  coplanar  waveguide  were  approximately  2-1/2  to  3 times  nar- 
rower than  those  required  for  inverted  and  trapped  inverted  microstrip.  This  feature  en- 
hances the  attractiveness  of  the  inverted  and  trapped  inverted  microstrip  lines  by  virtue  of 
the  mitigation  of  fabrication  difficulties.  By  concentrating  more  of  the  field  energies  in 
air  (with  correspondingly  lower  eeff),  wider  strips  are  possible  for  a prescribed  impedance 
level  (compared  to  the  other  two  lines).  These  advantages  should  be  realized  even  if  the 
dimensions  must  be  contracted  to  suppress  higher  order  moding  at  the  higher  frequencies. 
This  should  be  investigated  more  extensively. 

The  analyses  for  computing  conductor-  and  dielectric  loss  coefficients  and  for  relating 
these  parameters  to  terminal  electrical  characteristics  were  applied  to  edge-coupled  micro- 
strip with  a dielectric  overlay.  Computed  results  were  compared  to  measured  characteristics 
for  two  experimental  models  of  two-section  couplers,  having  nominal  coupling  values  of 
6.7  dB  and  20  dB  respectively.  Total  dissipation  values  computed  for  *hese  couplers  agreed 
with  experimental  evaluations  within  0.07  dB. 
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The  difference  between  the  even-  and  odd-mode  loss  coefficients  for  this  structure 
gives  rise  to  a finite  isolation  for  an  otherwise  perfect  coupled  line  section.  For  the  6.7- 
dB  coupler  considered  here  the  tight  and  loose  sections  are  limited  to  maximum  isolation 
levels  of  42  dB  and  59  dB  respectively  due  to  this  phenomenon.  The  tight  and  loose  sec- 
tions of  the  20-dB  coupler  are  similarly  limited  to  70-dB  and  92-dB  isolation  levels 
respectively. 

Several  factors  can  contribute  to  error  in  the  analyses  described  here.  One  is  the 
discrete  representation  of  charge-density  distributions  employed  in  the  quasi-TEM  model 
employed.  In  one  sense  this  discretization  can  be  viewed  as  a “surface  roughness”  which 
could  lead  to  high  estimates  of  losses.  Mathematically  smoothing  the  charge-density  dis- 
tributions might  lead  to  improvement.  Another  source  of  error  is  the  approximation  shown 
in  Eq.  (9).  A method  which  would  reduce  this  error  is  to  solve  the  magnetostatic  problem 
of  evaluating  the  magnetic  field  at  conductor  surfaces  for  these  transmission  media.  Cer- 
tainly the  difference-quotient  approximation,  shown  in  Eq.  (29),  introduces  error  in  the 
dielectric-loss  evaluations.  Also,  judicious  choices  must  be  made  for  the  dc  conductivity 
of  the  metal  system  and  the  loss  tangents  for  dielectrics. 

In  the  next  section  documentation  is  given  for  the  user-oriented  computer  programs 
which  can  provide  design  information  for  the  structures  treated  in  earlier  sections. 


COMPUTER  PROGRAMS 

This  section  de  bribes  the  computer  programs  MS,  CPW,  IM,  TIM,  and  LEMDOC, 
listed  in  Appendixes  A,  B,  C,  D,  and  E.  These  programs  compute  the  electrical  character- 
istics, including  loss  characteristics,  for  microstrip,  coplanar  waveguide,  inverted  microstrip, 
trapped  inverted  microstrip,  and  edge-coupled  microstrip  with  a dielectric  overlay 
respectively. 


Program  MS  (Microstrip) 

Program  MS,  which  computes  the  characteristics  for  micrcstrip  (Fig.  2),  is  written  in 
Fortran  IV  and  is  compatible  With  the  CDC  3800  digital  computer  system.  With  minor 
modifications  this  program  should  also  be  suitable  for  use  with  the  CDC  6000-series  com- 
puter systems,  with  the  IBM  360  or  370  systems,  and  with  any  other  computer  having 
comparable  storage  capability  and  compatibility  with  the  Fortran  IV  language.  The  core 
storage  necessary  for  executing  this  program  on  the  CDC  3800  system  at  the  NRL  Re- 
search Computation  Center  is  42,000  words.  Although  the  total  core  storage  of  this 
computer  is  now  75,536  words,  the  maximum  storage  available  for  a single  array  (without 
resorting  to  special  array -handling  techniques)  is  32,768  words.  It  is  this  system  require- 
ment that  constrains  the  largest  array  in  MS  to  be  181  X 181  or  32,761  words. 

Furthermore,  because  of  the  standard  loading  procedures  employed  by  the  CDC  3800 
system  at  NRL,  the  following  scheme  for  handling  large  arrays  was  incorporated  into  MS. 
The  two  relatively  large  arrays  in  MS  are  array  A,  181  X 181,  and  array  Al,  132  X 132. 
These  arrays  appear  in  different  subroutines  and  are  used  consecutively.  For  these  reasons 
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the  arrays  A and  A1  are  placed  in  a block  of  COMMON  storage  labeled  HELP.  During 
the  loading  of  the  program,  HELP  is  loaded  into  a storage  bank  with  32,768  storage  ad- 
dresses available  to  accommodate  array  A.  This  is  done  using  a BANK  card.  In  addition 
another  BANK  card  forces  main  program  MS  and  subroutine  MSCUPF  into  a different 
storage  bank.  The  storage  required  to  operate  using  A and  A1  is  minimized  by  forcing 
them  to  share  the  same  storage  locations  by  means  of  the  statement  EQUIVALENCE 
(A,  Al).  These  techniques  are  necessary  to  properly  load  MS  into  the  CDC  3800  system 
at  NRL. 

The  length  of  time  necessary  for  program  execution  is  approximately  9 minutes  for 
each  configuration  to  be  analyzed. 


Input  Information 


Each  configuration  of  microstrip  to  be  analyzed  by  program  MS  is  characterized  by 
specifying  the  cross  section  in  terms  of  the  quantities  W,  H,  L,  , tan  , and  o (Fig.  2). 
Accordingly,  the  first  data  card  for  any  execution  of  MS  contains  the  quantity  NSETS. 
This  is  an  integer  specifying  the  number  of  configurations  to  be  analyzed.  NSETS  is 
punched  on  the  first  data  card  according  to  the  formal;  110. 

Each  configuration  included  in  NSETS  is  specified  on  a separate  data  card  located 
consecutively  behind  the  first  data  card.  Each  card  lists  the  following  six  quantities 
punched  according  to  the  format  F10.6: 


W = width  of  the  conducting  strip  in  mils  (1  mil  = 25.4  pm)', 

H = height  of  the  conducting  strip  above  the  ground  plane  (the  substrate 
thickness)  in  mils; 

L = substrate  width  in  inches,  which  is  a variable  quantity  to  allow  for  nar- 
row substrate  widths,  with  a value  L = 1.0  inch  (2.54  cm)  having  been 
adequate  to  eliminate  effects  due  to  finite  substrate  width; 


ER  ~ relative  permittivity  of  the  substrate  material: 

TAND1  = loss  tangent  of  the  substrate  dielectric  (X  104); 

SIGMA  = dc  Conductivity  of  the  metallization  in  mhos  per  meter  (X  10“7). 


The  number  of  cards  behind  the  first  data  card,  each  containing  sets  of  these  six  parame- 
ters, should  equal  the  value  read  in  for  NSETS. 


Output  Information 

The  first  line  of  output  is  invariant  for  each  configuration  and  indicates  that  the  pro- 
gram performs  computations  taking  the  strip  conductor  to  have  a metallization  thickness 
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of  6.35  pm.  The  next  line  is  skipped,  and  the  third  line  lists  the  preceding  six  input 
parameters  for  the  configuration  to  be  analyzed.  After  skipping  another  line,  the  follow- 
ing quantities  are  tabulated  from  left  to  right  with  identifying  labels  printed  directly  above 
computed  values. 

ZO(OHMS)  = characteristic  impedance  of  the  transmission  line  in  ohms; 

V(M/SEC)  = phase  velocity  in  meters  per  second; 

ALPHAC(DB/M)/SQRT(F)  = conductor-loss  coefficient  in  dB  per  meter  divided  by 

the  square  root  of  the  frequency  in  hertz; 

ALPHAD(DB/M)/F  = dielectric-loss  coefficient  in  dB  per  meter  divided  by  the 
frequency  in  hertz. 

For  program  executions  in  which  more  than  one  configuration  is  to  be  analyzed 
(NSETS  > 1),  an  output  block  similar  to  that  described  is  printed  fo<-  each  configuration. 
Four  blank  lines  separate  the  output  information  blocks  for  each  configuration. 


Program  CPW  (Coplanar  Waveguide) 

Program  CPW  is  written  in  Fortran  IV  and  is  compatible  with  the  CDC  6700  com- 
puter system.  With  minor  modifications  this  program  should  also  be  suitable  for  use  with 
the  IBM  360  or  370  systems  and  with  any  other  computer  having  comparable  storage 
capability  and  compatibility  with  the  Fortran  IV  language.  The  storage  necessary  for  ex- 
ecuting this  program  on  the  CDC  6700  computer  at  the  Naval  Ship  Research  and  Develop- 
ment Center  (NSRDC)  is  86,143  words.  The  length  of  time  necessary  for  program  execu- 
tion is  approximately  15  minutes  per  configuration. 


Input  Information 

Each  configuration  of  coplanar  waveguide  to  be  analyzed  by  program  CPW  is  char- 
acterized by  specifying  the  cross  section  in  terms  of  the  geometry  shown  in  Fig.  5.  Ac- 
cordingly the  first  data  card  for  any  execution  of  CPW  contains  the  quantity  NSETS, 
specified  in  the  110  format.  This  is  an  integer  specifying  the  number  of  configurations 
to  be  analyzed. 

Each  configuration  included  in  NSETS  is  specified  on  a separate  data  card  located 
consecutively  behind  the  first  data  card.  Each  card  lists  the  following  quantities  punched 
according  to  the  format  F10.6: 

W = width  of  the  conducting  strip  in  mils  (1  mil  = 25.4  pm); 

S = spacing  between  the  strip  edge  and  grounded  conductors  in  mils; 

H = thickness  of  the  substrate  in  mils; 
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L = substrate  width  in  inches,  1.0  inch  (2.54  cm)  being  adequate  to  elim- 
inate effects  of  finite  width; 

ER  = relative  permittivity  of  the  substrate  material; 

TAND1  = loss  tangent  of  the  substrate  dielectric  (X  104); 

SIGMA  = dc  conductivity  of  the  metallization  in  mhos  per  meter  (X  10”7). 

The  number  of  cards  bc'*:"'d  the  first  data  card  should  equal  the  value  read  in  for  NSETS. 


Output  Information 

The  output  information  provided  by  CPW  is  similar  in  format  to  that  for  program  MS , 
The  first  line  of  output  for  each  configuration  indicates  the  metallization  thickness.  After 
a line  is  skipped,  the  input  parameters  are  listed  in  a line  of  output.  After  another  line  is 
skipped,  the  quantities  ZO(OHMS),  V(M/SEC),  ALPHAC(DB/M)/SQRT(F),  and  ALPHAD 
(DB/M)/F  are  printed  in  a manner  similar  to  that  described  for  program  MS.  For  program 
executions  in  which  more  than  one  configuration  is  to  be  analyzed  (NSETS  > 1),  an  out- 
put block  similar  to  that  already  described  is  printed  for  each  configuration. 


Program  IM  (Inverted  Microstrip) 

Program  IM  is  written  in  Fortran  IV  and  is  compatible  with  the  CDC  6700  computer 
system.  With  minor  modifications  this  program  should  also  be  suitable  for  use  with  the 
IBM  360  or  370  systems  and  with  any  other  computers  having  comparable  storage  capa- 
bility and  compatibility  with  the  Fortran  IV  language.  The  storage  necessary  for  executing 
this  program  on  the  CDC  6700  computer  at  the  Naval  Ship  Research  and  Development 
Center  is  78,698  words.  The  length  of  time  necessary  for  program  execution  is  approxi- 
mately 12  minutes  per  configuration. 


Input  Information 

Each  configuration  of  inverted  microstrip  to  be  analyzed  by  program  IM  is  character- 
ized by  specifying  the  cross  section  in  terms  of  the  geometry  shown  in  Fig.  8.  Accordingly 
the  first  data  card  for  any  execution  of  IM  contains  the  quantity  NSETS,  specified  in  the 
110  format.  This  is  an  integer  specifying  the  number  of  configurations  to  be  analyzed. 

Each  configuration  included  in  NSETS  is  specified  on  a separate  data  card  located 
consecutively  behind  the  first  data  card.  Each  card  requires  the  following  quantities 
punched  according  to  the  format  F10.6: 

W = width  of  the  conducting  strip  in  mils  (1  mil  = 25.4  pm); 

H = height  of  the  strip  above  the  ground  plane  in  mils; 
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T = substrate  thickness  in  mils; 

I,  = substrate  width  in  inches,  1.0  inch  (2.54  cm)  being  adequate  to  elim- 
inate effects  due  to  finite  width; 

ER  - relative  permittivity  of  the  substrate  material; 

TAND1  = loss  tangent  of  the  substrate  dielectric  (X  104); 

SIGMA  = dc  conductivity  of  the  metallization  in  mhos  per  meter  (X  10"7). 

The  number  of  cards  behind  the  first  data  card  should  equal  the  value  read  in  for  NSETS. 


Output  Information 

The  output  information  provided  by  program  IM  appears  in  the  same  format  as  that 
described  for  programs  MS  and  CPW. 


Program  TIM  (Trapped  Inverted  Microstrip) 

Program  TIM  is  written  in  Fortnm  IV  and  is  compatible  with  the  CDC  3800  computer 
system.  With  minor  modifications  this  program  should  be  suitable  for  use  with  the  CDC 
6000-series  computer  systems,  the  IBM  360  or  370  systems,  and  any  other  computer  having 
comparable  storage  capability  and  compatibility  with  the  Fortran  IV  language.  The  core 
storage  necessary  for  executing  the  program  on  the  CDC  3800  system  at  the  NRL  Research 
Computation  Center  is  42,000  words.  Although  the  total  core  storage  of  this  computer  is 
now  75,536  words,  the  maximum  storage  available  for  a single  array  (without  resorting  to 
special  array  — handling  techniques)  is  32,768  words.  It  is  this  system  feature  that  con- 
strains the  largest  array  in  TIM  to  be  181  X 181  or  32,761  words.  To  properly  load  this 
program  into  the  computer,  a procedure  similar  to  the  one  described  for  program  MS  is 
used. 

The  length  of  time  necessary  for  program  execution  is  approximately  7.8  minutes  per 
configuration. 


Input  Information 

Each  configuration  of  trapped  inverted  microstrip  to  be  analyzed  by  program  TIM  is 
characterized  by  specifying  the  cross  section  in  terms  of  the  geometry  shown  in  Fig.  12. 
Accordingly  the  first  data  card  for  any  execution  of  TIM  contains  the  quantity  NSETS, 
specified  in  the  110  format.  This  is  an  integer  specifying  the  number  of  configurations  to 
be  analyzed. 

Each  configuration  included  in  N5ETS  is  specified  on  a separate  data  card  located 
consecutively  behind  the  first  data  card.  Each  card  lists  the  following  quantities,  punched 
according  to  the  format  F10.6: 
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IV  = width  of  the  conducting  strip  in  mils  (1  mil  = 25.4  jim); 

S = spacing  between  the  strip  edge  and  the  channel  upper  comer  in  mils; 

H - height  of  the  strip  above  the  channel  bottom  in  mils; 

T = substrate  thickness  in  mils; 

L = substrate  (and  support  block)  width  in  inches,  1.0  inch  (2.54  cm) 
being  adequate  to  eliminate  effects  of  substrate  and  support  block 
finite  widths; 

ER  = relative  permittivity  of  the  substrate  material; 

TAND1  = loss  tangent  of  the  substrate  dielectric  (X  104); 

SIGMA  = dc  conductivity  of  the  metallization  specified  in  mhos  per  meter  (X  10~7). 
The  number  of  cards  behind  the  first  data  card  should  equal  the  value  read  in  for  NSETS. 

Output  Information 

The  output  information  provided  by  program  TIM  appears  in  the  same  format  as 
that  described  for  programs  MS  and  CPW. 


Program  LEMDOC  (Lossy  Edge-Coupled-Microstrip 

Dielectric-Overlay  Coupler) 

Program  LEMDOC  is  written  in  Fortran  IV  and  is  compatible  with  the  CDC  3800 
digital  computer.  With  minor  modifications  this  program  should  also  be  compatible  with 
the  CDC  6000-series  computers,  the  IBM  360  or  370  systems,  or  any  other  computer 
which  has  similar  storage  capability  and  compatibility  with  the  Fortran  IV  language.  The 
core  storage  necessary  for  executing  this  program  on  the  CDC  3800  system  at  the  NRL 
Research  Computation  Center  is  about  42,000  words.  Although  the  total  storage  of  this 
computer  is  now  75,536  words,  the  maximum  storage  available  for  a single  array  (without 
resorting  to  special  array-handling  techniques)  is  32,768  words.  It  is  this  system  require- 
ment that  constrains  the  largest  array  in  LEMDOC  to  be  181  X 181  or  32,761.  Special 
techniques  required  for  the  loading  of  this  program  are  similar  to  those  described  for 
program  MS. 

The  length  of  time  necessary  for  program  execution  is  approximately  15  minutes 
per  configuration. 


Input  Information 

Each  configuration  of  edge-coupled  microstrip  with  a dielectric  overlay  to  be  analyzed 
by  program  LEMDOC  can  be  characterized  by  specifying  the  cross  section  in  terms  of  the 
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geometry  shown  in  Fig.  12.  Accordingly  the  first  data  card  for  any  execution  of  LFMDOC 
contains  the  quantity  NSETS,  an  integer  specifying  the  number  of  configurations  to  be  ana- 
lyzed. NSETS  is  punched  on  the  first  data  cards  in  the  110  format. 

Each  configuration  included  in  NSETS  is  specified  on  a separate  data  card  located 
consecutively  behind  the  first  data  card.  Each  card  requires  the  following  quantities, 
punched  according  to  the  F10.6  format: 


W ~ conducting  strip  width  in  mils  (1  mil  = 25.4  pm); 

S = spacing  between  the  two  conducting  strips  in  mils; 

T = thickness  of  the  dielectric  overlay  in  mils; 

H = height  of  the  strips  above  the  ground  plane  in  mils; 

L - width  of  the  supporting  substrate  in  inches,  L - 1.0  inch  (2.54  cm)  being 
adequate  to  eliminate  effects  due  to  finite  width; 


ER  = relative  permittivity  of  the  substrate  and  overlay  dielectric  material; 
TAND1  = loss  tangent  of  the  substrate  and  overlay  dielectric  (X  104); 

SIGMA  = dc  conductivity  of  the  metallization  in  mhos  per  meter  (X  10"7). 


The  number  of  cards  behind  the  first  data  card  should  equal  the  value  of  NSETS. 


Output  Information 

An  output  information  block  is  printed  for  each  consecutive  input  data  block.  The 
first  line  of  output  stipulates  the  thickness  of  the  metallization.  After  a line  is  skipped, 
the  input  parameters  are  listed  from  left  to  right.  After  another  line  is  skipped,  the  fol- 
lowing output  quantities  are  tabulated  directly  under  an  appropriate  label:' 

C(DB)  = computed  value  of  midband  coupling  in  dB; 

ZO(OHMS)  = coupled  line  characteristic  impedance  in  ohms; 

ZOO(OHMS)  = odd-mode  impedance  in  ohms; 

ZOE(OHMP1  = even-mode  impedance  in  ohms; 

VO(M/SEC)  = odd-mode  phase  velocity  in  meters  per  second; 

VE(M/SEC)  = even-mode  phase  velocity  in  meters/second; 

VAVG(M/SEC)  = average  of  VO(M/SEC)  and  VE(M/SEC)  in  meters  per  second; 

ALFACO(DB/M)/SQRT(F)  = odd-mode  conductor-loss  coefficient  in  dB  per  meter 

divided  by  the  square  root  of  the  frequency  in  hertz; 
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ALFACE(DB/M)/SQRT(F)  = the  even-mode  conductor-loss  coefficient  in  dB  per 

meter  divided  by  the  square  root  of  the  frequency 
in  hertz; 


1 

f 


ALFADO(DB/M)/F  = odd-mode  dielectric-loss  coefficient  in  dB  per  meter  divided 
by  the  frequency  in  hertz; 


ALFADE(DB/M)/F  = even-mode  dielectric-loss  coefficient  in  dB  per  meter  divided 
by  the  frequency  in  hertz. 
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Appendix  A 

LISTING  OF  PROGRAM  MS 

program  ms 
real  l 

READ  1»  NSETS 

1 FORMAT  < 1 10 ) 

DO  2 MN=1 »NSETS 
READ  3*  W*H*L»ER*TAND1»SIGMA 
3 FORMAT ( 6F10»6) 

TAND1=TAND1*1.E-0A 
SIGMA=SlGMA*l.E+07 
PRINT  8 

ft  FORMAT ( 1 H *//// 1H  **STRIP  Mf.TALlZAT  ION  THICKNESS  = 250  MICROINCHES 
1*) 

PRINT  A » W»H*L»ER*TAND1»SIGMA 
A FORMAT! lH0t2HW=»F5.1*3X.2HH=*F5.1 »3am2HL=* 
1F10,6*3X*3HER=»F10.6*3X.6HIAND1=*F1u.6*3a*6H3IGMA=*F1u.6) 

FACTOR= (1.0E-05)*(1./o.3937uu) 

W=W*FACTOR 

H=H*FACT0R 

L=L*FACT0R*<l.E+03) 

E1=ER 

E2=ER+0.01*ER 

CALL  MSF(W»H*L*E1»SIGMA»C*PDDSQ) 

pddsqf=pddsq 

CALL  MSF(W,H»L»E2*SIGMA*CP*FDDSQ) 

CALL  MSE(W»H»L.CE) 

A=C#CE 

B=CE/C 

ZO=l./( ! 3.0E+08 ) #SQRT ( A ) ) 

V=3.0E+08*SQRT(B) 

ARG1=B 

ARG2=CE/CP 

ER1=1./ARG1 

ER2=1./AR62 

FACTC=10./2.3 

FACTD=27.3/3.E+08 

AC=FAC1 C*PDDSQF*ERl*Zu 

AD=FACTD*(F 1/SORT (FR1) )*( ( ER?“ER1 ) / ( E2-E1 ) ) * I AND1 
PRINT  5 

ft  FORMAT! 1hO»1X.8HZO (OHMS) > 5X »8HV( M/SEC ) *5X. 

1?0HALPHAC(DR/M)/SQRT(F) *5X*1AHALPHAD(DB/M> / F ) 

PRINT  6*  Zu  »V  * AC » AD 

6 FORMAT  ! 1H  ♦2X.F6.2»AX»E10.3»8X*E1'J.3»11X»E10.3) 

2 CONTINUE 
STOP 
END 

SUBROUTINE  MSF(W»H»L*ER»SIGMA»CHRG»PDDSQF) 

DIMENSION  N(A)*X(91»A)*Y(91*A) * ALPHA! 91 >A) »BETA1(91»A) » 

1BETA2  < 91 » A ) » GAMMA! 91 >A) *CH<91 >A) >SCH( A) »TSCD( 91 »A*A)»A(181»181)» 
2A1(13A»13A) 

REAL  L 

common/help/a»ai 

EQUIVALENCE ( A * A 1 ) 

NO=A 
N ( 1 )=A3 
N 1 2 ) =91 
N ( 3 ) =25 
N ( A ) =25 

E0=l./(A.*3.1A159*2.99776*2.99776E+09) 

AAA=3.1A159*A.E-07*3.1A159/SI6MA 
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RDSQF=SQRT(AAA> 

FACTN=4. *3. 14159*3. 14159*RDSQF*B.85E-12 
FACTD=4. *3. 14159*1. E-07 
FACTL=FACTN/FACTD 
FACTOR=«1.0E~05)*(1./<J.  393700) 
DELW=W/20. 

E=ER*E0 
X ( 1 * 1 ) =-W/2 • 

Y ( 1 ♦ 1 )=H 
DO  26  1=1*20 
IP1=1+1 

X( IP1»1)=X(I*1)+DELW 

26  Y( IP1 *1 ) =H 

X ( 22  *1 ) =X(2 1 »1 ) 

Y(22*D=0.25*FACTOR+H 
DO  27  1=22*41 
IP1=1+1 

X( I PI  * 1 ) =X ( I*1)-DELW 

27  Y(IP1*1)=Y(I»1) 

X ( 43 ♦ 1 ) =X( 1 * 1 ) 

Y 1 43  * 1 ) =Y  < 1 * 1 ) 

DELL=L/44. 

X<l*2)=~L/2. 

Y( 1 *2 ) =0. 

DO  20  1=1*44 
IP1=I+1 

X( I PI *2 ) =X ( I *2)+0ELL 

28  Y ( IP1 .2 ) =0. 

X ( 46  *2 ) =X ( 45  *2 ) 

Y ( 46  *2 ) =-0»25*FACTOR 
DO  60  1=46*89 
IP1=I+1 

X(  IP1 *2 ) =X( I .2 )-DELL 

60  Y ( I PI *2 ) =Y ( I .2) 

X ( 91 .2 ) =X ( 1 *2 ) 

Y(91*2)=Y(1*2) 

DELH=H/10. 

DELD=(L/2.-W/2. > / 14 . 

Xtl,3)=Xtl,2) 

Y(  1 *3 ) =Y ( 1 *2 ) 

X<1 *4)=-X(l *3) 

YC 1 *4 ) =Y ( 1 .3 ) 

DO  61  1=1*10 
IP1=I+1 

X( I PI .3 ) =X ( I .3) 

Y ( IP1 *3 ) =Y ( I *3 ) +DELH 
X ( I PI *4 ) =X ( I *4  ) 

61  Y( I PI *4 ) =Y ( I PI >3 ) 

DO  62  1=11*24 
IP1=I+1 

X( IP1*3)=X( I *3 ) +DELD 
Y ( I PI *3 ) =Y ( I *3) 

X( I PI *4 ) =X ( I *4 ) -DELD 

62  Y< I PI *4 ) =Y( IP1 *3 ) 

DO  63  1=1.42 
ALPHA ( I » 1 > = 1 . 

BETAK  I »1>=0. 

BETA2 (1*1) =0. 

63  GAMMA ( I *1»=1. 
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DO  64  1=1 .90 
ALPHA! I » 2 > =1 • 

BETA1 ( I *2)=0. 

BETA2!I*2)=0. 

64  GAMMA! I .2)=0. 

DO  65  1=1.24 
ALPHA! I »3>=0. 

BETA1 ! I »3 ) =E 
BETA2! I >3)=EO 
GAMMA! I »3)=0. 

ALPHA! I »4>=0. 

BETAH  1 .4 ) =EO 
BETA2< I »4 ) =E 

65  GAMMA! I .41=0. 

XMIN=0. 

XMAX=0. 

YMIN=0. 

YMAX=0. 

NX=0 

NY=0 

IDIM=91 

R=1.0E+05 

NAXDIM=181 

NAYDIM=181 

CALL  LPLACF !NO.N.X.Y. ALPHA. BE !A1 »BE I A2 .GAMMA »CH »aCH » IDIM .R » I iCD* 

ixmin.xmax.nx.ymin.ymax.ny.naxdim.naydim) 

DELF=W/20. 

DELE=0.25*FACT0R 

FFl=6.2831852*t*DELF 

FF2=6.2831852*E0*DELF 

FF4=6.2831852*E0*DELE 

CHRG=0. 

DO  50  1=1.20 

50  CHRG=CHRG+FF1*CH( I .1) 

DO  51  1=22.41 

51  CHRG=CHRG+FF2*CH( I »1 ) 

CHRG=CHRG+FF4*<CH<  21 »1 )+CH(42 * 1 > ) 

DO  1064  J=1  .NO 

I F ( ALPHA (l.J)«NE«l»)  GO  TO  64 
NJ1=N ( J > -1 
DO  1065  1 = 1 . NJ1 
1P1=I+1 

SEGX=X!IP1,J)-X(I»J) 

SEGY=Y( IPl.J)-Y(I.J) 

SEGARG=SFGX*SEGX+SEGY*SEGY 

SEGLEN=SQRT(SEGARG) 

1065  SUM=SUM+CH( I »J)*CH< I »J)*SEGLEN 
1064  CONTINUE 

PDDSQF=FACTL*SUM 
54  CONTINUE 
RETURN 
END 

SUBROUTINE  LPLACF < NO. N i X.Y. ALPHA. BETA1.BE I A2 ’ GAMMA > CH >SCH * IDIM.R. 
1 TSCD>XmIN»XMAX»NX.YMIN.YmAX  >nY .nAXDIm.nAYDIm) 

DIMENSION  X!  IDIM  .NO)  . Yl.IDIM.NO  ) > ALPHA  l I D IM  .NO  ) »faE  I A1 ! I DIM  .NO)  > 
1BETA2< ID IM.NO ) .GAMMA! IDIM.NO) »CH! I DIM. NO) .TSCD* IDIM.N0.4) >N«NO) . 
2A( 181 .181) .SCH(NO) »B ( 180 ) »Al < 134 » 134 ) 

C0MM0N/HELP/A.A1 
FOU I VALENCE ( A » A 1 ) 
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PI=3. 1415926 
RR=R*R 
DO  1 L=1 *NO 
NN=N(L)-1 

00  1 I = 1 *NN 

XI  = X( I+1»L)-X( I»L) 

Y1=Y(I+1,L)-Y( I»L) 

TSCD (I ♦ L ♦ 1 ) =ATAN2 ( Y I . X I ) 

TSCD< 1*L*2)=SIN<TSC0« I *L»1) ) 

TSCD ( I *L.3)=COS(TSCDl  I-*L»1) ) 

1 TSCDU»L*4)=SQRT(XI*XI+Yl*YI) 

JJJ  = 0 

DO  4 LJ=l*NO 
NJ=N ( LJ ) -1 
JAJ=JJJ 
J J J= JJJ+NJ 
DO  4 J=1 *NJ 
JJ= JAJ+ J 

X J= ( X ( J » L J ) +X ( J+l . L J ) 1 / 2 • 

YJ=(  Y<  J»l_J)+Y(  J+l  *LJ>  )/2. 

1 I 1=0 

DO  4 L 1=1 »NO 
N I =N ( L I )-l 
111=111 +N I 
I A I = 1 1 1 -N I 
DO  4 1 = 1 * N I 
I I = I AI+  1 

I F ( I I.EO.JJ)  GO  TO  3 
X1  = XJ-X ( I »L  I ) 

X2=XJ-X< 1+1.L1 ) 

Yl  = Yj-Y ( I » L I ) 

Y2= YJ-Y ( 1+1 » L I ) 

R1=X1*X1+Y1*Y1 

R2=X2#X2+Y2*Y2 

51  = 0. 

52  = 0. 

YT  = Y( I »L1 )*X2-Y( 1+1 »LI ) #Xl+Yj* ( X ( 1 + 1 »Ll )-X( I * L I ) ) 

XT=XI*X2+Y1*Y2 
TETA=ATAN2( YT.XT) 

! F ( ALPHA ( J ? L J ! t EQs 0 * ) GO  lo  2 

Sl  = TSCD<  I »LI  »4)*<  l.-u.5*ALOGtR2/'RK)  )+u.  5*ALOG  IK2/K1 ) * l Al*  I SCO  1 1 ♦ L I 
1 >3<'+Yi*TSCD<  I ’Ll  »2>  )+TETA*5Xl*lbCD(  I * L I * 2 > -Y  1*  I ^CD  i I *L  I >3  ) > 

? TeTa1  = TScD(J*LJ>1)-TSCD( I *L I * 1 > 

S2=0.5*SIN(TETA1)*A!.OG«R2/Rl)+C0ii  iEiA1)*iEH 

S3=-S2 

GO  TO  4 

* S1=TSCD< 1 »LI >4)*< l.-ALOG(TSCD< I »LI »4 ) /2./R > ) 

S2=-PI 

S3=-P1 

4 A ( J J * 1 1 ) = ALPHA (J»LJ)*S1+BETA1(J>LJ)#S2+BETA2(J*LJ)*S3 
M=0 

DO  5 1=1 »NO 

5 M=M+N ( L ) -1 
JJJ  = 0 

DO  6 L= 1 »NO 
NN=N ( L ) -1 
JAJ= JJJ 
JJJ=JJJ+NN 
DO  6 J= 1 »NN 
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JJ=JAJ+J 

6 81 JJ)=GAMMA( J»L) 

CALL  ARRAY(2*M*M»NAXDIM»NAYDIM*A»A) 
call  simq<a»b*m»ks> 

IF  IKS.NE.O)  PRINT  low 
10o  FORMAT « 1H0*18HSYSTEM  lb  SINGULAR) 

JJJ=0 

DO  7 L=1 »NO 

NN=N(L)-1 

JAJ=JJJ 

JJJ=JJJ+NN 

DO  7 J= 1 »NN 

JJ=JAJ+J 

7 CH( J»L>=B< JJ) 

DO  8 L=1 »NO 

NN=N(L)-1 

SCH(U=0. 

DO  8 1 = 1,  NN 

8 SCHI L > =SCH( L ) +TSCD I 1 ,L»4)*CH( I »L) 

IF  tNX-1 ) 17,9,10 

9 DX=0. 

GO  TO  11 

10  OX=(XMAX-XMIN)/FLOAT(NX-1) 

U IFINY-1 ) 17,12,13 

12  DY=0. 

GO  TO  14 

18  DY=(YMAX-YMIN) /FLOAT (NY-1) 

14  DO  16  I 1 = 1, NX 

XJ  = XMIN+FLOAT I I 1-1 > *DX 
DO  16  JJ=1’,NY 
YJ=YMIN+FLOAT< JJ-1 )*DY 
A ( II,JJ)=0. 

DO  16  L I =1  »N0 
NN=N ( L I )-l 
DO  16  1 = 1, NN 
Xl  = XJ-X( I , L I ) 

X2  = XJ-X( 1 + 1 »LI  ) 

Y 1 =YJ-Y ( I ,L  I ) 

Y2  = YJ-Y( 1 + 1 , L I ) 

R 1 = X1*X  l+Yl  -&Y  ] 

R2=X2*X2+Y2*Y2 

I F < <R1.EO.O.).OR.(R2.EO.O,) ) GO  TO  15 

YT  = Y ( I ,LI )*X?-Y( 1+1 ,LI ) *Xl+YJ* (XI  I + 1 »L J ) -X ( I , L I ) ) 

XT=X1*X2+Y1#Y2 
T£Ta=ATaN2( YT,XT) 

Sl=TSCD( I »LI*4)*(1.-0.5#ALOG(R2/RR) )+0.5*ALOG(R2/Rl)*(Xl*TbCDt I*LI 
1*8  >+Yl*TSCD< I »LI ,2) )+T£TA*( Xl*TSCD( I *LI »2>-Yl*TSCD( I ,LI ,3) ) 

GO  TO  16 

15  Sl=TSCD< I »L I *4)*( l.-0.5*AL0G( (R1+R2 ) /RR) ) 

16  AI 1 1 ,JJ)=A( 1 1 »JJ)+S1*CH< I ,LI) 

17  RETURN 
END 

SUBROUTINE  MSE(W»H»L*CE) 

DIMENSION  N I 2 ) *X ( 91 » 2 ) * Y ( 91 , 2 ) >ALPHA(91*2) ‘BETA! I 91 , 2 ) , 

1RETA2<91 *2) > GAMMA  I 9 1 *2) *CH<91 ,2  > *SCH<  2) * ISCD191 *2  *4 ) >A> 181 >181) » 
2A1 ( 134,134) 

REAL  L 

COMMON/HELP/A, Al 
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EQUIVALENCE! A* Al) 

N0=2 
N ( 1 ) =43 
N ( 2 ) =91 

EO  = l •/ ( 4 »*3« 14159*2*99776*2  *99 776E+Q9 ) 

FACTORS  1.UE-05HM1./0. 393700) 

DELW=W/20. 

X( l»l)a-W/2. 

Y ( 1 » 1 ) =H 
DO  26  1=1*20 
IPl=i+l 

X! I PI » 1 )=X( I » 1 ) +DELW 

26  Y ( I PI  * 1 ) =H 
X(22*1)=X(21*1 ) 

y (22*1 )=0.25*FACTOR+H 
00  27  1=22*41 
IP1=I+1 

X(  I P I * 1 ) =X ( I.ll-OELW 

27  Y(1P1*1)=Y(I*1) 

X ( 43  » 1 ) =X( 1 » 1 ) 

Y(43*1)=Y(1 * 1 ) 

DELL=L/44. 

X ( 1 * 2 ) =“L/2  * 

Y ( 1 1 2 ) =*J  • 

00  28  1=1*44 
IP1=I+1 

X(  IP1»2)=X( I *2! +DELL 

28  Y ( I PI »2 ) =Q • 

X ( 46  *2 ) =X ( 45  *2 ) 

Y(46»2)=-0.25*FACT0R 
DO  60  I =46  >89 
IP1=I+X 

X< I PI »2 ) =X ( 1*2 l-DELL 
60  Y( I PI *2 i =Y( 1*2) 

X ( 91 *2 ) =X( 1 *2) 

Y ( 91 *2 ) =Y  C 1 >2) 

DO  63  1=1*42 
ALPHAt I * 1 » = 1 . 

BETAK  I >1)=0. 

BET  A2  l I*D  =0. 

63  GAMMA ( I » 1 ) =1  • 

DO  64  I=1»9U 
ALPHA ( I *2 ) = 1 - 
BETAK  I *2)=0. 

BETA2 ( I *2  » =0. 

64  GAMMA < I >2)=0. 

XMIN=G. 

XMAX=0« 

YMIN=0. 

YMAX=0. 

NX=0 

NY=0 
ID IM=91 
R=1.0E+05 
NAXDIM= 134 

CALL ILPLACF(NO*N*X»Y*ALPHA*BETa1»(3EtA2 ’GAMMA *CH*SCH* ID IM»R*TSCD 

1XMIN*XMAX,NX*YMIN*YMAX*NY*NAXDIM.NAYDIM) 

CE=6.2831.852*Eu*SCH(l) 
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RETURN 

END 

SUBROUTINE  LPLACE I NO*N »X*Y » ALPHA* BE T A1 » BET A2 * GAMMA »CH*SCH* IDIM*R * 

l tscd*xmin*xmax,nx,ymin.ymax.ny,naxdim*naydim) 

DIMENSION  XI IDIM*NO) »Y( IDIM.NO) .ALPHA ( IDIMtNO) »BETaK IDIM.NO) » 
1BETa?UDIM.NO)  . GAMMA ( IDIM.NO)  »CH(  IDIM»NO)  *TSCD<  ip IM  *NO .4  ) *N  ( NO)  * 
2A< 181*181) »SCH( NO) *B l 132 ) .Alt  134. 134) 
common/help/a*ai 

EQUIVALENCE (A* A1 ) 

PI=3. 1415926 
RR=R*R 
DO  1 L=l*NO 
NN=N ( L > — 1 
DO  1 1 = 1. NN 
XI=X(I+1*L)-X( I»L> 

Y I =Y ( 1+1 *L ) -Y ( I *L ) 

TSCD ( I *L . 1 ) =ATAN2 ( Y I »X I) 

T SCD I I*L'.2)=SIN(TSCDl  I *L  * 1 ) > 

TSCDl I*L.3)=COS(TSCDl I »L*D  ) 

1 TSCDlI .L.4)=SQRT(XI*XI+Y1*YI) 

JJJ=0 

DO  4 1 1 = 1 *NO 

NJ=N(LJ)-1 
JAJ= JJJ 
J J J= JJJ+NJ 
DO  4 J= 1 »NJ 
J J=JAJ+J 

XJ=(X(J»LJ)+X(J+1*LJ) )/2. 

YJr=(  Y(J»LJ)+Y(  J+l.LJ)  )/2. 

I 11  = 0 

DO  4 L I =1 »N0 
N I =N ( L I ) -1 
1 1 1 = 1 1 1 +N I 
I A I = 1 1 I -N I 
DO  4 1=1 »NI 
I I = I AI+ I 

I F ( I I.EQ.JJ)  GO  TO  3 
X1=XJ-X ( I »L I ) 

X2=XJ-X ( 1+1 *L I ) 

Y 1 = YJ-Y  t I »L I ) 

Y2=YJ-Y< I + l.LI  ) 

R1=X1*X1+Y1*Y1 
R2=X2*X2+Y2*Y2 
S1  = 0. 

S2=0. 

YT=Y< I *LI )*X2-Y( 1+1 »LI )#Xl+Yj# (XI 1+1 *LI )-X C I *LI ) ) 

XT=X1*X2+Y1#Y2 

TETA=ATAN2(YT.XT) 

+F ( ALPHA ( J *1 I ) »EQ«  0 * ) GO  TO  2 

Sl  = TSCD<  I *LI >4)*<  l.-0.5*AL0G(R2/RR>  ) +0, 5*AL0G  l R2/R1  ) M Xi*TSCD l I *LI 
l»3)+Yl*TSCDt I *LI *2) >+TETA*lXl*TSCD( I’Ll .2>-Yl*TSCD< I * L I * 3 ) > 

2 TETA1  = TSCD(J.LJ*1)-TSCD( I *L I » 1 ) 
S2=0,5«SiN(TfTa1)*AL0G(R2/R1)+C0S(TETa1)*TETA 
S3=-S2 

GO  TO  4 

■»  Sl  = TSCD(  I »LI  *4)*(  l.-ALOG<  TSCDl  I *LI *4 ) /2 . /R > > 

S2=-P I 
S3=-P I 

4 At JJ. 1 1 )=ALPHA( J*LJ)«S1+BETA1< J*LJ)#S2+BETA2l J»LJ)*S3 
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M=0 

DO  5 L=1 *NO 

5 M=M+N ( L > -1 
JJJ=0 

DO  6 L= 1 »NO 
NN=N ( L ) -I 
JA J=JJ J 
JJJ= JJJ+NN 
DO  6 J-l.NN 
JJ= JAJ+ J 

6 81 JJ)=GAMMA( J*L) 

CALL  ARRAY (2*M«M*NAXDIM*NAYDIM*A»A) 

CALL  SIMQ(A*B*M*KS) 

IF  IKS.NE.O)  PRINT  100 

lOo  FORf'ATI  1H0.18HSYSTEM  IS  SINGULAR) 

JJJ=0 

DO  7 L=l,NO 
NN=N ( L ) -1 
JAJ= JJJ 
JJJ= JJJ+NN 
DO  7 J=1,NN 
J J= JAJ+ J 

7 CH( J»L)=B( JJ) 

DO  8 L=l,NO 
NN=N(L)-1 
SCH(L)=0. 

DO  8 1=1, NN 

8 SCH ( L ) =SCH ( L ) +TSCD ( I »L  *4 ) *CH ( I *L ) 

I F I NX-1 ) 17*9,10 

9 DX=0„ 

GO  TO  11 

10  DX=(XMAX-XMIN)/FLOAT(NX-l> 

11  I F (NY-1 >17*12*13 

12  DY=0. 

GO  TO  14 

U f)Y=  ( YMAX-YMIN)  / FLOAT  ( NY- 1 ) 

14  DO  16  I 1 = 1. NX 
XJ  = XMIN+FLOAT ( I 1-1 )*DX 
DO  16  J J= 1 *NY 
vj=YMIN+FLOAT< JJ-1 )*OY 
A ( I I ,JJ)=0, 

DO  16  L I = 1 *NO 
NN=N (LI ) -1 
DO  16  I =1 »NN 
X1  = XJ-X( I »L  I ) 

X2=XJ-X( 1+1 ,LI ) 

Y 1 = YJ-Y { I *L  I ) 

Y2=YJ-Y ( 1+1 . L I ) 

R1=X1*XI+Y1*Y1 
R2=X2*X2+Y2*Y2 
IF((R1.EQ.O,).OR.(R2*EQ.O,))  GO  TO  15 

YT=Y( I *Ll )*X2-Y( 1+1 *LI )*Xl+YJ*(X< 1+1 *LI ) -X ( I *L I ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT.XT) 

S 1 = TSCD I I *LI »4)*< l.-0.5*ALOG(R2/RR) ) +u. 5*AlOG ( R2/R1 ) * ( Xl*T6CD ( I * L I 
1 »3)+Yl*TSC0< I «L I * 2 ) )+TFTA*(Xl*TScD( I * L I *2  )-Yl*Ti>CD< I *LI *S> ) 

GO  TO  16 

15  Sl=TSCD( I »LI »4)*( l.-o,5*ALOG( (Rl+R2)/RR) ) 
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16  A( 1 1 *JJ)=A( 1 1 » JJ) +Sl#CH( I * L I ) 

17  RETURN 
END 

SUBROUTINE  ARRAY  (MODE*I »J»N*M*S.D) 

dimension  S ( 1 ) *D( D 

NI=N-I 

IF(MODE-l)  100*100*120 
100  IJ=I*J+1 
NM=N*J+1 
DO  110  K=1 ♦ J 
NM=NM-NI 
DO  110  L=1 * I 
IJ=IJ-1 
NM=NM-1 

110  D(NM»=S( IJ) 

60  TO  140 
120  IJ=0 
NM=0 

DO  130  K=1 * J 
DO  125  L = 1 » I 
IJ=IJ+1 
NM=NM+1 

125  S( IJ)=D(NM) 

130  NM=NM+NI 
140  RETURN 
END 

SUBROUTINE  SIMQ(A»B*N»KS) 

DIMENSION  All J »0 C 1) 

TOL=0.0 

KS=0 

JJ=-N 

DO  65  J=1»N 
JY=J+1 
J J=J J+N+l 
BIGA=0. 

I T=JJ-J 
DO  30  I = J *N 
I J=IT+I 

I F ( ABS ( B I GA ) -ABS ( A d J ) ) ) 20*30.30 
20  BIGA=A( IJ) 

IMAX= i 
30  CONTINUE 

IFIABS(BIGA)-TOL)  35.35*40 
35  KS= 1 
RETURN 

40  1 1=J+N#  < J-2 ) 

IT=IMAX-J 
DO  50  K=J.N 
11=1 1+N 
I 2= 1 1+ 1 T 
SAVE=A( II) 

A ( 1 1 ) =A ( 1 2 ) 

A ( 1 2 ) =SAVE 
50  A ( 1 1 ) =A ( II ) / B I GA 
SAVE=B( IMAX) 

B ( IMAX ) =B ( J ) 

B ( J ) =SAV£/B I GA 
IF(J-N)  55*70.55 
55  I0S=N*(J-1) 
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DO  65  IX=JY.N 
IXJ-IGS+IX 
1T=J-1X 
DO  60  JX=JY»N 
I XJX=N*  ( JX-'I ) + IX 
JJX= IXJX+IT 

60  Al  IXJX)=aUXjX)-(A«  lXJ)*At  JJXH 
66  B( I X ) =B t IX)-<B<J>*A< I X J ) 1 
70  NY=N-1 
I T =N#N 

DO  80  J=1»NY 
IA=IT-J 
IB=N-J 
IC=N 

DO  80  K=1 » J 

B( IB)=B( IB5~A( IA)*B< IC) 

I A= IA-N 
80  I C= IC-1 
RETURN 
END 
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PROGRAM  CPW! INPUT tOUTPUT) 

REAL  L 

READ  1.  NSETS 
FORMAT! I XO) 

DO  2 MN*1»NSETS 

RFAD  3 ♦ WtStHtLtERtTANDltSlGMA 

FORMAT ( 7F10.6) 

TANDl=TANDl*l.OE-0A 
SIGMA=SIGMA*1.0E>07 
PRINT  8 

FORMAT! 1H  »////lH  t*STRIP  METALiZATION  THICKNESS  = ij 0 MICROINCHES 
1*1 

PRINT  A»  WtStHtLtERtTANDltSlGMA 

A FORMAT! lH0t2HW= tF5.lt 3Xt2HS= tF5.lt 3X.2HH= tF5.lt 3X>2HL=t 
rFl0.6t3Xt3HFR=tF10o6t3Xt6HTANDl=tEl0,3t3Xt6HSlGMA=tE10.3) 

FACTOR=<  1.0E-U5)*!  l./'J, 393700) 

W=W*FACToR 

s=s«factor 

H=H*FACTOR 

L=L*FACT0R*ll.E+03) 

E1=ER 

E2=ER+0.01*ER 

CALL  CPWF ! W.S tHtLt El t SIGMA tCtPDDSQ) 

PDDSQF=PDDSO 

CALL  CPWF (WtS»HtLtE2  tSlGMA  tCP  tPDDSQ ) 

CALL  'CPWE(WtStHtLtCE) 

A=C*CE 

B=CE/C 

ZO=l./< !3.E+08)*SORT!A) ) 

V=3.E+08»SORT(B) 

ARG1=B 
ARG2=CE/CP 
ER1=1./ARG1 
ER2-1./ARG2 
FACTC=lvJ./2.3 
FACTD=27 .3/3.E+08 
AC=FACTC*PDDSQF*ERl*Zu 
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AD = FAC  T il* ! F 1 / SORT  !ERI)  )"  ! 1 ER1 1''  * E*2" 

PRINT  5 

A FORMAT! !H0tlXt8HZO(0HMS) t 5X t8HV (M/SFC ) t5Xt 
120HALPHAC<DR/M)/SQRT(F) t 5X  t 1 AHALPHAD ( DB / M ) / F ) 

PRINT  6 > ZO  t V t AC  * AD 

6 FORMAT  < 1H  t2XtF6.2»AX»El0.3«3X*F10p3tllX»E10.3) 

2 CONTINUE 
STOP 
END 

SUBROUT  I NE  CPWF (W.S.HtLtFRtSI GMA  tCHRG t PDDSQF ) 

DIMENSION  N<6) tX!59t6) tY(59t6) t ALPHA! 59 » 6 > tBETAl ( 59  1 5) t 
1BETA2 ! 59  t6 ) tGAMMA ( 59  t6 ) tCH<59»6) tSCH(6) t TSCD < 59 1 6 1 1, ) sA ( 252 1 252 ) 
REAL  L 

COMMON/HELP/Z ( 252 1 252 ) 

EQUIVALENCE ( Z t A ] 

N0=6 
N ( 1 ) =A1 
N ! 2 ) =59 
N ( 3 ) =59 
N ( A ) =2 1 
N ( 5 ) =2 1 

N ( 6 ) =51  48 
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E0  = 1. /< 4. *3. 14159*2. 99776*2. 99776E+09) 
AAA= 3 • 14 1 59*4 »E-07*3» 14 159/SIGMA 
RDSOF=SQRT(AAA) 

FACTN=4. *3. 14159*3. 14159*RDSQF*8«85E-12 
FACTD=4. *3. 14159*1. E-07 

factl=factn/facto 

FACTOR=(1.0E-05)*(1./U. 393700) 
DELW=W/19. 

E=ER*E0 

X(l»l)=~W/2» 

Y(1,1)=-0.25*FACT0R 
DO  26  1=1*19 
IP1=I+1 

X( I P 1 * 1 1 =X  C I »1)+DELW 

26  Y<  1 P 1 * 1 ) =Y (1*1) 

X(21»1)=X<20»1) 

Y(  21 *1 )=0. 

DO  27  I=21»39 
IP1=I+1 

X(  IP1  »1 1 =X( I >1 J-DELW 

27  Y i ! PI » 1 ) =Y ( I ,1) 

X(41*1>=X(1»1) 

Y ( 41  ♦ 1 ) = »'*  l 1) 

DELG=(L/2.-v.,'5.~S)/28. 

X ( 1 ,2)=-L/2. 

Y ( 1 »2 ) =0. 

X(30»2)=-(W/2.)-S 

Y(30,2)=-0.25*FACT0R 

X < 1*31  =L/2« 

Y d .3 ) =0 . 

X(30.3)  = (W/2.)+S 

Y ( 30 13 ) = -0. 25*FACTOR 
DO  28  1 = 1.28 
IP1=I+1 

I P30= 1+30 

IP30M1=IP30-1 

X(  I PI »2 ) =X  * I »2)+DELG 

Y ( I P 1 . 2 ) >"  0 . 

X< !P3u»?)=X( IP30M1 »2 )-DELG 

Y ( IP30  * 2 ) =Y ( 3u  » 2 ) 

Xi  I PI  .3 ) =X ( I .3 )-DELG 
Y ( I PI »3 ) =0. 

XI  IP3«J.3)-X(  IP30M1  »3)+DELG 

28  Y ( I P3^>  *3  ) =Y  ( 30 .3  ) 

X(59»2 )=X( 1.2) 

Y{ 59.2) =Y< 1.2) 

X ( 59 .3 ) =X( 1.3) 

Y ( 59 .3 ) =Y( 1 .3) 

DELS=S/20. 

X ( 1 »4)=-W/2*-S 
Y(1  ,4)=0. 

X ( 1 » 5 ) =W/2. +S 

Y ( 1 , 5 ) =0 , 

DO  29  1=1.20 
IP1=I+1 

X( I PI .4 ) =X ( I.4J+DELS 
Y ( I PI .4 ) =0. 

X ( IP1 .5 ) =X ( I *5 )-DELS 

29  Y ( IP1 .5) =0. 
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0ELH=H/2. 

D£LL=L/46. 

X ( 1*6)  =*•!./ 2 • 

Y ( 1*61=0, 

X ( 49 .6 ) =L/2 • 

Y!49»6) =H 
00  30  1=1*2 
IP1=I+1 
IP49=I+49 
IP49M1=IP49-1 
X( IP1»6)=X(| ,6) 

Y( IP1»6)=Y( I *6 )+DELH 
X( I P49  *6 ) =X ( IP49M1 *6) 

30  Y ( I P49 »6 ) =Y ( IP49M1 >6 )-DELH 
DO  31  1=1*46 

I P3= I +3 
IP3M1=IP3-1 

X«  I P;J * 6 ) =X ( IP3M1  *6 ) +DELL 

31  Y ( JP3»6)=H 
DO  60  1=] *40 
ALPHA! I »1)=1. 

BETAl! I *1)=0. 

BETA2!  ? *1  >=0<, 

60  GAMMA! I * 1 ) =1 . 

DO  61  I = 1 *. 58 
ALPHA! I *2 ) =1* 

ALPHA! 1 »3)=1. 

BETAl! I ,2)=0. 

BETAl! I *3i=0. 

BETA2< I *2  > =0. 

BETA2! I * 3 ) =0« 

GAMMA! 1*2) =0. 

61  GAMMA! I *3)=0. 

DO  62  1=1*20 
ALPHA! I ,4)=0. 

alpha,  i »5)=o. 

BC  Al! I *4)=E0 
Rt  -AH  I *5)=E 
BE  fA2 I I *4)=E 
BETA2! I *5)=E0 
GAMMA ! I i 4 5 =0 • 

62  GAMMA! I *5)=0. 

DO  63  1=1*50 
ALPHA! I »6)=0. 

BETAl! I *6 ) =E 
BETA2 ( I >6)=E0 

63  GAMMA! I *6)=0. 

XMIN=0. 

XMAX=0. 

YMIN=Ot 

YMAX=0. 

NX=0 

NY=0 


I D IM=59 
R=1.0E+05 
NAXD IM=252 
NAYDIM=252 


CALL  LPlACF(NO,n*X,Y» ALPHA ’BETa1»BETA2»G/MMA>CH*SCH»10IM»R*TSCD* 
1XMIN  >XMAX»NX*YMIN*YmAX  >NY  »NAXDIM*NAY0IM) 
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DELP=W/19* 

DELE=0.25*FACT0R 

FF1=6.2831852*E*DELF 

FF2=6.2831852*E0#DELE 

FF4=6.2831852*t0*0ELF 

CHRG=0. 

DO  50  1=1*19 

50  CHRG=CHRG+FF4*CH(  1*1) 

DO  51  1=21*39 

51  CHRG=CHRG+FF 1*CH (1*1) 

CHRG=CHRG+FF2*«CH( 20  * 1 ) +CH( 40 ♦ 1 > i 
SUM=0. 

DO  64  J=1 »NO 

I F < ALPHA  I 1*J)»NE.1»)  GO  TO  64 
NJl=N( J ) -1 
DO  65  1=1 *NJ1 
IP1=I+1 

SEGX=X( IP1 * J)-X( I *J) 

Se6Y=Y ( I PI  * J I -Y ( I » J ) 

SEGaRG=SeGX*SEGX+SEGY*SEGY 

SEGLEN=SQRT(SEGARG) 

65  SUM=SUM+CH(I *J)»CH< I »J)*SEGLEN 
64  CONTINUE 

PDDSQF=FACTL*SUM 
54  CONTINUE 
RETURN 
END 

SUBROUT INE  LPLACF(N0*N*X>Y»ALPHA»BETA1*BETA2*GAMMA*CH»SCH»IDIM»R* 
1 TScD»XmIN»XMAX.NX*YMIN*YMAX»NY*NaXDIM*NaYdIM) 

DIMENSION  X( IOIM.NO) -Y( id im*noj  * ALPHA  I IDIM*NO) *BETaK IDIM*NO) * 
1BETA2I I0IM»N0) »GAMMa( IDIM*NO) *CH( IDIM»NO) *TSCD< IDIM»N0>4) »N(NO) * 
2A( 252  *252  I *SCH(N0! * B ( 246  > 

COMMON/HELP/Z (252*252) 

EQU I VALENCE ( Z * A ) 

PI=3. 1415926 

RR=R*R 

DO  1 L=1 *N0 

NN=N ( L ) -1 

DO  1 I = 1 sNN 

X I = X ( 1 + 1,L)-X( I ,L) 

YI=Y( I+1»L)-Y( I»L) 

T SCD ( I*L*1)=ATAN2<Yi,XI) 

TSCD( I >L  » 2 ) =S IN ( T SCD( I » L » 1 > > 

TSCDl I »L  > 3 ) =COS ( T SCD l I *L»i>  > 

1 TSCD ( I *L  *4 ) =SORT (XI*XI+YI*YI ) 

J J J=0 

DO  4 LJ=l*NO 
N J=N ( LJ ) -1 
JAJ= J J J 
J J J= JJJ+NJ 
DO  4 J=1 *NJ 
JJ=JAJi-J 

XJ=lX(J*LJ)+X< J+1>LJ) )/2. 

YJ=( Y( J»UJ)+Y( J+l »LJ) ) /2. 

I 11=0 

DO  4 LI=l*NO 
N I =N ( LI ) -1 
I 1 1 = ! II+NI 
I A I = I 1I-NI 
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DO  A I = 1 * N I 
1 1 = I A 1 + 1 

I F ( II.EQ.JJ)  GO  TO  3 
X1=XJ-X( I »l I ) 

X2=XJ-X( I+1.LI ) 

Y1=YJ-Y( I »L I ) 

Y2=YJ-Y ( 1+1 *Ll ) 

R1=X1*X1+Y1*Y1 
R2=X2*X2+Y2*Y2 
S1  = 0. 

S2=0. 

YT  = Y( I * L I ) #X2-Y ( I +1 .LI )*X1+YJ*(X( 1 + 1 »LI ) -X ( I »LI  ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT.XT) 

IFIALPKA(J.LJ) .EQ.U. ) GO  TO  2 

Sl=TSCD« I *LI »A)*( l.-0.5*ALOG<R2/RR) ) +0. 5*AL0G ( R2/R1 ) * ( Xl*TSCD ( I »LI 
] »3>+Yl*TSCD< I »LI *2) )+TfTa*<Xi*TScd( I »LI »2 )-Yl*TScD( I *LI »3 ) > 

? TETA1=TSCD( J»LJ>l)-TSCDt I * L I » 1 ) 

S2=0.5*GlN(TETAl)*AL0G(R2/Rl )+COS(TETAl )*TETA 

S3=-S2 

GO  TO  A 

3 Si=TSCD< I »L  I »4)*( l.-AL06<TSCD( I ’Ll » A ) /2 • /R ) ) 

S2=-P I 

S3=-PI 

4 A ( J J . 1 1 ) =ALPHA  ( J »L  J ) *Sl+BETAl  < J »L  J ) *S2+8ET A2  < J * L J ) *S3 
M=0 

DO  5 L= 1 .NO 

5 M=M+N ( L ) -1 
J J J=0 

DO  6 L= 1 .NO 
NN=N ( L ) -1 
JAJ= JJJ 
JJ J= JJJ+NN 
DO  6 J=1.NN 
JJ= JAJ+J 

6 B( JJ)=GAMMA( J.L) 

CALL  ARRAY ( 2 .M.M.MAXD I M.NAYD I M . A . A ) 

CALL  SIMQlA.B’M.KS) 

IF  (KS.NE.O)  PRINT  lOu 
100  FORMAT ( lHO.IBHSYSTFM  IS  SINGULAR) 

JJJ  = 0 

DO  7 L = 1 .NO 
NN=N(D-1 
JAJ= JJJ 

jjj=jjj+nn 

DO  7 J= 1 >NN 
JJ= JAJ+J 

7 CH( J.L)=B( JJ) 

DO  8 L=1 .NO 
NN=N ( L ) -1 
SCH(L)=u. 

DO  8 1= 1 >NN 

n SCH(L)=SCH(L)+TSCD( I »L»A)*CH( I *L ) 

I F < NX-1 117.9.10 
9 DX=0. 

GO  TO  11 

10  DX= ( XMAX-XMIN) /FLOAT (NX-1 ) 

11  I F ( NY-1 ) 17.12*13 
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12  DY =0 • 

GO  TO  14 

13  DY=(YMAX-YMIN)/FLOAT<NY-D 

14  DO  16  11=1 .NX 
XJ=XMIM+FLOAT( I I-1)*DX 
DO  16  J J=1 *NY 

yj=ymin+float<jj-d*dy 

A ( I 1 *JJ)=0. 

DO  16  L I =1 .NO 
NN=N(Ll)-l 
DO  16  1=1 .NN 
X1  = XJ-X( I *L I ) 

X2=XJ-X ( 1+1 »LI ) 

Yl  = YJ-Y ( I »L  I ) 

Y2=YJ-Y( 1+1 »L I ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

IF((R1.F,Q.O.).OR.(R2.EQ.O.))  GO  TO  15 

YT  = Y( I »LI)*X2-Y( I+l,LI)*Xl+Yj*(X( I + 1 *L1 ) -X ( I * LI ) ) 

XT=Xl*X2+Yl*Y2 
TETA=ATAN2(YT ,XT) 

Sl  = TSCD< I »LI »4)*( l.-0.5*ALOG(R2/RR) >+0.5*aLOG(R2/R1 )*(-  1*TSCD( I + L I 
1 »3>+Yl*TSCD< I *L1 *2) )+TETa*<Xi*TSCD« I *LI .2)-Yl*TScDl I*LI »3> > 

GO  TO  16 

15  Sl=TSCD< 1 *L 1 .4)*! l.-0.5*ALOG! (R1+R2J/RR) ) 

16  A ( 1 1 »JJ) =A( 1 1 .JJ)+Sl*CH< I »L I ) 

17  RETURN 
END 

SUBROUTINE  CPWE < W »S .H» L .CHI > 

DIMFNSION  N(81  »X(7u,3)  .Y(70»3)  » ALPHA!  70  >3)  *BETA1(70>3)  » 

1 BF^ A2  ! 7^  . 3 ) »GAMMA<  70  »3  ) »CH<  70  >3)  »Sch!  3)  . T ScD ( 70  »3  »/>  1 .* A1  ( 181  > 181) 
REAL  L 

COMMON/ HE LP/Z (252.252) 

EOU I VALENCE l Z » A 1 > 

N0=3 
N ( 1 ) =41 
N ( 2 ) =70 
N ( 3 ) =70 

E0  = 1 •/( 4 »*3 • 14159*2 •99776*2 •9V776E 
FACT0R= ( 1 . oE-05 ) * ! 1 . /O „ 393700 ) 

DELW=W/40. 

DELG=(L/2.-W/2.-S)/69. 

X(l,l)=-W/2. 

Y( 1 . 1 ) =u, 

DO  26  1=1.40 
IP1=I+1 

X!  I PI . 1 ) =X  ( I .D+DELW 
Y ( I PI » 1 ) = Y ( I .1) 

ALPHA! I . 1 ) = 1 . 

BETAl! I » 1 ) =0. 

BETA2  < I ♦ 1 ) =0. 

26  GAMMA! I »1)=1. 

X ( 1 , 2 ) =-L/2  • 

Y( 1 *2 ) =0, 

X ( 1 » 3 ) =L/2 • 

Y ( 1 >3 ) =0 , 

DO  27  1=1.69 
I P 1 = 1+1 

X! I PI » 2 ) =X ( I .2 )+DELG 
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v 


f Y ( I PI »2 ) =Y( 1 .2 ) 

* X( I PI *3 ) =X ( I .3J-DELG 

j Y( I PI .3 ) =Y ( 1.3) 

! ALPHA! I »2)=1. 

\ BETA1(I*2)=0. 

I BETA2I I »2>=0. 

1 GAMMA ( I »2 1=0. 

t ALPHA! 1 .31=1. 

! BETAl t I .3 ) =0. 

BETA2! I »3)=0. 

27  GAMMA! I .3 1=0. 
j XM1N=0. 

| XMAX=0. 

I YMIN=0. 

YMAX=0. 

NX-0 

NY=0 

j ID IM=70 

I R=1 .0E+05 

I NAXDIM-181 

! NAYDIM=181 

CALL  LPLACEIN0»N*X*Y»ALPHA»BETA1>BETA2*GAMA»CH*SCH*ID1M»R>TSCD» 
1XMIN*XMAX»NX»YMIN*YMAX*NY»NAXDIM»NAYDIM) 

CH1=SCH ( 1)#6.2831852*E0 
RETURN 

, END 

SUBROUTINE  LPLACE<N0.N»X»Y»ALPHA»BETA1»BETA2»GAMMA»CH»SCH»IDIM»R» 

l tscd»xmin»xmax,nx,ymin»ymax»ny,naxdim»naydim> 

DIMENSION  X! IDIM»NO) »Y( IDIM.NO) .ALPHA! IDIM»NO) .BETAl ! IDIM.NO) » 
1BFTA2 ! IDIM.NO) .GAMMA* IDIM.NO) »CH« IDIM »N0 ) . TSCD < IDIM.N0.4) »N«NO) > 
2A1 ( 181.181 ) »SCHt NO) *8! 180) 

COMMON/HELP/Z! 252.252) 

EQU1VALENCE!Z.A1) 

PI=3.141592& 

RR=R*R 
DO  1 L= 1 .NO 
NN=N ( L ) -1 
DO  1 1=1, NN 
XI=X<  I + l.D-X!  I.L) 

YI=Y(  I + l.D-Yt  I.L) 

TSCD(I»L»1)=ATAN2(Yi,XI ) 

TSCD! I»L,2)=SIN(TSCD< I »L»1>  > 

TSCD ! I »L,3) =COS!TSCD< I »L»1 ) ) 

1 TSCDt I »L  .4 ) =SORT (XI#Xl+Yl*YI) 

JJJ  = 0 

DO  A LJ=1,N0 
NJ=N ( LJ ) -1 
JAJ=JJJ 
J J J= JJJ+NJ 
DO  4 J= 1 ,NJ 
J J= JAJ+J 

XJ=!X«J»LJ)+X! J+l.LJ) >/2. 

YJ=( Y( J»LJ)+Y( J+l.LJ) )/2. 

1 11=0 

DO  A LI =1 .NO 
N I =N ( L I )-l 
I I I = 1 1 1 +N I 
I A I = 1 1 1 -N I 
DO  A 1 = 1, NI 
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I 


L ■ 

W i 


i 


I I-IAI+I 

IF ( 1 1 .EQ.JJ)  GO  TO  3 
X1=XJ-X( I *LI > 

X2=XJ-X< 1 + 1 *LI ) 

Y1=YJ-Y(I*LI) 

Y2=YJ-Y< 1+1 *LI ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

S1=0. 

S2=0. 

YT=Y ( I .LI )*X2-Y( I+1,U)*X1+YJ*(X( 1+1 »L I ) -X ( I »L I ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT,XT) 

IF(ALPHA(J»LJ).EQ.u.)  GO  TO  2 

Sl  = TSC0l  I *LI  *4)*(1.-o.5*ALOG<R2/RR)  )+0.5*ALOG<R2/Rl  >*(Xl*TSCDl  I .LI 
1 »3>+Yl*TSC0(  I »LI  * 2 > )+TFTA*<Xl*TSCD(  I *LI  * 2 > -Yl*T$CD ( T *.L I * 3 ) ) 

2 TeTa1=TScDU.I.J»1)-TSCD(  I *L 1 . 1 ) 
S2=0.5*SlN(Te7Al)*ALOG<R2/Rl>+COS(TETAl)*TETA 
S3=-S2 

GO  TO  A 

3 S'1  = TSCD<  I » L I »A)*<  l.-ALOGi  TSCDl  I »L I »A ) /2 . /R ) ) 

■ S2=-P I 

S3=-P I 

A A( JJ. I I )=ALPHA< J»LJ)*Sl+BETAl ( J.LJ)*S2+BETA2( J»LJ)*S3 
M=0 

DO  5 1=1. NO 

5 M=M+N(U-1 
JJJ=0 

DO  6 L=l.NO 
NN=N ( L ) -1 
JAJ=JJJ 

jjj=jjj+nn 
DO  6 J=  1 »N.N 
J J= JAJ+ J 

6 B( JJ)=GAMMA( J.L) 

CALL  ARRAY<2.M>M»NAXDIM»NAYDIM»A»A) 

CALL  SIMO(A.B.M.KS) 

IF  (KS.NE.O)  PRINT  100 
100  FORMAT ( 1H0  * 18HSYSTEM  IS  SINGULAR) 

J JJ=0 

DO  7 L= 1 .NO 
NN=N ( L ) -1 
J AJ= J J J 
J J J= JJJ+NN 
DO  7 J=1.NN 
J J=JAJ+ J 

7 CHI J.L) =B< JJ) 

DO  8 L= 1 .NO 
NN=N I L ) -1 
SCH ( L ) =0, 

DO  8 1=1. NN 

8 SCH I L ) =SCH I L ) +TSCDI I .L»A)*CH(  I .L) 

I F I NX~1 ) 17 .9 » 10 

9 DX=0. 

GO  TO  11 

10  DX= I XMAX-XMIN ) /FLOAT  I NX-1 ) 

11  IFINY-I 117.12.13 

12  DY=0. 
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GO  TO  14 

13  DY=(YMAX-YMIN) /FLOAT ( NY-1 J 

14  DO  16  I I =1 »NX 
XJ=XM1N+FL0ATU  I-1)*DX 
DO  16  J J=1 *NY 
Yj=YMIN+FLOAT( JJ-1 )*DY 
A ( 1 I »JJ ) =0« 

DO  16  L I =1 »NO 
NN=N ( L I ) -1 
DO  16  1 = 1 *NN 
X1=XJ-X< I »L I ) 

X2  = XJ-X  ( 1 + 1 tl_  I ) 

Y1=YJ-Y( I »L I ) 

Y2=YJ-Y( I+1.LI ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

I F ( (Rl.EQ.O.) .OR. (R2.EQ.0.) ) GO  TO  15 

YT  = Y<  I»in*X?-Y<  1+1  »L  I ) #Xl+Yj*  ( XI  1+1  »L  I )~X  ( I »L  I ) ) 

XT=X1*X2+Y1*Y2 
TETA=ATAN2( YT.XT) 

SlsTSCD< I *UI *4)*< l.-0.5*AL0G(R2/RR) ) +0.5*AL0G l R2/R1 )* ( Xl*TSCD ( I ♦ L I 
1 ’ 3 ) +Yi#TSCD  ( I » l_ I * 2 ) ) +TETA* < Xl*TSCD  1 1 »LI  *2>-Yl*TSCD<  I *LI  » 3 > > 

GO  TO  16 

15  Sl=TSCD< I »LI »4)*< l.-0.5*ALOG( (R1+R2I/RR) ! 

16  A<  1 1 *JJ)=A( 1 1 »JJ)+Sl*CH< I »LI) 

17  RETURN 
END 

SUBROUTINE  ARRAY  < MODE » I » J »N »M *S *D ) 

DIMENSION  SHI »D( 1) 

NI=N-I 

I F (MODE-1 1 10U»1U0*12U 
lUo  IJ=I*J+1 
NM=N*J+1 
DO  110  K=1 i J 
NM=NM-N I 
DO  110  L=1 s I 
I J= I J-l 
NM=NM-1 

110  0(NMI=S( IJ) 

GO  TO  140 
12o  I J=0 

NM=0 

DO  130  K=1»J 
DO  125  L=1 » I 
I J=I J+l 
NM=NM+1 

125  S ( I J ) =D ( NM) 

■13o  NM=NM+N I 
140  RETURN 
END 

SUBROUTINE  SIMQ<A»B*N»K;S) 

DIMENSION  All) »8( 1) 

TOL=0.0 
KS=0 
J J=-N 

DO  65  J= 1 »N 
JY=J+1 
J J=J J+N+l 
BIGA=0« 
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. t 


t 


* . 


i- 


1 


E 


I T=JJ"J 
DO  30  I=J»N 
I J= I T+I 

IF(ABS(BIGA)-ABS( A( IJ) ) ) 20,30,30 
20  8IGA=A< I J) 

IMAX= 1 
30  CONTINUE 

I F ( ABS ( B IGA )-TOL ) 35,35  ,40 
35  KS=1 
RETURN 

40  1 1= J+N* I J-2 ) 

I T= IMAX-J 
DO  50  K=J»N 
11=1 1+N 
1 2= I 1+1 T 
SAVE=A(  ID 
A(  II) —A (12) 

A ( 1 2 ) =SAVE 
50  A( 1 1 ) =A ( 1 1 ) /B IGA 
SAVE=B( IMAX) 

B ( 1MAX)=B(J) 

B( J)=SAVE/BIGA 
IF(J-N)  55,70,55 
55  IQS=N*(J-D 
DO  65  I X=JY  ,N 
IXJ=IQS+IX 
I T=J- IX 
DO  60  JX=JV»N 
I XJX=N#  < JX-1 ) + IX 
JJX= IXJX+IT 

60  A(IXJX)=A< IXJXI-I A< IXJ)»A(JJX) ) 

65  B( IX)=B( 1 X ) - ( B ( J)*A(  IXJ) ) 

70  NY=N-1 
I T =N*N 

DO  80  J=1»NY 

ia=it-j 

I B=N-J 
I C=N 

DO  80  K=1  * J 

B ( I B ) =B ( I B ) -A ( I A ) #B ( IC) 

IA=IA-N 
80  IC= IC-1 
RETURN 
END 
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PROGRAM  IM( INPUT »OUTPUT ) 

REAL  L 

READ  1*  NSETS 

1 FORMAT (1 10) 

DO  2 MNel *MSETS 

READ  3*  W>H»  T »L>ER  *TAND1 »S IGMA 
3 FORMAT ( 7F10»6) 
tandi=tandi*i.oe-04 
6 16ma=SIGMA*1»OL+07 
PRINT  8 

b format«ih  *////ih  >*strip  metalization  thickness  = 250  microinches 
1*) 

PRINT  4»  W*H*T»L»ER»TAND1»SIGMA 
A FORMAT! 1hO*2HW=*F5.1»3X,2HH=*F5.1*3X»2HT=>F5.1*3X*2HL=> 
1F10.6»3X*3HER=»F1u.6*3X,6HTAND1=»E10.3»3X»6HSIGMA=*E10.3) 

FACTOR=( l.oE-05)*! l./u. 393700) 

W=W*FACT0R 

H=H*FACT0R 

T=T*FACTOR 

L=L*FACTOR*(1.E+03) 

E1=ER 

E2=ER+0.01*ER 

CALL  SUSmSF(W.H»T»L*E1 »S IGMA »C »PDDSQ ) 

pddsqf=pddsq 

CALL  SUSmSF(W>H>T.L»E2*SIGMA»CP»PDDS0) 

CALL  SUSmSE!W»H»T.L*CE> 

A=C*CE 

B=CE/C 

ZD=l./( (3.E+^8)*SQRT(A)  ) 

V=3.E+08*SQRT(B) 

ARG1=B 

ARG2=CE/CP 

ER1  = 1./'ARG1 

ER2=1./ARG2 

FACTC  = l>J»/2«3 

FACTD=27.3/3.E+o8 

AC  = FACTC*PDUSQF#i:Rl*Zv 

AD=FACTd*(F1/SQRT(ER1 ) ) * ! (EK2-ER1)/<C2-E1) )«TaND1 
PRINT  5 

5 FORMAT ( 1HU  » 1X»8HZ0 ( OHMS) » 5X > 8HV l M/SEC ) »5X» 

120HALPHAC<DB/M)/SQRT(F) » 5X  * 14HALPHAD ( DfcJ/M ) /F ) 

PRINT  6 » ZO  » V » AC  > AD 

6 FORMAT ( lH  *2X*F6.2»4X»E10.3»8X»E10.3*11X.E10.3) 

2 CONTINUE 
STOP 
END 

SUBROUT  I NE  SUSMSF ( W ,H ♦ T , L » ER  * S I GMA  *CHRG  jPDDSQF ) 

DIMENSION  N(6)*X(51*6)»Y(51>6) * ALPHA! 51 »6) >BETA1<51»6) » 

1BET  A2  < 5 1 »6 ) * GAMMA! 51 >6) *CH<51 >6) >SCH(6) »TSCD( 51 >6>4)«A<238»238) 
REAL  L 

COMMON/HELP/Z ( 238 » 238 ) 

EQUIVALENCE ( Z »A ) 


NO 

i=6 

N ( 

1) 

=43 

N ( 

2) 

=41 

N ( 

3) 

=41 

N ( 

A! 

= 31 

N ( 

5) 

= 31 

N ( 

6) 

= 51 
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•E0=1./ (4. *3. 14 J 59*2.99776*2. 99776E+09) 
AAA=3. 14 159*4 .E-07*3» 14159/SIGMA 
RDSQF=SQRT( AAA) 

FACTN=4. *3. 14159*3. 14159*RDSQK*8.85E-12 
FACTD=4.*3.14159*1.E-U7 

facti.=factn/factd 

FACT0R=(l.0E-05)*( l./O. 393700) 
DELW=W/20. 

E=ER*E0 

X(l»l)=-W/2. 

Y ( 1 * 1 ) =-0.25*FACTOR 
DO  26  1=1*20 
IP1=I+1 

XI  IP).»1)=X<  I *1)+DELW 

26  Y( I PI  * 1 ) =Y ( 1 * 1 ) 

X(22*1)=X<21*1) 

Y C 22  * 1 ) =0. 

DO  27  1=22*41 
IP1=I+1 

X< IP1,1)=X( I *1 J-DELW 

27  Y ( IP1 >1 )=Y( I *1) 

X ( 43  » 1 ) =X 1 1 * 1 ) 

Y ( 43  * 1 ) =Y ( 1 * 1 ) 

DELG=(L/2.)/4U. 

X(l,2)=-L/2. 

Y ( 1 *2 ) =“S 
X( l»3)=L/2. 

Y ( 1 * 3 ) =~S 
DO  28  1=1*40 
IP1=I+1 


X(  I PI *2 ) =X ( I *2 ) +DELG 

Y ( I PI *2 ) =-S 

X( I PI  * 3 ) =X ( I >3)-DELG 

28  Yt I PI *3 ) =-S 

DELSUB= (L/2.-W/2. ) / 30 • 

X ( 1 *4 ) =-L/ 2 » 

Y( 1 ,4)=u. 

X( 1 *b)=L/2. 

Y ( 1 *5 ) =0. 

DO  29  1=1*30 
IP1=I+1 

X< I PI *4 ) =X ( I *4 ) +DELSUB 
Y ( I PI *4 ) =0. 

XI IP1*5)=X(I . 5 ) -DELSUB 

29  Y ( I PI *5)=0. 

DELH=H/2 • 

DELL=L/46* 

X ( 1 *6 ) =-L/2 . 

Y( 1 ,6 ) =0. 

X ( 49  *6 ) =L/2 . 

Y ( 49  *6 ) =H 
DO  30  1=1*2 
IP1=I+1 
I P49= 1+49 
IP49M1=IP49-1 
X( IP1 *6 ) =X ( I ,6) 

Y ( IP1 *6)=Y( I *6 ) +DELH 
X( IP49»6)=X( IP49M1.6) 

30  Y( IP49*6)=Y( I P49M1 *6 > -DELH 


,:-4r 


A»f  1 


ti 


.V 


MS 
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DO  31  1=1*46 
IP3= 1+3 
1P3M1=1P3-1 

X( IP3»6)=X( IP3M1.61+DLLL 
31  Y ( I P3  *6 ) =H 
DO  60  1=1 *42 
ALPHA( 1 • 1 ) =1 • 

BETAK  1 *1)=0. 

BETA2U  *1)=0. 

60  GAMMA ( I * 1 1 =1 • 

DO  61  I =1 *40 
ALPHA ( I * 2 1 “ 1 • 

ALPHA ( 1 *3 ) -1 • 

BETAK  I*2>=0. 

BETAK  I *31=0. 

BETA2 ( I *21=0. 

BETA21 1*31=0. 

GAMMA< I *2 1 =0. 

61  GAMMA ( 1 *3 1=0. 

DO  62  1=1*30 
ALPHA ( 1*41=0. 

ALPHA( 1*51=0. 

BETAK  I»4l=E0 
BETA1 ( I » 5 1 =E 
BETA2 ( I *4 1 =E 
BETA2« I * 5 1 -EO 
GAMMA! I *4 1 =0, 

62  GAMMA ( 1*51=0.. 

DO  63  1=1*50 
ALPHA ( I *6 1 =0. 

BETAK  1,6>=E 
BETA2I I .6 1 =E0 

63  GAMMA ( I *6 1=0. 

XM1N*0, 

XMAX=0*. 

YMIN=0. 

YMAX=0. 

NX=0 

NY=0 

ID IM=51 

R=1.0E+05 

NAXDIM=238 

NAYDIM=238 

CALL  LPLACF(M0*N*X>Y»ALPHA»BETA1*BETA2*GAMMA*CH»SCH*IDIM»K*TSCD. 
1XMIN»XMAX»NX*YMIN.YMAX*NY*NAXDIM*NAYDIM> 

DELF=W/20» 

DELE=0.25*FACTOR 

FF1=6«2831852*E#DELF 

FF2=6.2831852*E0*DELE 

FF4=6.2831852*EO*DElF 

CHRG=0. 

DO  50  1=1.20 

50  CHRG=CHRG+FF4*CH(  I .1) 

DO  51  1=22.41 

51  CHRG=CHRG+FFl*CH( 1*1) 

CHRG=CHRG+FF2*  < CH ( 21  * 1 1 +CH < 42  * 1 1 1 
SUM=0. 

DO  64  Jsl.NO 

IF(ALPHA(l»J).NE.l.)  GO  TO  64 
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NJ1=N<J)-1 
DO  65  1=1 >NJ1 
IP1=I+1 

SEGX=X( IP1.J)-X( l»J) 

SEGY=Y(IP1.J)-Y(I,J) 

segarg=segx*segx+segy*segy 
3cGLEN=SQRT  ( SEGARG) 

65  SUM=SUM+CH< I .J)*CH<  I »J»*SEGLEN 
64  CONTINUE 

PDDSQF=FACTL*SUM 
54  CONTINUE 
RETURN 
END 

subroutine  lplacf <no.n»x»y.alpha»betai*beta2 .gamma .ch.sch. i dim. r> 

1 T SCO » XM I N » XMAX  >NX»YMIN» YMAX  *N Y *NAXD i M *NAYD I M ) 

DIMENSION  X( IDIM*NO) »Y( IDIM.NO) »ALPHA( IDIM»NO) .betak IDIM.NO) . 
1BETA2UDIM»N0)  » GAMMAt  IDIM»NO)  .CHI  IDIM.NO)  »TSCDt  IDIM»NO*4)  >N(NO)  » 
2A( 238*238) .SCH(NO) * B t 237 ) 

COMMON/HELP/Z 1 238  > 238 ) 

EQUIVALENCE^,  A) 

PJ=3.1A1592S 

RR=R>R 


DO  1 L= 1 .NO 
NN=N ( L ) -1 
DO  1 1=1. NN 
XI=X( I+1,L)-X( I .L) 

Y I = Y ( 1 + 1 » L ) -Y ( I ,U 
TSCDU»L.1)=ATAN2(YI,XI  ) 

TSCDI I »L  »2 ) =SIN ( T SCO ( I *L  * 1 ) ) 

TSCD ( I >L .3) =COS ( TSCDI I »L » 1 ) ) 

1 TSCDI  I»L»4)=SQRT(XI*XI+YI*YI) 

JJ  j=o 

DO  4 LJ=1 »NO 
NJ=N  I LJ  ) — 1* 

JAJ=JJJ 
JJJ=JJJ+NJ 
DO  4 J=  1 »NJ 
J J= JAJ+ J 

X J = ( X I J f L J ) +X  ( J+l  . i_  j ) ) /2  , 

YJ=(Y(J»LJ)+Y(J+1,LJ) ) / 2 • 

II  1=0 

DO  4 1.1  = 1 .NO 
N I =N ( L I ) -1 
111=111 +N I 
I A I = 1 1 1 >-N  I 
DO  4 1=1. NI 
I I = I A 1+ 1 

I F I 1 1 .EQ.JJ)  GO  TO  3 
X1=XJ-XU  »L I ) 

X2=XJ-X I I+l.LI ) 

Y1=YJ-Y(I,LI) 

Y2=YJ-Y I 1 + 1 .LI  ) 

Rl=Xl*Xl+Yl*Yl 
R2=*2*X2+Y2*Y2 
S1=0. 

S2=0. 

YT  = Y( I »LI )*X2-Y< I+1,LI)*X1+YJ*(X( I + l.LI )-X(  I .LI ) ) 

XT=Xi*X2+Yl*Y2 

TETA=ATAN2(YT,XT) 
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IF ( ALPHA (J»LJ)  . EQ«>J.  ) GO  TO  2 

Si  = TSC0l  I *LI  * A ) * C l.-0.5*ALOG(R2/RK)  )+U.5#ALOG*R2/Kl  J t X l «-T  i>CO  * I *LI 
l,3)+Yl*TSCO( I *LI »2> )+TETA*<Xi*TSCD( I»LI >2)-Yl*TSCD( I »LI »3) ) 

2 TETA1=TSCD l J »LJ* 1 )~TSCD( I »L I * 1 ) 
S2=0,5*SIN(TETA1)*ALOGIR2/R1)+COSITETA1)*TETA 
S3=-S2 

GO  TO  A 

3 Sl  = TSCD<  I »LI »4>*(  l.-AL0G(T6CDl  I »L I »4 ) /2 . /R > > 

S2=-P 1 

S3=-P I 

4 A(  JJ»II  J=ALPMA<J*LJ)*S1-*-BETA1(J»LJI»S2+BETA2«J»LJ)*.S3 
M=0 

00  5 L=1 *N0 

5 M=M+N(L ) -1 
JJJ=0 

DO  6 L= 1 »N0 
NN=N(L>-1 
JAJ= JJJ 
JJJ= JJJ+NN 
DO  6 J=1*NN 
JJ=JAJ+J 

6 B( jj)=GAMMA(J»L) 

CALL  ARRAY(2»M*M»NAXDIH.NAYDIM»A»A) 

CALl  SIMQ(A»B*M*KS) 

IF  (KS.NE.u)  PRINT  1<jo 
100  FORMAT! 1H0.1BHSYSTEM  IS  SINGULAR) 

JJJ=0 

DO  7 L= 1 *N0 
NN=N(L)-1 
JAJ= JJJ 
JJJ=JJJ+NN 
DO  7 J=1»NN 
„ J J=JAJ ►J 

7 CH( J iL ) =B ( J J ) 

DO  8 L=1»N0 

NN=N(L>-1 

SCHIL1-U. 

DO  8 1=1 »NN 

8 SCH«D=SCH(L)+TSCD(  I *L*4)#CH<  I »L > 

IFlNX-1 ) 17  * 9 * lo 

9 DX=0. 

GO  TO  11 

10  r>X=(XMAX-XMlN)/FLOAT(NX-l ) 

11  I F ( NY-1 )17»12*13 


12  DY=0. 

GO  TO  14 

13  DY=(YMAX-YMIN)/FL0AT(NY-1) 

14  DO  16  I I = 1 *NX 

X J=XM I N+FLOAT (II - 1 ) *DX 
DO  16  JJ=1 »NY 
YJ=YMIN+FL0AT(JJ-1)*DY 
A ( I I »JJ)=U. 

DO  16  L 1=1 »N0 
NN=N(LI )-l 
DO  16  1=1 »NN 
X1  = XJ-X(  I * L I ) 

X2=XJ-X( 1+1 »LI ) 

Y1  = YJ-Y( I >L I ) 
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Y2=YJ-Y( 1+1.L1 > 

Rl=Xl»Xl+Yl*Yl 

R2=X2*X2+Y2*Y2 

IF((Rl.EQ.U»).OR.(R2.EO.O.n  GO  TU  15 
YT=YU.Ln*X2-Y(I+l,LI)*Xl+Yj*<X(Hl»LI  l-X(I*LID 
XT=X1*X2+Y1*Y2 
TETA=ATAN2(YT»XT) 

Sl  = TSCD(  I»LI  »A>*(  1.-o(,5*aLOK(R2/RR)  ) +u. 5*ALOG l R2/R1  \*  !Xl*TSCD< 1 >LI 
l»3>+Yl*TScr>(  I *LI  *2)  J+TETa# l Xl*TSCDl  I >LI  >2l-Yl*TbCDt  1 »Li  *3)  ) 

GO  TO  16 

15  Sl=TSCD< !*Ll»4>*< l.-U.5*ALOG( IR1+R2I/RR) ) 

16  A(  II »JJ)=A( II >JJ)+Sl*CH! 1 >L  I ) 

17  RETURN 
END 

subroutine  susmse!w>h>t»l>chii 

DIMENSION  N«3) »Xt7o,3) »YC70»3) > ALPHA « 70 * 3 ) >8ETa1(7U>3) . 

1BE^A2 ( 70»3 ) »GAMMA(7w*3) >CHl70*3) .5CHI3! > T5CD 1 7U *3 » A ) *AH 181 » 181 1 
REAL  L 

COMMON/ HELP/2 ( 238*238) 

EQU IVALFNCE l Z ♦ A1 ) 

N0  = 3 
N(1)=M 
N ( 2 J =70 
N l 3 ) =70 

EO=l./!A.*3.1AI59*2.99776*H.9977t>fc.i-09) 

FACTOR^ tl.UE~O5)*(l./0.?937U0) 
delw=w/ao. 

OELG=(L/2.)/69. 

X ( 1 * 1 ) =-W/2 • 

Y ( 1 1 1 J =0 . 

DO  26  J-1»A0 
IP1=1+1 

XUP1*1)=X(I  .1I+DELW 
Y ( I PI » 1 ) =Y ( I » 1 } 

ALPHA! I>1)=1.  ' 

9ETA1U  »1)=0. 

BETA2' I » 1 ) =0. 

26  GAMMA  1 1 s 1 ) = 1 . 

X( ? •? ) =-L /2  • 

Y ( 1 * 2 ) =-S 

X ( 1 >3)=L/2. 

Y ( 1 ■?•)  ~~5 
DO  27  I =1 >69 
IP1=1+1 

XI IPJ  >2 J=X( I >2 I+DELG 
Y ( I PI >2 ) =Y ( I >?) 

XI  I PI >3) =X( I >3 ) -DELG 
Y ( 1 PI > 3 ) =Y ( I >3) 

ALPHA ! I >2 1 = 1 • 

BETAK  I >2)=0. 

BETA2S I,2)=0. 

GAMMA! I >2>=0. 

ALPHA! I * 3 ) = 1 - 
BETA1 ! I >3  > =0. 

BETA2< I>3)=0. 

27  GAMMA! I >3>=0. 

XMIN=0. 

XMAX=0. 

ymin=o. 
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YMAX=0. 

NX=0 

NY=0 

IDIM=70 

R=1 »0E+U5 

NAXDiM=181 

NAYD1M=181 

CALL  LPL ACE !NOtNtXtY  t ALPHA iBETAI t BE TA2  >GAMMA»CHtSCH » IDIMtR  t TSCD t 
lXM!NtXMAX*NX»YMIN»YMAX»NY*NAXDIM»NAYDIM) 

CHl  = SCHll)*6.2831ti52*EU 

RETURN 

END 

SUBROUTINE  LPLACE< NO tNtXtY /ALPHA iBETAI *BETA2 i GAMMA tCHtSCHt IDIN*R * 
1 TSCD»XMIN»XMAX,NX,  t vlINtYMAX  *NY»NAXDIM»NAYDIM) 

DIMENSION  XI IDIMtNO) >Y( IDIMtNO) t ALPHA! IDlM*NO>  *BETa1( IDIM*NO) » 
1BETA2! IDIMtNO) ? GAMMA! IDIMtNO) tCH! IDIMtNO) tTSCDl IDlMtNOt A) tN(NO) t 
2Al(181tl81)  tSCH(NO)  *Bd8y) 

COMMON/HELP/Z ! 238 1 238 ) 

EQUIVALENCE ! Z t A1 ) 

P 1*3.1415926 

RR=R*R 

DO  1 L=ltNO 

NN=N!L>-1 

DO  1 I = 1 tNN 

XI=X!  !UvL)-X!  ItL) 

Y I =Y! 1+1 »L ) -Y ! I tL) 

TSCD! I tLtl)=ATAN2(YI tXl  ) 

TSCD!  I tLt2)=SlN!TSCD!  I » l_  1 1 1 ! 

TSCD! I tL  1 3 ) =COS( TSCD! 1 t L » 1 ) > 

1 TSCD(I»Lt4)=SQRT!Xl*XI+YI»Y!) 

JJJ  = 0 

DO  4 LJ  = 1 tNO 
NJ=N!LJ)-1 
JAJ=JJJ 
J J J= JJJ+NJ 
DO  4 J=ltNJ 
J J= JAJ+ J 

XJ=!X( JtLJ)+X< J+l tLJ) )/Z. 

Y J = ( Y ( J t|_J ) + Y ! J+l  1 1 I ) ) / 2 . 

I I 1=0 

DO  4 L 1=1 tNO 
N I =N ! LI ) -1 
I I I = I I I+NI 
I A I = 1 1 1 ~N ! 

DO  4 1=1 tNI 
I 1 = IAI+I 

IF! II.EQ.JJ)  GO  TO  3 
X1  = XJ-X( I * L I ) 

X2  = XJ-X ( 1 + 1 tLl  ) 

Y 1 =YJ-*Y ( I tLl ) 

Y2=YJ-Y! 1+1 tLl ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

Sl=0. 

S2  = 0. 

YT=Y ( I tLl )*X2-Y( I+ltLI >*Xl+YJ*!X< 1+ltLI >-X( I tLl ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2!YTtXT) 

I F < ALPHA <J tLJ) .EQ»0 • ) GO  To  2 
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Si=TSCD< I >L  I *4)*( l.-o.5*ALOG(R2/RR) ) +0, 5*ALOG ( R2/R1 )* ( Xl*TScD l I » L I 
] *3)+Yi*TScn< I *L I * 2 ) )+TETA* (Xl*TSCD« 1 »LI >2)-Yl*TScD< I »LI *3) ) 

? TETA1  = TSCD' J»LJ»1)-TSCD( I » L I * 1 ) 

S2=0.5*SIN(T£TA1)*ALOG<R2/R1)+COS(TETA1)*TeTA 
S3  = -S2 
GO  TO  4 

^ Sl=TSCD( I *LI .4)#! l.-ALOGl TSCD< I *L I *4 > /2. /R) > 

S2=-Pi 
S3=-P i 

A A(  JJ> 1 1 )=ALPHAI J»LJ)#Sl+BETAl ( J >L J ) *S2+BETA2 < J.LJ)*S3 
M=0 

DO  5 L=l»NO 
b M=M+N  ( I. ) -1 
JJJ  = 0 

DO  6 L--l.NO 
NN=N ( L ) — 1 
JAJ= JJJ 
JJJ= JJJ+NN 
DO  6 J=1,NN 
J J= JAJ+ J 

6 B( JJ)=GAMMA( J.L! 

CALL  ARRAY(2»',»M»NAXDIM»MAYDIM»A.A) 

CALL  S I MO (A.B.M.KS) 

IF  (KS.NE.u)  PRINT  lou 
100  FORMAT (1H0.18HSYSTEM  IS  SINGULAR) 

JJJ  = 0 

DO  7 L= 1 » NO 
NN=N ( L ) — 1 
JA J= J J J 
JJJ=JJJ+NN 
DO  7 J= 1 »NN 
JJ= JAJ+ J 

7 CH! J»L)=B< JJ) 

DO  8 L= 1 .NO 
NN=N ( L ) -1 
SCH(L)=U. 

DO  8 I = 1 ,NN 

p,  SCH(L)=SCH(L)+TSCD(  I »L*4>*CH<  I » L ) 

I F ( NX-1 ) 17  * 9 » lu 
9 DX-O. 

GO  TO  11 

10  DX= ( XMAX-XM IN) / F* OAT(NX-l) 

11  IF(NY-1)17.12.13 

12  DY=0 • 

GO  TO  14 

13  DY=(YMAX-YMIN)/FLOAT(NY-l! 

14  DO  16  11  = 1 .NX 
XJ  = XM1N+FL0AT< I I - 1 ) *DX 
DO  16  JJ=1,NY 
YJ=YMIN+FL0AT( JJ-1 )*DY 
A ( II  »JJ)=0. 

DO  16  L I = 1 .NO 
NN=N ( L 1 )-l 
DO  16  I = 1 >NN 
X1=XJ-X( I ,LI ) 

X2  = XJ-X( 1 + 1 .LI  ) 

Y1  = YJ-Y  t I ,LI  ) 

Y2=YJ-Y( I+1.LI  ) 
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R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

IF((R1.EO.O.)*OR.(R2*EQ*0.))  GO  TO  15 

YT=Y( I*LI)*X2-Y(I+l»LI)#Xl+YJ*(X( 1+1 »LI )-X( I *LI ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT.XT) 

Sl=TSCD«I*LI*A)*(l.-0.5*ALOG(R2/RR)  )+0«5*AL0G«R2/R1)*U1*TScD(  I*LI 

l,3)+Yl*TScr>U»LI»2)  )+TETA*‘Xl*TSCB(  I*LI  »2)-Yl*TScD<  I ’Ll  *3)  > 

GO  TO  16 

15  Sl=TSCDl I *LI .4)*< l.-0.5*ALOG( (R1+R2)/RR) * 

16  Ad  I : JJ  ) =A(  1 1 * JJ)  +Sl*CH( I *LI ) 

17  RETURN 
END 

SUBROUTINE  ARRAY  < MODE » I » J »N »M»S»D ) 

DIMENSION  S(l) »D< 1) 

NI=N-I 

IF (MODE-1)  100.100.1^0 
100  IJ=I*J+1 
NM=N*J+1 
DO  110  K=1 t J 
NM=NM-NI 
DO  110  L = 1 . I 
IJ=IJ-1 
NM=NM-1 

110  D( NM ) =S ( I J ) 

GO  TO  140 
120  IJ=0 
NM=0 

DO  130  K=1 . J 
DO  125  L=1 * I 
IJ=IJ+1 
NM=NM+1 

125  S( IJ)=D«NM) 

130  NM=NM+NI 
140  RETURN 
END 

SUBROUTINE  SIMQ(A.B»N»KS) 

DIMENSION  All) »B( 1) 

TOL=0.0 

KS=0 

JJ=-N 

DO  65  J=1*N 
JY= J+l 
J J= JJ+N+1 
BIGA=0. 

I T = JJ-J 
DO  30  I = J » N 
I J=IT+I 

I F ( ABS ( B IGA ) -ABS ( A ( IJ> ) ) 20.30,30 
20  BIGA=A( IJ) 

IMAX=I 
30  CONTINUE 

I F ( ABS ( BIGA l-TOL ) 35.35.40 
35  KS=1 
RETURN 

40  1 1 = J+N* ( J-2  ) 

IT= IMAX-J 
DO  50  K-J.N 
11=1 1+N 
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1 2 = 1 1+1 T 
SAVE=A< Il> 

At  1 1 > =A ( 12) 

At  1 2 ) =SAVE 
50  At  1 1 » = A t 1 1 ) /RIGA 
SAVE=B( I MAX) 

B(  IMAX)=B< J) 

Bt J)=SAVE/BIGA 
IFtJ-N)  55  »70*55 
55  IQS=N*( J-l ) 

00  65  i.x=jy*n 
IXJ=IQS+IX 

1 T=J**IX 

DO  60  JX=JY»N 
IXJX=N*< JX-1)+IX 
JJX=IXJX+IT 

60  At  IXJX ) =>A(  I XJX ) - ( A t I XJ  > *A  ( JJX  ) ) 
65  Bt IX) =B t IX ) - ( B ( J ) *A ( I X J } ) 

70  NY=N-1 
I T=N*N 

DO  80  J= 1 tNY 
I A= I T-J 
I B=N-J 
I C=N 

DO  80  K = 1 1 J 

Bt IB ) =B  < IB ) -At  I A ) t I C ) 

I A= I A-N 
80  IC=IC-1 
• RETURN 
END 
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Appendix  D 

LISTING  OF  PROGRAM  TIM 


PROGRAM  TIM 
REAL  L 

READ  1»  NSETS 

1 FORMAT ( 1 10) 

DO  2 MN=1 *NSETS 

READ  3*  W*S»H*T*L*ER*TAND1»SIGMA 

3 FORMAT ( 8F10.6) 

YAND1=TAND1*1.0E-04 
SlGMA=SlGMA*1.0£+07 
PRINT  8 

8 FORMAT!  1H  *////lH  **STRIP  METALIZATION  THICKNESS  = 250  MICROINCHES 
1*) 

PRINT  4*  W*S*H*T»L*£R*TAND1»SIGMA 

4 FORMAT! 1hO»2HW=.F5.1.3X»2HS=*F5.1»3X.2HH=*F5.1*3X*2HT=»F5.1» 
13X»2HL=*F10.6*3X.3HER=»F10.6*3X»3HTaND1=»E10«3»3X>6HSIGMA=*E10,3) 

FACTOR^ (l.0E-O5)*(l./O. 393700) 

W=W*FACTOR 
S=S*FACT0R 
H=H*FACTOR 
T=T*FACTOR 
L=L*F ACTOR* ! 1 • E+03 ) 

E1  = ER 

E2=ER+0.01*ER 

CALL  TIMF(W»S»H*T»L»E1 *SIGMA*C*PDDSQ) 

PDDSQF=PDDSQ 

call  TlMF<W,S»H>T,L*E2*SIGMA*CP*PDDSQ) 

CALL  TIME<W.S.H*T.L»CE> 

A=C*CE 

B=CE/C 

Z0=l./< ( 3 • E+08 ) *SQRT ( A ) ) 

V=3 • E+08*SQRT ( B ) 

ARG1=B 

ARG2=CE/CP 

ER1=1./ARG1 

ER2= 1 • / ARG2 

FACTC=10./2.3 

FACTD=27. 3/3. E+08 

AC=FACTC*PDDSQF*ERl*Zu 

AD=FACTD*(E1/SQRT!ER1) )*! ! ER2~hRl ) / ! E 2 — E 1 ) )*TAND1 
PRINT  5 

5 FORMAT! 1H0.1X.8HZ0! OHMS) * 5X  *8HV ( M/SEC ) >5X» 

120HALPHAC<DB/M)/S0RT(F) »5X»14HALPHAD(DB/M)/F) 

PRINT  6 > ZO » V *AC  > AD 

6 FORMATdH  *2X»F6.2»4X»E10.3»-8X»E10»3,11X»E10.3) 

2 CONTINUE 
STOP 
END 

SUBROUTINE  TIMF(W»S»H»T»L»ER*SI GMA *CHRG  *PDDSQF ) 

DIMENSION  N«6) »X!43*6) *Yl43*6) *ALPHA<43*6) >BETA1 (43 *6) * 
1BE1A2<43»6)  » GAMMA (43*6)  *CH<43*6)  *SCH(6)  > TSCD ( 43  *6  >4  ) »A(  181*181)  * 
2A1 (181*181) 

REAL  L 

C0MM0N/HELP/A*A1 

EQUIVALENCE(A*A1) 


NO 

= 

6 

N ( 

1 

) 

=43 

N ( 

2 

) 

=41 

N ( 

3 

) 

= 41 

N ( 

4 

) 

= 11 

6 8 


m 
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N < 5 ) =11 
N<6)=34 

E0=1./ (4. *3. 14159*2. 99776*2. 99776E+09) 
AAA=3. 14159*4. E-07*3« 14159/SIGMA 
RDSQF=SQRT(AAA) 

FACTN=4.*3.14159*3.14159*RUSQF*B.85E-12 

FACTD=4.*3.14159*l.E-07 

FACTL=FACTN/FACTD 

FACTOR^  <1.0E-05)*U./U.  393700) 

DELW=W/20. 

E=ER*E0 
X( l.l)=-W/2. 

Y ( 1*1) =-0.25*FACT0R 
DO  26  1=1 »20 
IP1=I+1 

X( IP1 >1 ) =X ( I *1 J+DELW 

26  Y I IP1»1)=Y( 1.1) 

X(22.1>=X(21»1) 

Y ( 2 2 * 1 ) =0* 

DO  27  1=22.41 
IP1=I+1 

X( I PI *1 ) =X ( I » 1 J-DELW 

27  Y ( I PI  * 1 ) =Y ( 1,1) 

X < 43  » 1 ) =X ( 1 » 1 ) 

Y ( 43  » 1 ) =Y ( 1 » 1 ) 

DELG=(L/2.-W/2.-S)/20. 

DELT=T/10. 

DELCW2=«W/2.+S)/10. 

X(l,2)=-L/2. 

Y ( 1 »2 ) =u# 

X( 2 1 .2 ) =~(W/2. )-S 
Y { 2 1 » 2 ) =0. 

X(31,2)=X(21»2) 

Y ( 3 1 ,2 ) =-T 
X ( 4 1 .2 ) =0. 

Y ( 41 » 2 ) =-T 
X ( 1 ♦ 3 ) =L/2» 

Y(l,3)=0. 

X(21»3)=(W/2.)+S 
Y( 21 ,3 ) =0. 

X ( 3 1 .3 ) =X ( 2 1 .3 ) 

Y ( 3 1 » 3 ) =-T 
X ( 41 .3 ) =0. 

Y ( 41 .3 ) =-T 
DO  28  1=1,19 
IP1=!+1 

X' I PI ,2 ) =X ( 1 ,2 ) +DELG 
Y ( I PI ,2 ) =0. 

X< IP1,3)=X( I ,3 ) -DELG 

28  Y ( I PI »3 ) =0. 

DO  228  1=1,9 
IP2 1= 1+2 1 

I P2 1M1= I P21~l 
X( IP21,2)=X( IP21M1.2) 

Y ( IP21,2)-YUP21M1,2)-DELT 
X( IP21,3)=X( IP21M1.3) 

Y ( I P21 »3 ) =Y ( IP2 1 , 2 ) 

IP31=I+31 
IP31M1= I P31-1 
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X(  IP31*2)=X! IP31M1 * 2 1+DELCW2 
Y ( I P31 *2 ) =“T 

X(  I P31 *3 ) =X ( IP31M1 *3 )“DELCW2 
228  Y ( I P31 *3 ) =Y  t IP31»2) 
DELS=S/10. 

X ! 1 *4')  =~W/2  •~S 
Y(l*4>=0. 

X(l»5)=W/2.+S 

Yll*5)=0. 

DO  29  1=1*10 
IP1=I+1 

XI IPl*4)=X( I *4 ) +DELS 
Y ( IP1,4)=0. 

X! IP1*5)=X( I *5) -DELS 

29  Y( IP1*5)=0. 

DELH=H/2. 

DELL=L/29. 

X( l*6)=-L/2. 

Y ( 1 *6 ) =0 , 

X ( 32  *6 1 =L/2  • 

Y ( 32  *6 ) =H 
DO  30  1=1*2 
IP1=I+1 
I P32= 1+32 
I P32M1= I P32~l 
X( I PI -6 ) =X ( I *61 
Y( IP1 *6) =Y< I *6 ) +DELH 
X( IP32*6)=X( IP32M1*6) 

30  Y(  I P32 * 6 ) = Y < IP32M1  *6)**DELH 
DO  31  1=1*29 

IP3= 1+3 
1P3M1=IP3-1 

X ( I P3 » 6 ) =X ( I P3M1 *6 ) +DELL 

31  'i  (1 P3 *6 ) =H 
DO  60  1=1*42 
ALPHA! I *1>=1# 

betai ( I * 1 ) =o« 

BETA2 ( 1 > 1 ) =0. 

60  GAMMA! I *1>=1. 

DO  61  1=1*40 
ALPHA! I *2)=1. 

ALPHA! I > 3 ) =1  • 

BETAI! I *2>=U. 

BETAI! I *3»=0. 

BETA2 ( I *2)=0. 

BETA2! I *3>=0. 

GAMMA! I >2)=0. 

61  GAMMA! I >3)=0. 

DO  62  1=1*10 
ALPHA! I »4>=0. 

ALPHA! I *5>=0. 

BETAI! 1 *4 ) =E0 
BETAI! I *5>=E 
BETA2! I »4)=E 
BETA2! I * 5 > =E0 
GAMMA! 1 *4 ) =0. 

62  GAMMA! 1*5) =0» 

DO  63  1=1*33 
ALPHA!  l.*6)=0. 
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BETAK  I »6 ) =E 
BETA2U*6>«E0 

63  GAMMA! I »6>=0. 

XMIN=U. 

XMAX=0. 

YMIN=0. 

YMAX=0. 

NX=0 
NY=0 
ID IM=43 
R=1.0E+U5 
NAXDIM=181 
NAYDIM=lfc 1 

call  lplacf <no »n. x. y . alpha. beTai *beta2 »gamma»ch .sen.  idim.r.tscd* 

1XMIN>*MAX*NX*YMIN»YMAX»NY>NAXDIM»NAYDIM> 

DELF=W/20. 

DELE=0.25*FACT0R 

FF1=6.2831852*E*DELF 

FF2=6»2831852*EO*DELE 

FF4=6.2831852*EO*DELF 

CHRG=0. 

DO  50  I =1 .20 

50  CHRG=CHRG+FF4*CH< 1.1) 

DO  51  1=22.41 

51  CHRG=CHRG+FF1*CH( 1*1) 

CHRG=CHRG+FF2* ( CH ( 2 1 . 1 ) +CH « 42  * 1 ) ) 

SUM=0. 

DO  64  J = 1 .NO 

I F ( ALPHA (1»J)«NE»1« ) GO  TO  64 

NJ1=N(J)-1 

DO  65  1=1. NJ1 

IP1=I+1 

SEGX=X( IP1.J)-X(I * J ) 

SEGY=Y( IPl.J)-Y(I.J) 

segarg=segx*segx+segy*segy 

SEGLEN=SQRT(SEGARG) 

65  SUM=SUM+CH( I »J)*CH( I . J ) *SEGLEN 

64  CONTINUE 
PDDSQF=FACTL*SUM 

54  CONTINUE 
RETURN 
END 

SUBROUTINE  LPLACF 1 NO *N . X .Y .ALPHA .BETA1 .BETA2 .GAMMA »CH .SCH . ID IM »R . 

i tscd»xmin*xmax.nx,ymin»ymax»ny»naxdim»naydim) 

dimension  X(  I DIM.NO ) >Y(  IDIM.NO)  .ALPHA  UD  IM  .NO  > . BE  i A1 UDIM  .NO ) . 
1BETA2<  I DIM.  NO)  .GAMMA  t IDIM.NO) .CH\ IDIM.NO)  . io<_D«  IDIM.NO. 4)  *N»NO)  » 
2A(  181.18D  .SCH (NO)  »BU8v ) »A1 » 181 . 18 i ) 

COMMON/HELP/A. A1 
EQUIVALENCE ( A » A 1 ) 

PI=3. 1415926 

RR=R*R 

DO  1 L=1»N0 

NN=N(L)~1 

DO  1 1=1. NN 

X I=X( I + 1,L)-X( I >L ) 

Y I =Y ( I+1,L)-Y( I.L) 

TSCD< I.L.1)=ATAN2(Yi,XI) 

TSCDU  >L.2»=SIN(TSCD(  I»I.*D  ) 

TSCD< I.L.3)=C0S(TSCD( I *L  * 1 ) ) 
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1 TSCD< I * L * A ) =SQRT( XI*XI+YI*YI ) 

J JJ*0 

DO  4 LJS1 »NO 

NJ=N(LJ)-1 

JAJ=JJJ 

JJJ=JJJ+NJ 

DO  4 J=1 *NJ 

JJ=JAJ+J 

XJ= ( X ( J *LJ ) +X( J+l *LJ ) ) / 2 • 

YJ*(Y(J.|_J)+Y(  J+1*LJ)  ) / 2 • 

111  = 0 

DO  4 Ll  = l *NO 
N I =N ( LI ) -1 
1 1 1 = 1 1 1 +N I 
I A 1 = 1 1 I-NI 
DO  4 I = 1 *N I 
I 1= IAI+ I 

I F C I I.EQ.JJ)  GO  TO  3 
Xl  = XJ-X( I * L I ) 

X2=XJ-X( 1 + 1 »LI ) 

Y1=YJ-Y(I.LI) 

Y2=YJ-Y( I + l.LI  ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

S1=0. 

S2=0. 

YT=Y( I *LI)*X2-Y( 1+1»LI)*X1+YJ*(X( t+l»LI)-X( I*LI) ) 

XT-X1*X2+Y1*Y2 

TETA=ATAN2(YT.XT) 

I F l ALPHA (J  »LJ ) «EQ«0« ) GO  TO  2 

51  = TSCD(  I *LI  >4)*<  l.-'J.5*ALOGiR2/'RR)  ) +0. 5*AL0G  ( R2/R1  )*  l Xl*T  ScD  * I * L 1 
1 »3 )+Yl*TSCD< I »LI *2) ) +TET  A* l Xl*TSCD< 1 » L I >2)-Yl*TScD< I *LI »3> ) 

2 TETA1  = TSCD( J.LJ»1)-TSCD< I * L I * 1 ) 

52  = 0.‘5*SlN(TETAl)*ALOG(R2/Rl)+COS(  I E I A1 ) * I E I A 
S3=-S2 

GO  TO  4 

T Sl  = TSCDll»LI»4)*U.-AL0G(TSCD(l»LI»4)/2./R)  ) 

S2=-P 1 
S3=-P  I 

4 AUJ.II  )=ALPHA(J*LJ)*Sl+BETAlU»LJ)*S2+BElA2l  J*LJ)»S3 

M=0 

DO  5 L= 1 »N0 

5 M=M+N ( L ) -1 
JJJ  = 0 

DO  6 L= 1 *N0 
NN=N(D-1 
JAJ= JJJ 
JJJ=JJJ+NN 
DO  6 J=1»NN 
JJ= JAJ+J 

6 B ( J J ) =GAMMA i J »L ) 

CALL  ARRAY(2»M*M»NAXDIM*NAYDIM>A*A) 

CALL  SIMO(A»B*M»KS) 

IF  (KS.ME.O)  PRIRT  100 
100  FORMAT  1H0  *18HSYST>IM  IS  SINGULAR) 

JJJ=0 

DO  7 L=1 »NO 
NN=N ( L ) “1 
JAJ=JJJ 
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JJJ=JJJ+NN 
DO  .7  J=1»NN 
JJ=JAJ+J 

7 CH( J*t)=B( JJ) 

DO  8 L=1 *NO 
NN=N ( L ) -1 
SCH(L>=0. 

DO  8 1=1 *NN 

8 SCh' ( L ) =SCH ( L ) +TSCD  ( I »L.A)*CHI  I *D 
IF(NX-l) 17*9.10 

9 DX=0. 

GO  TO  11 

10  DX=(XMAX-XMIN!/FLOAT(NX-l> 

U I F ( NY-1 ) 17  * 12  * 13 

12  DY=0. 

GO  TO  1A 

18  DY=(YMAX-YMIN)/FLOAT«NY-l) 

1 A DO  1-6  I 1 = 1. NX 

XJ  = XMIN+Fl OAT ( I 1-1 ) *DX 
DO  16  JJ=1.NY 
YJ=YMIN+FLOAT( JJ-1)*DY 
A ( II.JJ)=0. 

DO  16  L I =1 .NO 
NN=N ( L I ) -1 
DO  16  1=1. NN 
X1=XJ-X< I .LI ) 

X2=XJ-X( I+1.LI) 

Y1=YJ-Y( I .LI) 

Y2=YJ-Y< I+l.LI ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

IF( (Rl.EQ.O. > .OR. (R2.EO.O. ) ) GO  10  15 

YT  = Y(  I . L I ) *X2-Y  ( 1 + 1 >L  I ) *-Xl  + Yj»  l X l 1 + 1 . l I ) —A  ( I *l  I ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT.XT) 

Sl=TSCD( I »L I »A) *< l.-0.5*AL0G(R2/RR) ) +0. 5*AL0G l R2/R1 )* l XI* I i>CD l I .LI 
l*3)+Yi*T.Scn<  I »L  I » 2 ) )+TETA*lXi*i6cD>  I’Ll  »2)-U*i*CD‘  I » L I >3)  ) 

GO  TO  16 

15  Sl=TSCD( I >L I » A ) * ( l.-^.5*AL0G< (R1+R2>/RR> ) 

16  A( I I »JJ)-A( I I *JJ)+Sl*CH( I .LI ) 

17  RETURN 
END 

subroutine  timeiw.s.h.t.l.chi ) 

DIMENSION  N<3)  »X(70»3)  >Y(70.3)  .ALPHAI70.3)  .BETAH70.3)  > 
1BETA2I7U.3) . GAMMA ( 70.3) *CHl7U*3) *SCH13>  > I iCD i 7o * 3 * A ) »A» 181  * 181 ) . 
2A1 (181.181) 

REAL  L 

COMMON/HELP/A. A1 
EQUIVALENCES. Al) 

N0=3 
Nil) =A1 
N ( 2 ) =70 
N < 3 ) =70 

E0=l,/(A.*3«1A15 9* 2.99776*2.99776 E+09) 

FACTOR=( l.UE-05)*( 1./O.3937U0) 

DELW=W/AO. 

DELG=(L/2.-W/2.-S)/39. 

DELT=T / 10. 
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DELCW2=(W/2.+S)/20. 

X ( 1 » 1 ) =-W/2 • 

Y(1»1J=0. 

DO  26  I=1»A0 
IP1=I+1 

X( IP1»1)=X( I.ll+DELW 
Y ! I PI » 1 > =Y( I *1) 

ALPHA! I »1)=1. 

BETA1 ( I » 1 ) =0» 

BETA2 ( I * 1 ) =0. 

26  GAMMA ( I * 1 ) = 1 • 

X ( 1 *2 ) =-L/2 • 

Y ( 1 »2 ) =0. 

X(l»3)=L/2. 

Yl  1 *3) =0« 

X ( AO .2 ) =-W/2«“S 
Y(A0.2)=0. 

X!  50.21=X(A0.2) 

Y( 50  *2 ) =-T 

X(A0*3)=W/2.+S 

Y!A0.3)=0. 

X( 50.3)=X(A0,3) 

Y ( 50  *3 ) =-l 
00  27  1=1.38 
IP1=I+1 

X( IP1 *2 )=X( I *2 ) +DELG 
Y ( IP1 .2 ) =Y( 1 .2) 

X ( IP1 »3)=X( I .31-DELG 

27  Y ( I PI »3 ) =Y( I .3) 

DO  227  1=1.9 
IPAO= I+AO 

I PA0M1= I P40-1 
X( IP40»2)=X( IPAUM1 »2 ) 

Y ( IPA0,2)=Y< IPA^Ml .21-DELT 
X(  I PAU  » 3 ) =X ( IPAOM1.3) 

227  Y( I PAO . 3 ) =Y ( I PAUM1 »3 ) -DELT 
DO  228  1=1.20 

IP50=I+50 
I P50M1= I P50-1 

X( IP50.2>=X(IP50M1.2)+DELCW2 
Y! IP50.2)=Y( IP50M1.2S 
X( IP50,3)=X( IP50M1.3)-DELCW2 

228  Y( IP5U»3)=Y( IP50M1.3) 

DO  229  1=1,69 
ALPHAt I .2  > = 1 • 

BETAK  I .2  5 =0, 

BETA2 ( 1 »2)=u. 

GAMMA ( I .2 ) =0« 

ALPHA1I »3)=1. 

BETA1U  >3>=0, 

BEiA2( I ,3)=0. 

229  GAMMA! I * 3 > =0. 

XMIN=0. 

XMAX=0. 

YMIN=0. 

YMAX=0. 

NX=0 

NY=0 

IDIM=70 
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I 


c 


tt 

[ 


0 


r 

|; 


R=1.0E+05 

NAXD1M=181 

CALL  LPLACE«NO»N.X.Y»ALPHAtBETAl.BETA2.GAMMA.CHt!.CH.IDIM »R*  I 5CD* 
1XMIN.XMAX.NX.YMIN.VMAX.NY .NAXDIM.NAYDIM) 

CH1=SCH(1»*6.2831852*EU 

RETURN 

END 

SUBROUTINE  LPLACEt NO.N.X.Y .ALPHA.BE  I A1 .BE  I A2 .GAMMA.CH.SCH. IDIM.R. 

1 TSCD.XMIN.XMAX.NX.YMIN.YMAX.NY.NAXDIM.NAYD1M) 

DIMENSION  XI IDIM.NO) »Y( IOIM.NO) .ALPHA I IDIM.NO) .BEIA1I ID1M.NO) » 
1BETA2(IDIM.N0) .GAMMA! IDIM.NO! »CH( IDIM.NO) . ISCD! IDIM.NO. A.)  .NINO) . 
2AI 181.181 ) *SCH< NO) »B( 180) .All  18 1.181) 

COMMON/HELP/ A.A1 
EQUIVALENCE! A. Al) 

PI=3.1A15926 

RR=R*R 

DO  1 L=l.NO 

NN=N(L)“1 

DO  1 1 = 1. NN 

X I=X! I+1.L1-X! I .L) 

YI=Y!I+1.L)“Y!I.L) 

TSCD ! I .L  » 1 ) =ATAN2 « Y I .XI ) 

TSCD!  I «L » 2 ) =SIN ! T SCO < I *1. * 1 ) ) 

TSCD!  I.L.31=COS!TSCD!  I .1. *1 ) > 

1 TSCD! I»L.A)=SQRT!XI*XI+YI*YI) 

JJJ=0 

DO  A I — 1 = 1 .NO 
NJ=N!LJ)-1 
JAJ=JJJ 
JJJ=JJJ+NJ 
03  A J=1.NJ 


J J= JAJ+J 

X J = ( X t J »L«* ) +X  ( J4l » LJ  ) )/2. 
YJ=(Y!  J»(.J)+Y«J+1.LJ)  ) / 2 • 
I ! 1=0 

DO  A L I =1 .NO 
N I =N  ! LI  ) -1 
III=III+NI 
I A 1= 1 1 1 ~N I 
DO  A 1 = 1 . N I 
I I = I A 1 4-1 

IF!  I I.EQ.JJ)  GO  TO  3 
X1=XJ-X! I .LI ) 

X2=XJ-X« 1+1 »LI ) 

Y1=YJ-Y( I .LI ) 

Y2=YJ-Y ! 1+1 .LI ) 

Rl  = Xl*Xim*Yl 
R2  = X2*X2-»  i'2*Y2 
S1=0. 


c 2 ^ o i 

YT=Yt I .LI )*X2-Y( 1+1 »LI )*Xl+YJ* (XI 1+1 »LI)-AlI»LI)> 

XT=X1*X2+Y1*Y2 
TETa=ATAN2!YT»X1  ) 

I F (ALPHA (J»LJ).EQ«U«)  GO  TO  2 

Sl-TSCD  ( I >LI.A)*(l.-0.5*ALOG(R2/RR)  )+0. 5*ALOG  ! K2/K-1  )*  I XI*  I *CD  1 1 »L  i 

lt3)+Yi#TSCDlI.L1.2))+TETA*«Xl*T5»CD(I»LI»2>-Yl*liCDiI»LI.3)) 

2 TeTa1=TScD«J»LJ.1)-TScD( I.LI.l) 

'S2=0.3*SIN(TETA1)*AL0G(R2/RI)+C0S! I E I Al » * • E I A 
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S3=-S2 
GO  TO  4 

3 Sl  = TSCD(l»LI*4)»(l.-AL0G(TSCDtl»LI  »/>|/2»/R)t 
S2=-PI 

S3=~P I 

4 A ( J J . 1 1)  =ALPHA  ( J *L  J )>S1+BETA1  (J  »LJ ) *S2+BETA2  ( J ,LJ ) *S3 
M=0 


00  5 L=1 *N0 
5 M=M+N(L.)-1 
JJJ=0 


DO  6 L=1 *N0 
NN=N(L)~1 
JAJ=JJJ 
JJJaJJJ+NN 
DO  6 J=1 »NN 
JJ= JAJ+-J 

6 B( JJ)=GAMMA(J»L) 

CALL  ARRAY(2*M»M»NAXDIM»NAY0IM»A»A) 
call  S1MQ<A,3*M»KS) 

IF  (KS.NE.O)  PRINT  10b 
100  FORMAT (1  HO »18HSYSTEH  IS  SINGULAR) 
JJJ=0 

DO  7 L-liNO 
NN=N(L)-1 
JAJ= JJJ 

jjj-jjj+nn 

DO  7 J-l »NN 
J J= JAJ+ J 

7 CH( J »L 5 =B I J J ) 

DO  8 L=1.N0 
NN=N ( L ) -1 
SCH(L)=0. 

DO  8 1=1, NN 

8 SCH 1 1, ) =SCH(  L ) +TSCD  I I »L»4)*CKl  I ,L) 
IFtNX-l)17,9,10 

9 DX=0  • 

GO  TO  11 

10  DX=(XMAX-XMIN) /FLOAT (NX-1) 

11  I F ( NY~1 )17»12>13 


12  DY=0. 

GO  TO  14 

13  DY=(VMAX-YMIN)/FLOAT(NY-l) 

14  DO  16  II=1*NX 

XJ  = XMIN+FLOAT( I 1-1 ) *DX 
DO  16  JJ=1,NY 
YJ-YMIN+FLOAT (JJ-1 )*DY 
A! I I ,JJ)=0. 

DO  16  L I =1 »N0 
NN=N(LI ) -1 
DO  16  I "1 »NN 
X1=XJ-X( I » L I ) 

X2=XJ-X( I+ltLI ) 

Y1  = YJ-Y ( I » L I ) 

Y2=YJ-Y ( 1+1 »LI ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

I F ( ( R 1 . EQ»0 • ) • OK . ( K2 • Lu» 0 . ) ) bu  lu 

YT=Y(  1 ,LI  )*X2-Y(  I+1,lI  )*Xl+Yj-  ( XU+1  >t-I ) -A  i 1 , u I) ) 
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XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT»XTJ 

Sl  = TSCD< I »LI *AJ*( l.-0.5»AL0G(R2/RR) )+0.5*ALOG(R2/Rl>*«Xl*T6COU  »LI 
l*3)+Yl»TSCD<  I*LI*2)  ) + TETA*(Xl*TScD(  I *LI  *2>-Yl#li»CDl  I*LI  »3U 
GO  TO  16 

15  Sl=TSC0( I »LI»4)*(1.-0.5*AL0G( (R1+R21/RR)) 

16  A(  1 1 »JJ)=A( 1 1 *JJ)+Sl*CH( I »Ll ) 

17  RETURN 
END 

SUBROUTINE  ARRAY  (MODEM  *J*N»M»S*D) 

DIMENSION  SCI) »D( 1 ) 

NI=M-1 

I F (MODE-1 ) 10O.1OO.12U 
10o  IJ=I*J+1 
NM=N*J+1 
DO  11U  K=1,J 
NM=NM-N I 
DO  110  L=1 » I 
IJ=IJ-1 
NM=NM-1 

110  D(NM)=S( I J ) 

GO  TO  140 
120  IJ=0 
NM=0 

DO  130  K=1 , J 
DO  125  L=1 . 1 
IJ=IJ+l 
NM=NM+1 

125  S( I J!=D(NM) 

130  NM=NM+NI 
lAo  RETURN 
END 

SUBROUTINE  SIMO(A,B,N»KS) 

DIMENSION  A ( 1 ) »B ( 1 ) 

T0L=0.0 

KS=0 

JJ=-N 

DO  65  J=1»N 
JY=J+1 
JJ=JJ+N+1 
BlGA=0. 

I T = JJ-J 
DO  30  I=J»N 
I J=IT+I 

if(abs(biga)-abs(a<  urn  20.30,30 
20  biga=a;ij) 
imax=i 

30  CONTINUE 

IF(ABS(BIGA)-T0L)  35,35,40 
35  KS  = 1 
RETURN 

40  Il=J*N*(J-2) 

IT=IMAX-J 
DO  50  K=J»N 
1 1=1 1+N 
12=1 1+ I T 
SAVE=A< II) 

A(  I1')=A(  12) 

A( I2)=SAVE 
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50  AtIl)=A(  IU/BIGA 
SAVE=B« IMAX) 

B( IMAX)=B(J) 

B( J)=SAVE/BIGA 
IF(J-N)  55  »70 155 
55  IQS=N*(J-1) 

DO  65  I X=JY  »N 
IXJ= IQS+IX 
IT=J-IX 
DO  60  JX=JY»N 
IXJX=N*( JX-ll+IX 
JJX= IXJX+IT 

60  A<  I X JX ) = A (I X JX ) - ( A ( IXJ)^A(JJX) ) 
65  BUX)=B!  IXjt— <B<  J)*AC  IXJ)  ) 

70  NY-N-1 
IT=N*N 

DO  80  J=1»NY 
I A= I T-J 
IB=N-J 
IC=N 

DO  80  K=1  * J 

B( I B ) =B ( IB}-A(IA>«B< IC) 

I A= IA-N 
80  IC=IC-1 
RETURN 
END 


78 


SI 


t*r***&«  jn*w*frh,y»Wi—  y y „ 


Appendix  E 

LISTING  OF  PROGRAM  LEMDOC 

PROGRAM  LEMDOC 
REAL  L 

READ  1*  NSETS 
1 FORMAT  ( 1 10 ) 

DO  2 MN=1*NSETS 

READ  3»W»S»T»H»L*ER»TAND1»SIGMA 

3 FORMAT ( 8FIO. 6 ) 

TANDl=TANDl*l.E-04 
S I GMA=S I GMA* 1 • E+0  7 
PRINT  8 

8 FORMAT  ( 1H  t////lH  »*STRIP  METALUATION  THICKNESS  = 250  MICROINCHES 
1*) 

PRINT  4*W*S»T»H»L*ER»TAND1  *SI'GMA 

4 FORMAT (1H0*2HW= ^5. 1 » 3X. 2HS= ?F5. 1 »3X»2HT=» 
lF5.1»3X'.2HH=tF5.1»3X*2HL=»F10.6»3X»3HER=»F10.6*3X» 
26HTANDI=»E1U.3»3Xj6HSIGMA=*E10.3) 

FACTOR=! 1.0E-05)*! 1./0.3937UO) 

W=W*FACTOR 
S=S*FACTOR 
T=T*FACTOR 
H=H*FACTOR 
L=L*FACTOR*< l.E+03) 

E1=ER 

E2=ER+G.01*ER 

CALL  MSCUPF(W,S»T*H»L>E1 *SIGMA  »C0  »CE  *PDDSQO  »PDDSQE ) 

PDDSF0=PDDSQ0 

PDD$FE=PDDSQE 

CALL  MSCUPF(W»S»T,H»L»E2  >SI  GMA  tCOP  *CEP  >PDDS(JO  >PDDSUE ) 

CALL  MSCUPE ( to  >S  »T  »H*L  > 1 1 *CtO»CEE ) 

AE=CE*CEE 

AO=CO*CEO 

BF-’CEE/CE 

BO-CEO/CO 

Z0F  = 1./ I ( 3<OE+u8 ) #SQRT ( At ) ) 

ZOOr-i./  ( (3.0E+o8)*SQRT<  A0> ) 

AA-ZOLvZOO 

ZO=SQRT(AA) 

VE=3.0E+08*SORT (BE) 

VO=3.0E+08»SQRT(BO) 

VAVG= ( VE+VO ) / 2 . 

RHO=ZOE/ZOO 

C= ( RHO-1 « ) / ( RHO+1 • ) 

CDB--20*  X-ALOGIO  ( C ) 

ARGlt=BE 

ARG10=80 

ARG2E=CEE/CEP 

ARG20=CE0/30P 

ERXE=1./ARG1E 

ER10=1./ARG10 

ER2E=1./ARG2E 

ER20=1./ARG20 

FACTC=10./2.3 

F AC TD= 27.3/3* E+0 8 

ACE=FACTC*PDDSFE*ERlE*ZOE 

ACO=FACTC*PDDSFO*ER10*ZOO 

A0E=F^CT0*(E1/SQRT(ER1E) )*( ( ER2E-ER1E ) / ( E2-E1 ) >*TAND1 
AD0=FACTd*<E1/SQRT(ER10> !*( (ER20-ER10)/ (E2-E1 ) )*TAND\ 

PRINT  5 

5 FORMAT ( 1hO*1X»5HC<DB) » 3X » 8HZ0 ( OHMS ) * 2X»9HZOO< OHMS ) » 
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19H^0E«0HMS) *5X»9HVO(M/SEC) »5X * 9HVE (M/SEC ) »5X* 

2 1 1HVAV6 (M/SEC ) ) 

PRINT  6 *CDB  *20 .ZOO .ZOE  »V0  * VE  * VAVG 

6 FORMATIlH  *F5.2*5X»F6.2»5X.F6.2»5X*F6.2*5X*E10.3»5X.E10.3» 
15X.E10.3) 

PRINT  7 

7 FORMAT! 1H<J»20HALFAC0(DB/MJ/3QKT(F) ♦ 5X » 20HALFACE (DB/M) /SORT t F ) »5X, 
1 14HALFADO( DB/M ) /F  * 5X  * 14HALFADE ( DB/M ) /F ) 

PRINT  8 * ACO.ACE  .ADO  *ADE 

8 FORMAT ( 1H  . 5X »E1U . 3 . 15X , ElO. 3 > 12X *E10 . 3 >9X »E 10. 3 ) 

2 CONTINUE 

STOP 

ENO 

SUbROUT  i ne  MScUPF  ( v.  ,S  » T ,h  »L » ER » S I GMA  *C0  »CE  » PDDSQO » PDDSQE  ) 
DIMENSION  N<6) »X(46.6) *Y(46»6)  *ALPHAU6>6)  »BETA1  ( 46 *6 ) . 
1BETA2<A6*6)  » GAMMA ( 4 6 » 6 ) »CH(46»6)  »SCH(6)  >TSCD<  46  r-6  »4  ) *At  131  * 181 ) » 
2A1( 132.132) 

REAL  L 

common/help/a. A1 

EQU I VALENCE ( A.A1 ) 

NO=6 
N ( 1 ) =43 
N ( 2 ) =43 
N(3) =46 
N ( 4 ) = 13 
N ( 5 ) =13 
N ( 6 ) =23 

EC=I,/ < 4. *3. 141 59*2. 99776*2. yy 776L+09) 

AAA=3 .141 59*4. t-07*3, 141 39/SJOMA 
RDSQF=SQRT ( AAA ) 

FACTN=4. *3. 14159*3. 141 59*RDSQF*0 . 85E-1 2 
FACTD=4. *3. 14159*1. E-07 
FACTL=FACTN/t-ACTD 
FACTOR=(  1.0E-U5)*!  1./v^.3937vjO) 

DE'-W=W/20. 

E=ER*E0 

X( 1 »l)=-(S/2. )-W 
Y ( 1 » 1 ) =H 
X ( 1 »2 ) =-X (1*1) 

Y ( 1 * 2 ) =H 
DO  26  1=1,20 
IP1=I+1 

X( IPl.l )=X( I » 1 ) +DELW 
Y ( I PI .1 ) =H 

Xt I PI .2 ) =X( 1*2 '-DELW 

26  Y ( I PI >2 ) =H 

X(  22  *1 ) =X ( 2 1 » 1 ) 

Y ( 22  » 1 ) =0, 25*FACT0R+H 
X( 22  *2 ) =X ( 21 .2 ) 

Y(22*2)=Y(22»1) 

DO  27  1=22.41 
IP1=I+1 

X( IP1 *1 )=X( I ,1)-DELW 
Y( IPl.l )=Y( I ,1) 

X< I PI . 2 ) =X ( I »2 ) +DELW 

27  Y ( I PI .2 ) =Y( I .2) 

X ( 43 . 1 ) =X( 1 . 1 ) 

Y ( 43  » 1 ) =Y ( 1 . 1 ) 

X ( 43 .2 ) =X 1 1 . 2 > 
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Y ( 43  >2 ) =Y ! 1 ♦ 2 ) 

DELL=L/45. 

X ( 1 » 3 ) =-L/2 • 

Y( 1 ♦3)=U. 

00  28  1=1.45 
IP1=I+1 

X(  I PI .3 ) =X ( I .3 ) +DELL 

28  Y ( I PI .3 ) =0. 

DELH=H/2. 

X ( 1 .4 ) =X ( 1 »3 ) 

Y( 1,41=0. 

X<  >.5)=L/2. 

Y(1.5)=0. 

DO  29  1=1.2 
IP1=I+1 

X( I PI ,4 1 =X ( 1 .4 ) 

Y ( IP1 .4 ) =Y( I .4 ) +DELH 
X( I PI .5 ) =X ( 1 .5  ) 

29  Y ( I PI »5 ) =Y ( I .5 l+DELH 
DEL=(L/2.-<S/2.+W) >/1u. 

DO  30  1=3.12 

IP1=I+1 

X( I pi ,4 ) =X ( I .41+DEL 
Y ( I PI .4 ) =H 
X ( I PI . 5 ) =X ( I .51-DEL 

30  Y ( I P 1 » 5 ) =H 
X(l»6)=-S/2.-W 

Y l 1 .6 1 =H+0. 2 5* FAC TOR 
X(21*6>=S/2.+W 
Y ( 2 1 » 6 > = Y ( 1 » 6 ) +T 
DELT  = T / 2 « 

DO  31  1=1.2 
1P1=1+1 
I P2 1= 1+21 
1P21K1=IP21~1 
X( I PI .6 ) =X ( I .6) 

Y ( I PI .6 ) =Y ( I .61+DELT 
X( I P21 » 6 ) =X ( IP21M1  .6) 

31  Y( IP21»6)=Y( IP21M1 .61-DELT 
DELOV=(2*W+S)/18. 

DO  32  1=3.20 
IP1=I+1 

X ( IP1 >6 ) =X ( I .61+DELOV 

32  Y ( I PI .6 1 =Y ( I *6) 

DO  47  LLL--1.2 
IF(KKK-O)  38.37.38 

C ODD  MODE 
37  CONTINUE 

DO  332  J=1 .2 
DO  33  1=1.42 
ALPHA ( I »J)=1. 

BETAK  I » J)  =0. 

BETA2U  »J)=0. 

GO  TO  (34.35) .J 

34  GAMMA ( I » J ) = 1 • 

GO  TO  33 

35  GAMMA ( I »J>=-1. 

33  CONTINUE 
332  CONTINUE 
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KKK=KKK+1 
GO  TO  39 
C EVEN  MODE 

38  CONTINUE 

DO  40  J=l»2 
DO  41  1=1.42 
ALPHA! I »J)=1. 

BETAK  I »J)=0. 

BETA2 ( 1 * J) =0. 

42  GAMMA! I * J>  = 1« 

41  CONTINUE 
40  CONTINUE 

39  CONTINUE 

DO  36  1=1.45 
ALPHA! I .3>=1. 

BETAl! I .3>=0. 

BETA21 I »3>=0. 

36  GAMMA! I ,3»=0. 

DO  45  1=1.12 
ALPHA! I »4 ) =0* 

BETAl! I . 4 ) =E 
BETA2! I .4 ) =EO 
GAMMA! I .4) =C. 

ALPHA! I . 5 > =0. 

BETAl! I »5>=Eu 
BET  A2 ( I . 5 1 =E 

45  GAMMA! I .5)=U. 

DO  46  1=1.22 
ALPHA! I .6>=0, 

BETAl! I .6)=E 
BETA2! I » 6 ) =E0 

46  GAMMA! I >6)=0. 

XMIN=0. 

XMAX=o. 

YMIN=0. 

YMAX=0. 

NX  = 0 
NY  = 0 
I D IM=46 
R=1  •0E+'-<5 
NAXDIM=181 
N AYD I M= 181 

CALL  LPLACF (NO.N.X.Y.ALPHA.BETA1.BETA2 .GAMMA .CH.bCH. IDiM.R.TSCD. 
1XMIN.XMAX.NX.YMIN.YMAX.NY.NAXDIM.NAYDIM! 

DELF=W/20. 

DELE=0.25*FACTOR 

FF1=6«2831852*E*DELF 

FF2=6.2831852*E*DELE 

FF4=6.2831852*EO*DELt 

CHRG=0. 

DO  50  1=1, 2U 

50  CHRG=CHRG+FF1*CH< I »1) 

DO  51  1=22.41 

51  CHRG=CHRG+FF1*CH< I ,1) 

CHRG=CHRG+(FF2*CH<  21,1 )+FF4*CH( 42.1 > ) 

SUM=0. 

DO  64  J=l»NO 

IF (ALPHA! 1»J)«NE«1.)  GO  TO  64 
NJ1=N( J)-l 
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DO  65  1=1 »NJ1 
IP1=I+1 

SEGX=X( IPl.Jl-XI I.J) 

SEGY=Y( IP1»J)-Y( I *J) 

segarg=segx*segx+segy*segy 

SEGLEN=SQRT(SEGARG) 

65  SUM=SUM+CH< I *J)*CH< I *J)*SeGLEN 
64  CONTINUE 

PDDSQF=FACTL*SUM 
I F ( GAMMA  1 1 * 2 ) - 1 • ) 1001 , 1000 » 1001 

1000  PDDSQO=PDDSQF 
CO=CHRG 

GO  TO  47 

1001  PDDSQE=PDDSQF 
CE=CHRG 

47  CONTINUE 
54  CONTINUE 
RETURN 
END 

SUBROUTINE  LPLACF ( NO*N  »X » Y .ALPHA  * BETA l *BETA2  > GAMMA »CH»SCH* IDIM.R » 
i tscd.xmin.xmax.nx.ymin.ymax  *NY»NAXDIM*NAYDIM) 

DIMENSION  X( IDIM.NO) »Y( IDIM.NO) » ALPHA! IDIM.NO) .BETaK I DIM* NO) * 
1BETA2!  IDIM»N0)  .GAMMA < IDIM.NO)  »CH<  IDIM.NO) »TSCD<  IDIM*N0»4)  »‘N(NO) » 
2A( 181 *181 ) »SCH(N0) *B  < 1 75 ) .AK  132  *132) 
common/help/a.ai 

EQUIVALENCE  C A * A 1 ) 

F 1=3.1415926 

RR=R*R 

DO  1 L= 1 »N0 

NN=N ( L ) -1 

DO  1 1=1 »NN 

XI=X( I+1*L)-X( I.L) 

Y 1 - Y ( I + 1,L)“Y( I .L) 

T SCO ( I *L  * 1 ) =ATAN2 ( Y I » X I ) 

T SCD ( I »L  » 2 ) =S IN ! T SCD ( I *L»1 ) > 

TSCDI I » L * 3 > =COS  t T SCD l I »L  » 1 > > 

1 TSCDt I»L»4)=SQRT(XI#XI+YI*Y1 ) 

JJJ=0 

DO  4 I 1=1  >N0 

N J=N ( L J ) -1 
JAJ=JJJ 
JJ J= JJJ+NJ 
DO  4 J= 1 >NJ 
J J= JAJ+ J 

XJ=(X(J»LJ)+XU+1  »LJ)  ) / 2 • 

YJ=(Y(J»LJ)+Y(J+1.LJ) )/2. 

I 11=0 

DO  4 Ll=l *NO 
N I =N ( LI )-l 
1 1 1 = 1 1 1 +N I 
I A I = 1 1 1 —N I 
DO  4 1=1 >NI 
I I = I A 1+  I 

IF(II.EQ.JJ)  GO  TO  3 
X1=XJ-X( I »L I ) 

X2=XJ-X ( 1+1 *LI ) 

Y 1= YJ-Y ( I »L I ) 

Y2=YJ-Y ( 1+1 *LI ) 

R1=X1*X1+Y1*Y1 
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R2=X2*X2+Y2#Y2 

S1=0. 

S2=0. 

YT=Y( I»LI )*X2-Y( I+l.LI )*Xl+YJ*(X( I+ltLI } —X ( I *LI ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT»XT) 

IF ( ALPHA (J*LJ) . EQ.O. ) GO  TO  2 

Sl=TSCD< I »LI *4)*( l.-0.5*ALOG(R2/RR> ) +0. 5*AL0G ( R2/R1 ) * t X1*TSCD < I »LI 
1 »3)+Yl*TSC0( I »LI»2> )+TETA*(Xl*TSCO( I *LI *2 ) -Y1*TSCD t I *LI  *3) > 

2 TETA1=TSCD< J*LJ»1)-TSCD( 1*LI»1) 
S2=0.5*SIN(TETA1)*AL0G<R2/R1)+C0S(TETA1)*TETA 
S3=-S2 

GO  TO  A 

3 Sl=TSCD< I >LI *A)*«l.-ALOG(TSCDl I *L I >4 I /2 ./R ) ) 

S2=-PI 

S3=-PI 

U A ( J J * 1 1 ) =ALPHA I J »L J ) *S 1+BETA1 < J »L  J ) *S2+BET A2 ( J ♦ L J ) *S3 
M=0 

DO  5 L=1.N0 

5 M=M+N(L)-1 
JJJ  = 0 

DO  6 L= 1 »N0 
NN“N(L)-1 
JAJ=JJJ 
JJJ=JJJ+NN 
DO  6 J= 1 *NN 
J J-JAJ+J 

6 B< JJ)=GAMMAI J.L) 

call  array(2»m.m,naxdim*naydik»a»a) 

CALL  SIMQ(A»B*M*KS) 

IF  (KS.NE.O)  PRINT  luv 
100  FORMAT ( 1H0.18HSYSTEM  IS  SINGULAR) 

JJJ  = 0 

DO  7 L = 1 .NO 
NN=N I L ) -1 
JAJ= JJ J 
jjj=jjj+nn 
DO  7 J= 1 »NN 
JJ= JAJ+J 

7 CHI J.L)=6< JJ) 

DO  8 L= 1 *N0 
NN=N(L)-1 
SCH ( L ) =0 , 

DO  8 1 = 1. NN 

8 SCH(L)=SCHIL)+TSCD( I .L.4)*CH( 1 » L ) 

I F ( NX-1 ) 17 .9 » 10 

9 DX  = 0. 

GO  TO  11 

10  DX= ( XMAX-XM IN i /FLOAT ( NX-1 ) 

11  IF(NY-1)17.12.13 

12  DY=0. 

GO  TO  14 

13  DY=(YMAX-YMIN)/FL0ATINY-1) 

14  DO  16  1 1 = 1 .NX 

XJ  = XMIN+FL0AT ( 1 1 - 1 ) *DX 
DO  16  JJ=1.NY 
YJ=YMIN+FL0AT( JJ-1)*DY 
AIII . JJ ) =0. 
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DO  16  L I =1 .NO 
NN=N ( L I ) -1 
DO  16  1=1 >NN 
X1=XJ-X ( I »LI ) 

X2=XJ-X ( 1+1 .LI > 

Yl  = YJ-Y ( I * L I ) 

Y2=YJ-Y( I+l.LI ) 

R1=X1*X1+Y1*Y1 

R2=X2*X2+Y2*Y2 

I F ( (Rl.EQ.'J.  ) .OR.  (R2.EQ.0.)  ) GO  TO  15 

YT=Y( I ♦L1)*X2-Y< I+l.LI )*X1+YJ*(X( 1+1 *LI ) —X ( I • L I ) ) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2<YT,XT) 

Si  = TSCD( I »LI »4)*( l.-0.5*AL0G<R2/RR) >+0. 5*ALOG ( R2/R1 > •> ( Xl*TSCD U * L I 
1»3>+Y1*TSCD< I *L1 .2) ) + TETA*Ul*TSCD( I .LI »2)-Yl*TSCD( I *LI >3 ) ) 

GO  TO  16 

15  S1=TSCD( I *L  I .41*1 l.-0.5*AL0G< (R1+R2)/RR> ) 

16  A(  1 1 .JJ)=A< 1 1 »JJ)+S1*CH< I ,L1) 

17  RETURN 
END 

SUBROUTINE  mscupeiw.s.t.h.l.er.ceo.cee) 

DIMENSION  N<3)  »X(46*3)  »Y(46>3)  .ALPHA (46. 3)  >BETA1(46>3)  * 
1BETA2U6*3)  * GAMMA  t 46  »3 ) »CHU6*3)  »SCH<3>  > TSCD ( 46 *3 *4 ) *AU81»181)  » 
2A1  ( 132  * 132 ) 

COMMON/ HELP/ A.A1 
EQUIVALENCE  (A. All 
REAL  L 

EO=l. /< 4. *3. 14159*2. 99776*2. yV776E+Oy) 

FACTOR=(1.UE~05)-M  l./u. 3937UO) 

NO=3 
N ( 1 ) =43 
N ( 2 ) =43 
N ( 3 ) = 46 
KKK  = 0 

DEI W=W/20. 

E=ER*EO 

XI 1 ,l)=-(S/2.)-W 
Y ( 1 » 1 ) =H 
X ( 1 *2 ) =**X( 1 » 1 ) 

Y( 1 4?) =H 
DO  26  1 = 1 .20 
IP1=I+1 

XI  I PI  * 1 ) =X ( I » 1 ) +DELW 
Y ( I PI » 1 ) =H 

XI  I PI >2 ) =X ( I *2 ) -DELW 

26  Y ( J PI >2 ) =H 

X ( 22 » 1 ) =X ( 2 1 * 1 ) 

Y ( 22  * 1 > = 0.2  5*F ACTOR+H 
X(22»2)=X(21>2) 

Y(22.2)=Y<22»1) 

DO  27  1=22.41 
IP1=I+1 

XI IP1»1)=X( I » 1 ) -DELW 
Y ( I PI > 1 ) =Y ( 1*1) 

XI 1P1 »2 ) =X( I.21+DELW 

27  Y I IP1 .2 ) =Y( I .2) 

X ! 43  > 1 ) =X ( 1 » 1 ) 

Y(43.1)=Y(1  .1) 

X I 43  * 2 ! =X 1 1 . 2 ) 
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YC43»2)=Y(1.2) 

DELL=L/45. 

X( 1.3)=-L/2. 

Y(1*3)=U. 

DO  28  1=1*45 
IP1=I+1 

X(  IP1 *3 ) =X( I *3 ) +DELL 
28  Y ( IP1*3)=0. 

DO  47  LLL=1.2 
IF(KKK-U)  38.37*38 
C ODD  MODE 

37  CONTINUE 

DO  332  J=1 . 2 
DO  33  1=1.42 
ALPHA ( I »J)=1. 

BETA1(I*J)=U. 

BETA2 ( I . J ) =0. 

GO  TO  (34*35) *J 

34  GAMMA ( I » J ) = 1 • 

GO  TO  33 

35  GAMMA! I ,J)=-1. 

33  CONTINUE 

332  CONTINUE 
KKK=KKK+1 
GO  TO  39 
C EVEN  MODE 

38  CONTINUE 

DO  40  J= 1 *2 
DO  41  1=1*42 
ALPHA! I , J ) = 1 . 

8ETA1! I »J)=0. 

BETA2! I *J)=U. 

42  GAMMA! I * J ) = 1 . 

41  CONTINUE 
40  CONTINUE 

39  CONTINUE 

DO  36  1 = ] .45 
ALPHA! I * 3 ) = 1 • 

BETAl! I »3)=U. 

BETA2 ( I . 3 • =0. 

36  GAMMA! I ,3)=0. 

XM!N=Ot 

XMAX=0. 

YMIN=0. 

YMAX=U. 

NX=0 
NY=0 
I D IM=46 
R= 1 . 0E+U5 

NAXDIM=132 

NAYDIM=132 

CALL  LPLACE  !N0*N»X*Y*  ALPHA*  WEl  A1  > BET  A2*  GAMMA  *CH*i>CH.»  IDlM*K*Ti>CD* 

ixmin»xmax*nx*ymin*ymax*ny,naxdim*naydim) 

IF  (LLL.EQ.2)  GO  TO  5U 

CEO= 2**3.14159*8. 8 55E-12*SCH( 1 ) 

GO  TO  47 

50  CEE=2.*3.14159*8.855E-12*SCH! 1 ) 

4.7  CONTINUE 
RETURN 
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END 

SUBROUTINE  LPLACE(NO>N*X*y»ALPHA»BETAl»BETA2*GAMMA»CH*SCH* idim»r» 
lTSCDiXMIN»XMAXtNX*YMIN»YMAX»NY,NAXOIM*NAYDIM) 

DIMENSION  XUDIM*NO)  *Y(  IDIM»NO)  * ALPHA!  IDIM*NO)  »BETA1  ( IDIMtNO ) t 
1BETA2UDIM»N0) »GAMMAl IDIMjNO) *CH( IDIM.NO) »TSCD( IDIM»N0»4) »NlNO) ♦ 
2A1 1 132  > 132 ) »SCH( NO ) »B«129) ,A<181*181) 

COMMON/ HELP/ A *A1 
EQUIVALENCE  (A,A1> 

PI=3. 1415926 

RR=R#R 

DO  1 L= 1 »NO 

NN=N(L)-1 

DO  1 1=1, NN 

XI=X( I+1*L)-X( 1 »L) 

Y I=Y( 1+1 ,L)“Y ( I ,L ) 

TSCDI I »L,1)=ATAN2(YI »XI  ) 

TSCD( I ,L,2)=SIN(TSCD(I ,L*1) ) 

TSCDI I , L ,3 ) =COS I T SCD I I ,L , 1 ) ) 

1 TSCDI I » L * 4 ) =SQRT I X I #X I + Y 1 * Y I ) 

JJJ=0 

DO  4 LJ=l»NO 
N J=N ( LJ ) -1 
JAJ= JJJ 
JJJ= JJJ+NJ 
DO  4 J= 1 »NJ 
JJ= JAJ+ J 

XJ=(XIJ»LJ)+X(J+1,LJ) )/2. 

YJ=< Y(J»LJ)+Y(J+1,LJ) )/2. 

1 1 1=0 

DO  4 LI =1 ,NO 
N I =N I L I ) -1 
1 1 1 = 1 1 1 +N I 
I A I = 1 1 1 — N I 
DO  4 1=1, NI 
I I = I A 1+  I 

IF!  II.EQ.JJ)  GO  TO  3 
X1=XJ-X ( I ,L I ) 

X2=XJ-X ( 1+1 ,LI ) 

Y1  = YJ-Y I I ,LI  ) 

Y2  = YJ-Y I 1 + 1 »LI ) 

R1=X1*X1+Y1*Y1 
■ R2=X2*X2+Y2#Y2 
51=0. 

S2=0  • 

YT=Y(I ,LI )*X2-YI I+1*LI )*X1+YJ*IXI I +1 ,L I ) -X ( I , L I)) 

XT=X1*X2+Y1*Y2 

TETA=ATAN2(YT*XT) 

I F I ALPHA  I J ,1 I ) .EQ.O • ) GO  TO  2 

Sl=TSCDl I »LI »4)*( l.-0.5*ALOG(R2/RR) ) +0. 5*ALOG I R2/R1 }* I Xl*TSCD I I » L I 
1»3>+Y1*TSCDU  »LI ,2) )+TETA*<Xl*TSCDl I »LI *2 )-Yl*TSCDt 1 ,LI *3) ) 

? TETa1=TScD(J»LJ»1) -TSCDI I »LI ,1) 

S2=0.5*SIN(TETA1)*AL0G(R2/Rl)+C0S<TETA1>*TETA 

S3=-S2 

GO  TO  4 

3 Sl=TSCD< I »LI »4)*( l.-ALOG(TSCD< I »LI »4)/2./R) > 

S2=-P I 

S3=-P I 

4 A1 I JJ»I I )=ALPHA( J»LJ)*Si+BETA1( J*LJ)*S2+BETA2 I J»LJ)*S3 
M=0 
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DO  5 L=l.NO 

5 M=M+N ( L ) -1 
JJJ=0 

DO  6 L=l.NO 
NN=N(L)-1 
JAJ=JJJ 
JJJ=JJJ+NN 
DO  6 J=1»NN 
J J= JAJ+J 

6 B( JJ)=GAMMA( J*L) 

CALL  ARRaY(2»M»M»NAXDIM»NAYDIM»A1*A1) 

CALL  SIMO(Al*U*M*KS) 

IF  (KS.NE.O)  PRINT  100 
100  FORMAT ( 1H0 * 18HSYSTEM  IS  SINGULAR) 

JJJ=0 

DO  7 L = 1 » NO 
NN=N ( L ) -1 
JAJ=JJJ 
JJJ=JJJ+NN 
DO  7 J=1 *NN 
JJ= JAJ+J 

7 CH(J»L)=B(JJ) 

DO  8 L= 1 *N0 
NN=N ( L ) -1 
SCMIL i 

DO  8 1=1 *NN 

8 SCH(L)=SCH(L)+TSCD( I *L  *4 ) *CH( I *L ) 

IF (NX-1 ) 17  »9»lo 

9 DX=0 • 

GO  TO  11 

10  DX=(XMAX-XMIN)/FL0AT(NX-1> 

11  I F ( NY-1 ) 17*12*13 
1?.  DY=0. 

GO  TO  14 

18  DY=< YMAX-YMIN) /FLOAT INY-l ) 

14  DO  16  11=1 »NX 
XJ=XMIN+FL0AT ( I I-1)*0X 
DO  16  J J=1 *NY 
YJ=YMIN+FL0AT < JJ-1 ) *DY 
A 1 t I I *JJ)=0. 

DO  16  L I =1 *N0 
NN=N ( L I )-l 
DO  16  1=1 »NN 
X1  = XJ-X ( I *L I ) 

X2=XJ-X( 1+1 »LI ) 

Y1=YJ-Y( I *L I ) 

Y2=YJ~Y(  I + l.LI  ) 

Rl=Xl*Xl+Yl*Yi 

R2=X2*X2+Y2*Y2 

IF((R1.EQ.O.).OR.(R2.EQ.O.))  GO  TO  15 

YT  = Y( I *LI )*X2-Y( I+l ,LI )*Xl+Yj*(X( I + 1 * L I ) -X ( I * U I > ) 

XT=X1*X2+Y1#Y2 
TETA=ATAN2( YT*XT) 

Sl=TSCD< I »LI »4)*( l.-0.5#AL0G(R2/RR) ) +0. 5*AL0G < R2/R1 ) * ( X1*TSCD U * L I 
1 » 3 )+Yl*T SCO  ( I *L  I * 2 ) ) + TETA*Ul*TSCD(  I »L  I *2  )-Yl*TSCD<  I * L I * 3 > ) 

GO  TO  16 

15  Sl  = TSCD( I »LI »4)*( l.-u.5*AL0G( tRl+R2>/KK>  > 

16  Ai< 1 1 >JJ)=A1(II »JJ)+S1*CH( I * L I ) 

17  RETURN 
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SUBROUTINE  ARRAY  (MODE»I *J*N»M»S*D) 
DIMENSION  S 1 1 > »D( 1 ) 

N I =N-I 

IF (MCDE-1 ) 10U  * 100  » 120 
100  IJ=I*J+1 
NM=N*J+1 
DO  110  K=1»J 
NM=NM-N1 
DO  110  L=1 » I 
IJ=IJ“1 
NM=NM-1 

110  DtNM)=S( IJ) 

GO  TO  140 
120  IJ=0 
NM=0 

DO  130  K=1 * J 
DO  125  L=1 * I 
I.I  = IJ+1 
NM=NM+1 

125  S( IJ)=D(NM) 

130  NM=NM+NI 
140  RETURN 
END 

SUBROUTINE  SIMQ<A*B*N*KS) 

DIMENSION  AC  1) >B( 1) 

T0L=0.0 
KS=0 
J J=-N 

DO  65  J-1»N 
JY= J+l 
JJ=J J+N+l 
BIGA=0. 

I T =JJ-J 
DO  30  I = J >N 
I I T+ I 

IF(ABS(0IGA)-ABS(A( IJ) ) ) 20*30*30 
20  B1GA=A( IJ) 

IMAX=I 
30  CONTINUE 

IF(ABSCBIGA)-TOL)  35*35*40 
35  KS--1 
RETURN 

40  1 1= J+N* ( J-2 ) 

IT=IMAX-J 
DO  50  K=J*N 
11=1 1+N 
I 2= 1 1+1 T 
SAVE=A(  ID 
AC  I1)=A( 12) 

A ( 1 2 ) =SAVE 
50  AC  U)=A<  ID/BIGA 
SAVE=B« IMAX) 

B( IMAX ) =B  C J ) 

B ( J ) =SAVE/B I GA 
IFCJ-N)  55*70.55 
55  IQS=N*  < J~1 ) 

DO  65  I X=JY  *N 
IXJ= IQS+IX 
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Appendix  P 

IMPEDANCE  APPROXIMATION  FOR  LOW-LOSS  TRANSMISSION  LINES 


To  show  how  the  characteristic  impedance  of  a low-loss  uniform  transmission  line 
depends  on  the  loss  coefficients,  consider  the  following.  The  characteristic  impedance  of 
such  a line,  represented  in  terms  of  the  so-called  primary  constants  ??,  L,  C,  and  G,  is 


7 -|/-R  + juL 

Zq  G +;'coC 


(FI) 


or 


Zq  ~ (zo  )ii 


tl-jJL 

t oL 

i -/  -4-  ’ 

OiC 


(F2) 


where 


(zq)ll  =]/“ 


(F3) 


(ZQ)LL  is  the  characteristic  impedance  of  the  lossless  line  with  primary  constants  L and 
C.  Equation  (F2)  can  be  approximated  by 


Zq  ^ (Z0)LL 


/i  j sA 

2gj  L 

1-j  JL)' 

2 wC/ 


(F4) 


provided  R < wL  and  G < c oC.  By  simple  algebraic  manipulation  ZQ  in  Eq.  (F4)  can  be 
written  as 


where 


(Zq)ll 


l' 


1_;  Hl 


1 -j 


ad 
$LL , 


R 


“c~‘2  r I 


ft- 


(F5) 
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where 


To  see  the  effect  of  this  representation,  consider  Eq.  (5)  from  the  main  text: 

Pf  = [Re  Z0) 

From  Eq.  (F5),  Re  Z0  can  be  expressed  as 

Re  Zq  « {Zq)h 

To  assess  the  significance  of  the  bracketed  factor,  consider  the  following  example  using 
raicrostrip  parameters.  Let  ac  = 0.215,  ad  = 0.132,  and  f}LL  - 52.36,  corresponding  to 
an  operating  frequency  of  1.0  GHz.  For  this  case  the  bracketed  factor  in  Eq.  (F10)  has 
a value  of  1.000004  ~ 1.00.  At  f - 100  GHz  the  bracketed  factor  is  again  approximately 
1.00.  Hence  the  approximation 


tvb’2 

\Z0\2 


»-20z 


1 + 


acad 

0LL2 


( ad\ 
\0LL  ) 


lF9) 


(FlO) 


i?c  Z0  - (Z.0)LL 

is  quite  adequate  for  most  design  purposes.  Similarly 

\Zq\  ~ (ZQ)LL. 


(Fll) 


(F12) 
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Appendix  G 

DERIVATION  OF  EQ.  (23) 


The  derivation  of  Eq.  (23)  given  here  is  similar  to  that  given  in  Ref.  11.  Consider  a 
transmission  line  composed  of  Nc  conductors  and  ND  different  homogeneous  dielectric 
regions.  The  total  time-averaged  energy  stored  the  electric  field  is  given  by 

WJ.  - (1/2)  f ~ dr,  (Gl) 


where  r is  the  volume  obtained  by  taking  the  region  enclosed  by  two  cross  sections  per- 
pendicular to  the  axis  of  the  tranmsission  line,  the  cross  sections  being  separated  by  a unit 
length.  The  quantities  e and  D in  Eq.  (Gl)  are  functions  of  position. 

Consider  a slight  perturbation  of  the  dielectric  constant  5e.  Then 


= fSSS.  dT. 

J « •;  2s2 

T T 

Equation  (G2)  can  be  rearranged  using  the  equations 


l> 

1 

II 

LU 

(G3) 

VD  =p, 

(G4) 

V-(06  D)  = - E • 5 D + 0yS D. 

(G5) 

Hence  Eq.  (G2)  can  be  expressed  as 


8We 


r 


dr  - 1/2  /EHe  dr. 

T 


(G6) 


By  use  of  Gauss’s  theorem 
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c r c 

Jvi<f>8D)  dV  = ^ <p:  <f  SD-n  dl  » - s <t>i  t>Qi> 


i=i 


where  C,-  is  the  contour  encircling  the  ith  conductor  and  n is  a unit  vector  normal  to  and 
directed  inward  toward  the  conductor.  For  zero  space  charge,  5 p = 0.  Holding  the  po- 
tentials 4>i  constant  and  letting  6e  represent  a perturbation  of  only  , which  is  the  dielectric 
constant  of  the  hth  homogeneous  dielectric  region,  the  following  results: 


8We  = 


= <*,•  SQ,-  - (1/2)  f £2  8ek  dr, 


where  7/j  is  the  volume  enclosed  by  the  kth  homogeneous  dielectric  region.  Equation  (G8) 
can  be  reexpressed  as 


5eh  — 
2,  Wek, 


where 


Wek=(H2)  J ekE*dr. 


(G10) 


The  summation  2 0,-  6 Q{  can  be  rewritten  in  terms  of  a coefficient  of  capacitance  as 
i=l 


Nc 


2^  <PiSQi  = Msscti  = 28We 


i- 1 J-l 


(Gil) 


By  use  of  Eqs.  (G9)  with  (Gil)  the  following  result  is  obtained: 


dek  ek 


, k = 1,  2,  . . , Nrf  • 


(G12) 
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